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in Liquid Propellants 


offers outstanding advantages ,-| 
to designers of rocket motors ¢@ 


HIGH SPECIFIC IMPULSE: Nitrogen Tetroxide exceeds many other well- 


known oxidizers in pounds of thrust developed per pound 
of fuel consumed per second. 


EASY T0 HANDLE: Nitrogen Tetroxide may be shipped, piped and - J 
stored in ordinary carbon steel equipment. It possesses - i 
high chemical stability, high density, low freezing point, *. bs 
and a reasonably low vapor pressure. 


Nitrogen Tetroxide is available at tow cost in 125-pound I.C.C. approved 
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lands with ease on aircraft carriers 


“Straight up’’ bat-wing interceptor 


Douglas F4D Skyray 


Problem: to find a Navy interceptor to 
operate from aircraft carriers and meet 
the threat of today’s fast jet bombers. 
Answer: the Douglas F4D Skyray. 
With its mighty power plant and 
radical swept-back wings, Skyray zooms 
into action at blazing speed. Minutes 


after radar warning, it’s off the deck and 
on station—ready to intercept approach- 
ing aircraft with a lethal load of bullets 
and rockets. Yet for all its power and 
speed, this agile interceptor lands at low 
speeds—is perfectly adapted to require- 
ments of present carriers. 


Performance of the U. S. Navy’s F‘D 
Skyray is another example of Douglas 
leadership in aviation. Developing both 
military and civilian planes that can be 
produced in quantity—to fly faster and 
farther with a bigger payload—is the basic 
rule of Douglas design. 
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Servo-Stahilization of Low-Frequency Oscillations in a 
Liquid Bipropellant Rocket Motor 


FRANK E. MARBLE! and DALE W. COX, JR. 


Daniel and Florence Guggenheim Jet Propulsion Center, 


The recent work of H. S. Tsien concerning the servo- 
stabilization of rocket motors is extended to the liquid bi- 
propellant rocket motor. It is shown that by use of a 
feedback system containing a device to sense the combus- 
tion chamber pressure, a suitably designed amplifier, and 
aservomechanism which governs the propellant flow, the 
low-frequency oscillations which occur in the rocket con- 
figuration may be stabilized for any value of combustion 
time lag. A method is given for determining a transfer 
function of the feedback loop which will assure stable 
operation. The technique of the Satche diagram is em- 
ployed in stability analysis. 


Introduction 


Spee low-frequency instability of liquid propellant rocket 
motors has been attributed by Gunder and Friant (1),* 
Yachter (2), Summerfield (3), and Croeco (4) to the coupled 
oscillations of rocket chamber pressure and of flow rates of 
propellant in the feed lines. The instability arises from the 
fact that propellant combustion lags injection, due to the heat 
transfer required for propellant evaporation, time required 
for chemical reaction, ete., thereby introducing a phase dif- 
ference between the oscillations of the two systems. The 
most general of these analyses is that of Croeco which treats 
both monopropellant and bipropellant systems and gives some 
physical basis for the introduction of a time lag. This latter 
aspect has since been generalized by Tsien (5). 

Various possibilities have been discussed for removing the 
low-frequency instability through suitable choice of rocket 
configuration parameters, such as the propellant pressure 
drop across the injector nozzle. In some cases these changes 
may be undesirable or impossible; for example, in a rocket 
required to have a variable thrust, it may prove difficult to 
stabilize for all conditions of operation. Tsien (5) has in- 
vestigated the more general concept of removing combustion 
chamber instability by use of a feedback servo control which 
senses the oscillating pressure in the combustion chamber and 
changes the propellant flow by the proper amount and in the 
proper phase to damp the oscillation. In the same paper the 
feasibility of such a system is demonstrated for the monopro- 
pellant rocket. For the purpose of this demonstration, Tsien 
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California Institute of Technology, P. 1 > Calif. 


introduced a graphical method of stability analysis, based on 
a suggestion by Satche (6), which is the appropriate technique 
for analyzing systems with time lag in the same sense that the 
familiar Nyquist diagram is appropriate for systems without 
time lag. Using the Satche diagram, Tsien was able to place 
quite general restrictions on the feedback loop which would as- 
sure stable rocket operation. 

The present paper aims to extend these considerations to 
the bipropellant rocket and to find criteria on the feedback 
loop to provide stable operation. The authors are deeply 
grateful to Professor H. S. Tsien for initiating this problem 
and for his active interest during its progress, 


Dynamics of the Bipropellant Rocket 


The dynamics of the bipropellant rocket chamber, Fig. 1, 
are more complicated than those of the monopropellant rocket, 
due to the additional degree of freedom associated with the 
second propellant line and to the periodic changes of combus- 
tion gas temperature caused by variations in mixture ratio. 
Because only low-frequency oscillations are under considera- 
tion, the pressure at any time will be assumed uniform over 
the combustion chamber. 

If the dependent variables of the system are taken to be the 
chamber pressure, the oxidizer flow rate, and the fuel flow rate, 
three equations may be written relating these quantities to each 
other and to the physical configuration. They are the con- 
tinuity relation for the combustion chamber, and two condi- 
tions of dynamic equilibrium for the fluid in the oxidizer and 
fuel lines. 

Following Crocco’s notation (4), let M, be the mass of 
burned gas in the combustion chamber, 7, the mass rate of 
burning, and 7, the mass rate of ejection from the discharge 


nozzle. Then the relation d/dt(M,) = m, — m, expresses 
2 
| 


CHAMBER 


FUEL PUMP 
FIG. | SCHEMATIC DIAGRAM OF BIPROPELLANT ROCKET AND PRO- 
PELLANT FEED SYSTEM 
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pane 
¢ 
= 


the conservation of mass. Suppose the state of the dynamic 
ge’ is always ¢ close to an equilibrium value defined by 
9, ™, = m, = m, then the residence time which, on the 
average, the burned gas spends in the chamber, is #, 
M,/m. Anticipating the use of perturbation analysis, the 
fractional variation of burning rate and discharge rate may 


— m m, 
be introduced, = and pe = together 
m m 


with a dimensionless time z = t/d, to write the mass conser- 
vation as 


d M, 


The propellant burning rate m, differs from the total pro- 
pellant injection rate m; because of the combustion time lag. 
The time 7 required for the transformation of propellant into 
products of combustion may be represented by means of the 
integral 


f(p, T,) dt’ = C = const ............ [2] 


where f(p, T,) is a function depending upon the mechanics of 
the process. The chief quantities upon which this function 
depends are p, the uniform chamber pressure, and 7, the 
ambient gas temperature in the neighborhood of the injector. 
The meaning of this integral may be clarified somewhat by 
thinking of f to be, for example, the rate of heat transfer to the 
liquid propellant, a certain total quantity C of heat being re- 
quired before the transformation is complete. The princi- 
pal quantities associated with the heat transfer process are 
the ambient temperature and pressure so that the time 7 
required for transformation will be large or small depending 
upon whether the heat transfer rate is small or large. 

Now if the time lag were independent of time, the frac- 
tional] variation of burning rate u, would be identical with the 
fractional variation of injection rate measured at the time 7 
earlier. Denoting the dimensionless time lag 7/3, = 4, 
constant time lag would imply u,(z) = yi(z — 6). However, 
if the time lag is increasing during the process, the burning 
rate is depressed below the injection rate in proportion to the 
product of the propellant flow rate and the rate of change of 
time lag. Therefore, using dimensionless quantities and neg- 
lecting a term of second order in the small variations 


w(z — 6) =0.............. [3] 
The dependence of the time lag upon time is, in reality, 
governed by the integral, Equation [2], and differentiating 
this with respect to time gives 


flp(t), — flo(t — 7), Tt — (1 ai): 


Since it is assumed that the chamber pressure and ambient 
temperature vary only slightly from their steady-state values 


p and T,, the values of f may be approximated, followin, 
Tsien (5), through a Taylor expansion about p = , T, wy of 


AD, Ty) + ) tp(t) — pl + 
oT, [T(t) T,] 


TA) 


f{p(t — 7), = f(p, T,) + 


oT, 
Upon sila these relations into Equation [4] and de~ 
leting second-order differences, the variation of time lag may 
be expressed 


For a monopropellant rocket the combustion gas temperature 
is nearly constant so that the second term on the right side of 
Equation [5] vanishes. Calling 0 log f(p, T,)/d log p = n, 
Tsien’s generalization of Crocco’s treatment is obtained. 
For the bipropellant rocket the combustion gas temperature 
may experience considerable variation, resulting from the 
variation of propellant mixture ratio, and consequently 
another parameter 0 log f(p, T,)/O log = m enters the 
problem. Knowledge of combustion lag mechanics is so 
meager, however, that rational estimate of its importance is 
difficult. Consequently, as in previous investigations, the 
dependence of combustion time lag upon gas temperature 
will be neglected. Then, denoting the fractional variation 
in chamber pressure as ¢ = (p — p)/p, the expression /or 
dimensionless time lag variation is 


ine —S) — ofe)]..... 6] 


The fractional variation in total propellant injection rate 4;, 
which occurs in Equation [3], is most conveniently expressed 
in terms of the individual fractional variations of oxidizer 
flow, uo = (to — mo)/mo, and of fuel flow, uy =| (my — my)/ ry, 
and the steady-state mixture ratio 7? = tio/my. Again fol- 
lowing Crocco in defining H =? — 1/2(R + 1), depending on 
the steady-state mixture ratio, it is an elementary calculation 
to show that, to the first order, 


mi = (1/2 + A) wo + (1/2 — A) wy 


Consequently, the propellant burning rate is just 


= (5 H) poz — 6) + H) 
n[g(z — 6) — ¢(2)].... 


In calculating the variation of exit nozzle flow, m., it is 
convenient to consider the nozzle flow as quasi-steady, that is, 
the flow corresponds to a time sequence of steady equilibrium 
flow patterns. Then the ratio of instantaneous to steady mass 


discharge rate is 
[8] 
m p T Ss 


where 7, is the gas temperature at the exhaust end of the 
combustion chamber. Tsien has shown (7) that the quasi- 
steady approximation is adequate for low-frequency instabil- 
ity under discussion here. 

For the monopropellant rocket, where 7’, is invariable, the 
mass discharge rate depends directly upon the combustion 
chamber pressure p, and the fractional variation of mass dis- 
charge rate is simply 4, = ¢. Where the variation of mix- 
ture in a bipropellant rocket causes appreciable variation of 
temperature, this too must, according to Equation [8], be 
accounted for in computing the rate of gas discharge. Sup- 
pose the gas temperature at the exhaust end is 7’, at the time 
t. Reealling that the gas temperature depends upon only 
the propellant mixture ratio, the temperature 7, can be re- 
lated to the injection rates of propellant which formed the gas 
volume. The generation of gas took place at a time #, earlier 
than the gas appeared at the nozzle end of the chamber, that 
is, at a time ¢ — 8). Furthermore, the time required between 
injection of the propellants and generation of the gas is the 
time lag r. Therefore the gas temperature at the exhaust end 
of the chamber is that corresponding to the mixture ratio at 
the time t — 3, — 7, or at the dimensionless time z — 1 — 6. 

Now if r is the instantaneous mixture ratio, the change of 
gas temperature from its equilibrium value is, approximately, 
(dT,/dr); - (r — Furthermore, the variation of mixture 
ratio from its equilibrium value follows directly from the defi- 
nition r = M/m, ~ 71 + wo — wy). Denoting the di- 


log T, 


dt dlogp 


Tt — +) — 5] 
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mensionless slope of the gas temperature curve, (7/7,)(dT,,/ 
dr);, by 2K, the effect of combustion gas temperature upon frac- 
tional variation of mass discharge is K [uo(z — 6 — 1) — wy (z — 
§—1)]. Thus the complete expression is 


ue(z) = oz) — — — 1) — — 1)) .. [9] 


Finally, the time variation of gas storage within the com- 
bustion chamber is shown by Crocco to be 


d(M.\ _ de 

2K} — — 1) — — 6 — 1)] — [wolz — 8) — — 8)]} 
. [10] 


It ix easy to trace the origin of the terms entering the right 
side of this expression. When the temperature is constant, 
the variation of mass in the chamber is directly proportional 
to the fractional variation of the uniform pressure, that is, to 
dg dz. When the temperature is not constant, a term must 
appear to account for the fact that the gas leaving the dis- 
charge, generated at z — 6 — 1, is of different temperature and 
density from that generated at the injector, at z — 6. Con- 
version of these temperature differences to the corresponding 
difierences in oxidizer and fuel rates leads to the second 
term on the right side of Equation [10]. 

Introducing Equations [7], [9], and [10] in the original 
expression for the mass conservation in the combustion 
chamber, Equation [1], gives this relation in terms of the de- 
sired dependent variables 


+ +(1 + ng(z — 6) + 
— 6 — 1) — 2Kyo(z — 6) (5+ — 


[ Kus 6 = 1) 0) + (5 #) =0 


Description of the system may be completed by relating the 
individual propellant flow rates to the mechanical properties 
of the feed system. Consider each propellant feed system to 
consist, Fig. 1, of a feeding mechanism (e.g., turbopump), 
a propellant supply line, an injector nozzle, and a variable 
capacity situated in the line so close to the injector nozzle 
that the inertia of fluid between the capacity and the injector 
nozzle may be neglected. Since the process is the same in 
oxidizer and fuel lines, notation differentiating the two lines 
will be postponed. 

Denote by subscript 1 the conditions in the propellant line 
just downstream of the pump, and by 2 the conditions just 
ahead of the injector. Following Tsien (5), the variation of 
m,, the pump discharge rate, from its equilibrium value may 
be written 


When a = 0, this relation indicates no response of mass flow 
to a change of line pressure, that is, a positive displacement 
system. Conversely, when a = ~, no change of line pressure 
may be effected by changing the flow rate, which is charac- 
teristic of the constant pressure feed or pressurized tank. In- 
termediate values of a correspond to particular types of pump, 
and if the response time of the turbopump system is very long 
with respect to the period of chamber pressure oscillation, the 
pump will operate at constant speed and —a@ corresponds to 
the slope of the pump characteristic (Fig. 2). 

The pressure difference p; — p along the propellant lines is 
due primarily to the liquid acceleration through the injector 
nozzle and the inertial reaction of the fluid column in the 
propellant line. The frictional loss will be neglected. If S 
is the area of the propellant jet issuing from the injector, and 
p the propellant density, the pressure difference p. — p across 

he injector is m?/2pS*. Furthermore, if / and A are the 
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length and cross-sectional area of the connecting line, the in- 
ertial pressure difference follows from the momentum equa- 


tion (p, — p)A = d/dt(lm,). The mass flow m through the 
injector may differ from the mass flow ni, in the connecting 
line because of the variable capacitance situated between the 
two. Since the propellant is considered incompressible, the 
difference of the two values of mass flow may be expressed 
m — m, = dC /dt, where C is proportional to the instantaneous 
volume of the capacitance. This capacitance has the dimen- 
sion of mass per unit time X time and is conveniently repre- 
sented in the dimensionless form k = C/mvd,. Then if the 
pressure drop across the propellant injector under steady 
operation, m?/2pS?, is denoted Ap, two dimensionless parame- 
ters may be introduced: the inertia of the connecting line 

= Im/2APAB,, and the injector pressure ratio P = p/2Ap. 
Finally eliminating variable conditions at the pump and at 
the capacitance, the dynamic equilibrium of a propellant line 
may be written 


An equation of this type holds independently for the oxidizer 
and fuel lines. 
Stability of the Bipropellant Rocket 


When the propellant lines are of fixed capacity, the system 
is described by three independent variables, ¢(z), uo(z), and 
uy(z) representing the fractional variations in chamber pres- 
sure, oxidizer, and fuel flow. The three differential equa- 
tions are, for the combustion chamber 


+ ng(z — 6) + 


[ Kus — — 1) — 2Kpfz — 8) + (5 — =0 


and for the oxidizer and fuel lines 
A= + E + >) wot 
dz a 


The sian or instability of the system, that is, whether the 
amplitude of particular oscillatory solutions will diminish or 


an 
3 
tone 
| 
CONSTANT SPEED 


| PUMP CHARACTEFYSTIC 


PUMP MASS FATE, ™, 
FIG. 2. APPROXIMATION TO CONSTANT SPEED PUMP CHARACTER- 
ISTIC 


65 


ire 
of 
n, 
re 
he 
1e 
SO 
q 
he 
re 
m 
or 
| 
d 
n 
- 
: 
— 
= 
| 


grow, is most easily discussed through introducing the Laplace 
transform, defined as 


for the variable g(z). Similarly defining M,(s) and M;,(s) as 
the Laplace transforms of yo(z) and u,(z), each of the differen- 
tial equations may be transformed, as shown in Appendix 1, 
into inhomogeneous algebraic equations in the complex va- 
riable s. The inhomogeneous portion depends upon the initial 
conditions of the problem and hence is associated with the 
solution of transient problems. The corresponding homo- 
geneous equations 


[s + 1 —n + ne—S6](s) + 
e—sb [ Ke-s (2x 3) ates) 


[ Ke-s (2x +H >) | M,(s) = 0 


Py®(s) + ( + >) + 1| Mis) = 0 
ao 
+ [ys + = + >) 1 =0...[17] 


describe the free oscillations which are of interest in the stabil- 
ity problem. For nontrivial solutions the determinant of the 
coefficients of Equations [17] must vanish, that is 


| s+(1—n) +ne-s 

| 
ao 


P; 0 


Equation [18] is a transcendental equation in the complex 
variable s, the roots of which determine the stability of rocket 
chamber pressure oscillations: If any root of G(s) = 0 
posseses a positive real part, the system is unstable. To 
check on the existence of such roots, it is possible to apply the 
familiar Nyquist criterion which traces the course of the com- 
plex function G’(s) as s describes a counterclockwise contour 
in the right half of the s-plane consisting of the imaginary axis 
and a large semicircle in the right half plane. Then the nega- 
tive change in argument of G’(s) is the difference in number of 
zeros and poles of G’(s) in the right half s-plane. Graphically 
this number appears as the number of complete revolutions 
made about the origin by the trace of G’(s)._ If the number of 
poles is known independently, the number of zeros in the right 
half s-plane is known and the stability of the system is deter- 
mined. 

The present problem is complicated by the time lag 4, 
which not only makes calculation of G’(s) laborious but ac- 
centuates the fact that the value of 6 is usually not known ac- 
curately. To improve this situation, Tsien (5) has intro- 
duced another technique based upon a suggestion of Satche 
(6), which separates the function G’(s) into two parts, the first 
containing e~** as a factor, the second independent of the 
time lag. In the present example the separation may be af- 
fected by noting that the Equation [18] may be expressed as 
the combination of two determinants 


0 0 
Jos +1 + 
ao 2 
+1 + 
0 1 1 
= +3) 


[ (2K +H+ 3) 


such that Equation [18] may be written 

D,(s) + = [21] 
Following Tsien in calling g;(s) = e~** and defining g2(s) = 
—(D,(s)/D2(s)), the characteristic equation becomes 

— go(s) =G(s) =0.............. [22] 


where the complex function G(s) is represented as the dif- 
ference of the two vectors g;(s) and g2(s). 

Assume for the moment that G(s) has no poles in the right 
half plane. Then if s describes its previous contour, tiie 
number of zeros which G(s) has in the right half plane is equal 
to the number of complete clockwise rotations made by a vec- 
tor having its head on the trace of g;(s) and its tail on tlie 
trace of go(s). Since all values of g;(s) lie on or within tlie 
unit circle, the vector G(s) can make a complete rotation if tlie 


=G'(s) = 0 


Is +1 +2 (P +5) [18] 
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| 
| 
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trace of go(s) either passes through the unit circle or encircles 
the origin. Conversely, if the trace of go(s) lies outside the unit 
circle and does not encircle the origin, the system is definitely 
stable. 

Actually, the function G(s) may possess poles in the right 
half plane, due to zeros of the determinant D.(s) which con- 
stitutes the denominator of g.(s). Thus the number of 
clockwise turns which g2(s) makes about the origin must be 
modified by adding the number of poles induced through the 
roots of D.(s). This number of poles is found by applying 
the Nyquist criterion to D.(s) directly. 

A few examples will serve to illustrate this technique. To 
discuss the stability of any bipropellant rocket, two sets of 
parameters must be known; H, A, and n which depend upon 
the propellants and their mixture ratio; and Po, P;, Jo, Js, a, a 
which are associated with the propellant pumps, feed lines, 
and injector. In particular, the parameters H, K, and proba- 
bly m are not independent of each other; H depends upon 
the propellant mixture ratio, and K upon the slope of the 
combustion temperature vs. mixture ratio curve at the steady 
mixture ratio under consideration. For the examples to fol- 
low, these values will be taken as given in Table 1 for various 
mixture ratios. 


TABLE | 
7 TAPP) aT, H K 
d7 
2.15 5065 0 0.2335 0 
2.50 5010 360 0.2125 0.09 
2.25 4880 640 0.1925 0.1475 
2.00 4706 850 0.167 0.189 


For the first examples, choose the mechanical properties of the 
system such that 
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FIG. 8 SATCHE DIAGRAM FOR BIPROPELLANT ROCKET, 7 = 2.5, 
n= 03: ag = a; = 1.0, Po = P; = 1.0, Jo = 2:0; J; => 1.5 


= af 1.0 
Po = P; 1.0 
Jo = 2.0; J; = 1.5 


and consider four cases corresponding to two values of the 
steady-state mixture ratio and two values of the coefficient 
n = d log f/d log p. 


n ? 
1 0.2 2.5 
2 0.2 2.0 
3 0.6 2.5 
4 0.6 2.0 


Then for calculating the functions D,(s) and D.(s) for values 
of son the imaginary axis, that is, for s = tw 


Di(iw) n) + 31 — ur! + 
+ Pa — n)] 
D(iw) = +3 = 
iw (2K + H) +i +5K fiw cos w + w sin w} 


The curves of go(iw) for n = 0.2, calculated using these rela- 
tions, are shown in Fig. 3 for 7 = 2.5, and in Fig. 4 for 7 = 
2.0. They are nearly identical despite the considerable dif- 
ference in mixture ratio; neither of them intersect the unit 
cirele. To close the diagram, take s = Re’ where 3 de- 
creases from 2/2 to —(x/2). For large values of R, go(s) 
behaves as = —(R/n)e®. Thus 
the go(s) curve closes by a large are in the left half of the 
Satche diagram, progressing clockwise. Since g.(s) neither 
crosses the unit circle nor encircles the origin outside the unit 
cirele, instability can occur only if a zero of D.2(s) exists which 
cancels a root of ge(s) in the right half s-plane. The presence 
of such a zero may be detected by applying the Nyquist cri- 
terion to D2(s). The Nyquist diagram is shown in Fig. 5 for 
n = 0.2 andr = 2.5. For large values of R, D.( Re‘) be- 
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FIG. 4 SATCHE DIAGRAM FOR BIPROPELLANT ROCKET, 7 = 
n = 0.2; ao = ay = 1.0, Po = Py = 1.0, Jo = 2.0, Jy = 
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FIG. 5 NYQUIST DIAGRAM FOR 


haves as nJoJ;R%e?*”, and therefore the Nyquist diagram 
closes with a large are of 27 radians starting from the upper 
branch in Fig. 5 and progressing clockwise to the corre- 
sponding lower branch (not shown). Since the origin is not 
enclosed, the system is certainly stable. The same conclusion 


may easily be reached for the lower mixture ratio, 7 = 2.0. 
n = 0.2. 
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FIG. 6 SATCHE DIAGRAM FOR UNSTABLE BIPROPELLANT ROCKET, 
7 = 2.5, n = 0.6; ao = ay = 1.0, Py = Py = 1.0, Jo = 2.0, 
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FIG. 7 SATCHE DIAGRAM FOR UNSTABLE BIPROPELLANT ROCKET, 
r= 2.0, So = 0.6; a - a = 1.0, Py = P; = 1.0, Jo = 2.0, 
Jy = 1.5 


The stability curves for the two corresponding examples 
with n = 0.6, shown in Fig. 6 for 7 = 2.5, and in Fig. 7 for 
? = 2.0, both intersect the unit circle. Therefore a rocket 
propellant combination having parameters of these values 
will exhibit a low-frequency instability for some range of time 
lag 6. The method for determining the critical time lag has 
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FIG. 8 SATCHE DIAGRAM FOR BIPROPELLANT ROCKET, 7 = 2.0, 
n = 0.2; ao = ay = 1.0, Po = 1.0, Py = 4.0, Jo = 4.0, Jy = 1.0 


been shown by Tsien (5) in his original discussion of this 
method of analysis. Again the curves for the two values of 
mixture ratio are very similar. As in the first two examples, 
the closing ares in the Satche diagram indicate no further in- 
stability of the system. 

Because the mechanical] properties of the two propellant 
systems are nearly identical in the foregoing examples, the 
change in mixture ratio associated with an oscillation in 
chamber pressure is very small and the behavior is much like 
that of a monopropellant rocket without line elasticity. 
Consequently the Satche diagrams are particularly simple. 
The influence of mixture ratio oscillation becomes important 
when the two lines have different inertial constants and injec- 
tion pressures. As an example, consider a rocket for which 


The equations for D,(iw) and D.(iw) may be written 


49 
¢ 


(49°47 
io ~ 15(2K + H)t + 
9g 


5} K(cos w — isin w) + 15K(iw cos w + w sin w) 


The Satche diagram for this example, Fig. 8, exhibits large 
loops. These loops originate in the trigonometric terms of 
D.»(iw) which, in turn, arise because the change in combustion 
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temperature associated with a change 
in mixture ratio does not affect the 
chamber discharge rate until the resi- 


dence time J, later. Consequently these 
terms are associated with a reduced 


time lag of unity. 
Physically the influence of these 
terms is easy to see. When a gas 


mass of low temperature, caused by a 
fluctuation in mixture ratio, flows to the 
chamber discharge nozzle, the mass flow 


rate increases (cf. Equation [8]), causing 
the chamber pressure to drop. This re- 
duction of chamber pressure causes a 


variation of mixture ratio which may 
leal to another negative fluctuation in 
ga~ temperature. Clearly, at certain 


frequencies, these oscillations will rein- 
force each other so that they will, as in 
the case of the present calculation, be 


slowly damped. (As pointed out by 
T-:en, the Satche diagram indicates pos- 
sille slowly damped oscillations when 
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the ge(iw) curve passes close to the unit 
circle.) Furthermore, it is clear that 
there wil) be an infinite number of fre- 
quencies at which such a reinforcement will take place. In 
the present example only two show up to any extent. The 
alplitude of an oscillation which originates in this manner 
depends to a considerable extent upon the response of the 
propellant feed system. If the change in pressure produces a 
Jarge change in propellant flow, the oscillations may be 
severe; if the response is poor, the oscillation will disappear 
quickly. Now at high frequencies the inertia of the propel- 
lant supply causes the amplitude of its response to decrease, 
aiid consequently this particular resonance disappears. Mathe- 
matically this may be seen in the fact that for large values of 
w the term nJoJ;w? governs the behavior of D.(iw), whereas 
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e1G. 10 DETAILS OF SATCHE DIAGRAM NEAR UNIT CIRCLE, 
> = 2.0, n = 0.2; a = ay = 1.0, Po = 1.0, Py = 4.0, Jn = 4.0, 
Jy = 1.0 
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FIG. 9 NYQUIST DIAGRAM ASSOCIATED WITH SATCHE DIAGRAM OF FIG. 8 


the trigonometric terms are at most proportional to w. 

The closing are of the go(s) curve is the same as those dis- 
cussed in previous examples. The Nyquist diagram of the 
denominator D.(s) is shown in Fig. 9, indicating no zeros of 
this determinant in the right half of the s-plane. The details 
of the Satche diagram near the unit circle, Fig. 10, show that 
the curve lies outside the unit circle and, as a consequence of 
this and the above considerations, the system is conclusively 
stable for all values of the combustion time lag. 

The influence of mixture ratio variation may be further in- 
creased by considering a propellant combination for which the 
curve of combustion temperature vs. mixture ratio has a 
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slope greater than that considered in the previous 
example. Let all values of the previous example 


be retained with the exception that K is arbitra- 
rily increased from 0.189 to 0.30. The resulting 
Satche diagram, Fig. 11, shows that the loops 


have now grown so large as to encircle the origin. 
By drawing the closing are and checking the 


course of the complete curve, it is simply found 
that the g.(s) curve encircles the origin twice in 
a counterclockwise direction. The Nyquist dia- 


gram of 9.(s), Fig. 12, shows a corresponding 
double encirclement but in the clockwise direction. 


Consequently, this cancellation indicates that no 
instability can arise from encirclement of the ori- 


gin in the Satche diagram. Furthermore, the 
details of the Satche diagram in the neighborhood 
of the unit circle, Fig. 13, show no intersection 


| | 40 


and, thus, the system is definitely stable against 
all values of the combustion time lag. 
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Servo-Stabilization of the 
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Bipropellant Rocket 


Where possible unstable operation of a rocket is 
indicated, stability may often be assured by modifying the me- 
chanical design of the propellant supply system or injector, or in 
the extreme, by change in the chemical propellants. Alterna- 
tively the stability may be assured by introducing a feedback 
loop connecting the combustion chamber and the variable 
capacitances in each propellant line. In the simplest case the 
loop will consist of a pressure measuring device attached to 
the rocket chamber, an amplifier, and an appropriate servo- 
mechanism to actuate the capacitances in the propellant lines. 
For the bipropellant rocket, Fig. 14, two separate circuits may 
be employed, one for each propellant line. However, in the 
interest of simplicity it may be preferable (Fig. 15) to control 
only a single propellant line, or to control both through a 
single amplifier, building any required differences into the 
servomechanisms or into the capacitances.‘ 

In general, then, the capacitances xp and x; will be made to 
vary with time in a manner prescribed according to the varia- 
tion in chamber pressure. Employing dimensionless forms, 
the correspondences between the capacitances and the cham- 
ber pressure may be written symbolically 


where F,(d/dz) and F,(d/dz) are linear differential-integral 
operators. These relations simply indicate that the two 
capacitances are related to the chamber pressure through 
linear differential equations with constant coefficients, the 
forms of which depend upon the pressure pickup, the ampli- 
fiers, and the servomechanisms. 

The problem is then to design the feedback circuits, that is, 
to determine the forms of Fo(d/dz) and F,(d/dz) such that 
stable rocket operation is achieved. This resolves, practi- 
cally, to the design of amplifiers to work with given servo- 
mechanisms and with a given pressure pickup so as to pro- 
duce the required forms of Fy and F;.. The proper forms of 
these functions can be determined only through stability 
analysis of the complete set of five differential equations, the 
original set consisting of Equation [11], and two of the type of 
Equation [13] augmented by the two equations arising from 


‘ It has come to the attention of the authors that Y. C. Lee, 
M. R. Gore, and C. C. Ross, of the Aerojet Engineering Corpo- 
ration, have considered a somewhat different arrangement in their 
paper, “Stability and Control of Liquid Propellant Rocket Sys- 
tems” (see pp. 75-81 in this issue of the JouRNAL). 
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FIG. 12 NYQUIST DIAGRAM ASSOCIATED WITH SATCHE DIAGRAM OF FIG, 11 
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FIG. 13 DETAILS OF SATCHE DIAGRAM NEAR UNIT CIRCLE, 
K = 0.3,7 = 2.0, n = 0.2; av = ay = 1.0, Py = 1.0, Py = 4.0, 


Jo = 4.0, Jy = 1.0 


the feedback loops, Equations [23] and [24]. Two additional 
dependent variables xo(z) and «;(z) have been introduced by 
allowing the capacitances to vary. The Laplace transforms oi 
these equations are given in Appendix 1, and retaining only 
the homogeneous portion for stability analysis, the set of 
algebraic equations is 


[s + 1 — n + + 
eal Ke-s (2x >) | 
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wl ere Ko(s) and K;(s) are the Laplace transforms of the varia- 
bl- capacitances; Fo(s) and F;,(s) are the transfer functions 
of the feedback loops. The stability of the system is now de- 
peident upon the roots of the fifth-order determinant con- FIG. 14. SCHEMATIC DIAGRAM OF BIPROPELLANT ROCKET INCORPO- 
st ucted from the coefficients of Equations [25] RATING FEEDBACK SYSTEM 
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| F;(s) 0 0 0 [26 ] 
This determinant is most conveniently expressed through ex- F,(s), appear as multiplicative factors; each of the determi- 
pansion by second-order minors, employing the last two rows nants depends upon only the original rocket configuration and 
for the expansion. Defining two new third-order determi- is independent of the feedback circuit. 
nants The new form of the function g»(s) is 
D(s) = 
! 
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The characteristic equation obtained by setting the stability 
determinant (Equation [26]) to zero is 


Di(s) + e—55[Do(s) + Fo(s)Do(s) + Fy(s)D;(s)] = 0. . [29] 
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FIG. 15 TWO FEEDBACK CIRCUITS EMPLOYING ONLY A SINGLE 
The two transfer functions of the feedback loops, F'o(s) and AMPLIFIER AND PICKUP 
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Stable rocket operation is to be assured by choosing realiza- 
ble forms of the feedback transfer functions Fo(s) and F;(s) 
so that the Satche diagram satisfies the criteria previously 
discussed. One convenient procedure is to choose an ap- 
propriate Satche diagram and solve for the feedback transfer 
functions. For the present system the result is usually not a 
realizable or convenient feedback circuit. 

When the instability arises from intersection of the g2(iw) 
curve with the unit circle, as it appears to do in practical cases 
of low-frequency oscillation, a proper choice of F(s) and 
F;(s) may be effected quite simply. The intersection, when 
it exists, occurs for value of w considerably less than unity; 
if this were not true, the oscillations would be of rather high 
frequency considering residence times of the usual magnitude. 
Hence the modification of the g2(iw) function may be restricted 
to low values of w, this is to say, the gain of the amplifier in 
the feedback system should drop to a low value as the fre- 
quency w increases. The appropriate analytic procedure is 
then (a) to find the approximate behavior of Fy(s) and F;(s) 
for s = iw<<i such that go(iw) does not intersect the unit 
circle, and (6) to determine the functions Fy(s) and F;(s) so 
that they have the prescribed behavior for small values of w 
and vanish for large w, thus preserving the original Satche 
diagram in this range. 

This program may be carried out in the following manner. 
Each of the third-order determinants involved in Equation 
[30] is regular at s = 0, and hence may be approximated by a 
power series 


[30] 


D,(s) = + + dis? + dys? 
= de + dos + dos? +. 

Do(8) = do + + dos? + dos? + dost + .. 

Ds) = de + + dys? + dys? + + ... [31] 


The coefficients d;4 depend only upon the parameters of the 
original (unstabilized) rocket and hence are known numerical 
values for any particular case. The expressions for these 
coefficients are given in Appendix 2. In particular, it ap- 
pears that dg = d;® = 0, a result which is clear from in- 
spection of Do(s) and D,;(s) as given by Equations [27] and 
[28]. The transfer functions for the feedback loops Jikewise 
possess power-series expansions about s = 0, although they 
need not be regular at that point. From the fact that the 
determinants Do(s) and D,;(s) behave as s for small values of 
s, the expansions of the feedback transfer functions must be- 
have as 1/s in order that the value of g.(0) be changed from 
that of the original system. That is, the transfer function 


must correspond to a simple integrating circuit. Thus it is 
appropriate to choose 
. . [82] 


= + + fps? +... 


where the fo and the f;% are unknown constants. Clearly, 
upon substituting the approximations of Equations [31] and 
[32] into Equation [30], a power-series representation for 
g2(s) is obtained where the unknown coefficients fy, f; 
appear in the coefficients of the expansion. 

From the viewpoint of insuring against low-frequency in- 
stability, the Satche diagram must be such that the go(iw) 
curve lies outside the unit circle. For small w, which are in- 
volved near the unit circle, go(iw) may be approximated by the 
polynomial 


go(iw) = yo + yiliw) + yo(tw)? + ya(tw)* 


where the 7; are constants chosen so that this condition is satis- 
fied. Then there are two expressions for g2(s), valid for small s 
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g2(8) = yo + + 28? + (33) 


and that which followed from expansion of Equation [380]. 
Since they are both valid power-series expansions, they «re 
equal term by term. The result is four algebraic equations 
relating the constants fy‘ and f;‘) to the rocket parameters 
and to the chosen form of the Satche diagram near the unit 
circle. After some simplification these relations are 
d, 
Yo 
dy fo do fo + + + = 


yo Yo dy 
dy fo. dy fy + dy fy 4. + + 
d,; _ 4) 
dg + de vo | 
[86] 


+ +d = — 


dy fo 4 dy fy + do fy) + dy fo?) dy + 
+ 4+ 4+ = 


[37] 


Four more relations among the fy? and f;“ may be chosen 
arbitrarily so as to fix their values uniquely. Two possibili- 
ties have been mentioned previously: (a) the feedback loop 
operates only on a single capacitance so that either fy) = 0 
or f; = 0; (b) the two feedback loops are identical so th:t 
fo” = f,; j = —1,0,1,2. Ineach case only one pressure- 
sensing device and amplifier are required. 

With the values of fo, f,9 known, the behavior of the 
feedback circuit is determined for small s; the task is now to 
obtain general representations for F(s) and F;(s) which agree 
with the known expressions for small s and vanish for large s 
so as to leave the original Satche diagram unchanged for large 
s. If, for instance, it were decided to control the oxidizer 
flow and leave the fuel flow uncontrolled, f,4) = 0 and the 
fo’? are determined by successive application of Equations 
[34], [35], [36], and [37]. Then choose 


1 — cs ; 
It is clear by reference to Equation [32] that the singular be- 
havior as s — 0 is correct. Furthermore for large s, Fo(s) 
vanishes as 1/s* which is probably strong enough. If the 
constants a, 6, and ¢ can be found so that the expansion of 
Equation [38] yields the first few coefficients in agreement 
with the known values fo, fo“, and fo, then Equation [36] 
will constitute a transfer function of the feedback loop, satis- 
factory for stabilizing the rocket against low-frequency os- 
cillations regardless of the propellant time lag. 
The pret expansion of Equation [38] is 


(a? — bla — c)s? +... }... [89] 


which must be equal, term by term, to the first of Equations 
[32]. This condition determines the constants a, b, c directly 


as 
So? fo fo 
fo 
fo 
D 


where the fy‘? are known. Substitution into Equation [38] 
then supplies the transfer function of the feedback loop to 
stabilize the system. 
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To illustrate application of the stabilization procedure, con- 
sider the unstable rocket for which the Satche diagram is 
shown in Fig. 6. For this example, the parameters have 


values 
a = ary = 1.0 
Po = P; = 1.0 
Jo = 2.0 Jy = 1.5 
7 = 2.50 n = 0.6 


The relations of Appendix 2 allow computation of the co- 
effivients in the power expansions of the third-order deter- 


minants 
d,® = 2.50 d, = 6.25 di? = 3.01 
= 9.75 = 6.85 dy) = 6.16 
d,\®) = 9.95 d,® = 1.76 do? = 2.89 
= 3.00 dx = 0.02 dy = 0.01 


Now if, for small values of w, the Satche diagram of the sys- 
ten: may be moved to that indicated by the broken line in 
Fiy. 16, the rocket motor will be stable against low-frequency 
oscillations. The curve g2(iw) will lie as shown if the y; are 


chosen 
—2.0 = —1.0 
By application of Equations [34] and [37] 
= —1.41 fo? = —3.87 
fo = —2.95 fo? = 5.05 


which fixes the feedback loop for small values of w. If further 
the transfer function is chosen to be of the form of Equation 
[3s], the constants a, b, and c may be calculated from Equa- 


tions [40] 
a = —0.8 a = —1.29 
b= 1.08 b= 0.06 
ec = 1.29 c= 0.8 


the two sets arising through the quadratic equation determin- 
ing a. Choosing the first set (arbitrarily), the transfer func- 
tion for the feedback loop in the oxidizer line is 
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11G. 16 COMPARISON OF STABILIZED AND UNSTABLE SATCHE DIA- 
GRAMS, 7 = 2.5, n = 0.6; ao = ay = 1.0, Po = P; = 1.0, Jy = 2.0, 
Jy = 1.5 
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—1.41 (1 — 1.298) 
(1 + 0.88)(1 — 1.088?) 


= 


The Satche diagram for the stabilized motion is now easily 
computed through use of Equation [30] and the feedback 
transfer function found above. The result is shown in Fig. 16 
along with the original Satche diagram and the power-series 
approximation for small values of w. It is of interest to ob- 
serve from Fig. 16 that the stabilized curve of go(iw) and the 
corresponding series approximation diverge appreciably only 
as the value of w exceeds unity. Thus the character of the 
Satche diagram is determined with sufficient accuracy by a 
few terms of its power-series expansion. 


Concluding Remark 


The present analytical results together with the initial 
work of Tsien (5) demonstrate conclusively that the low- 
frequency oscillations of a liquid propellant rocket motor 
may, in all cases, be stabilized by use of a feedback servo 
loop. Through use of the Satche diagram and_ the 
associated methods which have been described, a projected 
rocket configuration may be analyzed for possible unstable 
operation and the requirements for a stabilizing feedback 
loop determined. The application of these methods lies 
within common engineering practice. The concept of feed- 
back stabilization removes restrictions on the propellant feed 
system and injector which would otherwise limit the freedom 
of the designer, with the consequence that he should be able 
to focus attention on improving combustion, reducing pump- 
ing loss, ete. This freedom is particularly important in the 
design of a rocket with variable thrust. 


APPENDIX 1 


Laplace Transform of the Rocket Stability Equations 


The equation describing the combustion chamber pressure 
oscillation is 


dg 
(1 + ng(z — 6) + 


[ Kutz (2k +H +5) - = 


while the equations of motion for the liquid in the oxidizer and 
fuel lines are 


J, +5 3) + 1 + Pog + 


dko 
Jo Ge +5) [42] 


Gi + (P+ 3) = 0... [43] 
Denoting the Laplace transforms of ¢(z), uo(z), «o(2), 
and «;(z) by ®(s), Mo(s), M;(s), Ko(s), and K;(s), respectively, 
the transformation of Equations [41], [42], and [43] is 
straightforward except for the terms involving a time lag. 
In general, the transform of f(z — 6) will be 


f(z — d)e—s2 dz 


where F(s) is the transform of f(z). As is well known, when 
the variable is retarded the transform is modified through 
multiplication by an exponential factor e~** and by the addi- 
tion of an integral extending over the period of the time lag. 


73 


Be 4 
= 
= 
| 


The equations then transform to 
[s +1 —n + + 
(0) + f o(é)e— dé + 
+6 —st dt — 


[ + >) | Mis) + + 
+ 2 (Po +5) Kuls) = + (0) + 


[vs + re 3) [45] 


[4s +1++ | Mj(s) + P)®(s) + 
[yt + 2 + 5) 8] = + Jy GE (0) + 
[4s += + >) | «(0)... . [46] 


Stability analysis requires only the left-hand sides of these 
equations. Transient behavior of the rocket system, how- 
ever, is determined by the complete equations where the right- 
hand sides, in usual problems, consist in the initial conditions 
of the dependent variables. For systems with time Jag, 
such as the iiquid propellant rocket, Equations [44], [45], 
and [46] show that in addition to initial conditions evaluated 
at z = 0, certain features of the behavior of the system must 
be known for a period preceding z = 0 equal in duration to 
the time lag. Physically this is clear, for due to the combus- 
tion time lag, the variation of chamber pressure ¢(z) cannot be 
described for z>0 unless the behavior of the propellant flow 
is known for the period —6 < z < 0. 


APPENDIX 2 


Power Series Representation of the Determinants 


Dy(s) = dh + dis + + dys! 


Since this determinant, Equation [19], is a polynomial, 


this expression is exact. 
= (1 —n) +1 ++ (P; 
2)] + 


+2 (>, + 
[1 + >) | 
= (1 +Ju[1+ +3)] + 
| 


Jy 


Do(s) = do + dos + dys? + +... 
Referring to Equation [20], the coefficients are 


+2 (P+ 5) | [i+ i (e+ >) + 
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= Joly 


1 


2 | 


d,\)) 


R 
dy 


dy 


dy = 


Jy (p, + a) +5 5 (Poly + + 


(K + H){ Poy — 


= nJoly + K {Poly — Pro} + 
K \p 1 1 
{PyJo — + 
K \p 1 1 1\7) 
= iP +35)]- [1 +5 (+3) 


Do(s) = dy + dys + dy‘ + dy s4 + 


eferring to Equation [27], the coefficients are 


dy 


+K [1 +2 +5)]- 


0 = +2 (P, + - >) | + 


DAs) = ds + + + dys? + ast +... 


Referring to Equation [28], the coefficients are 


dy 


d; 


J 


= (K +H Joly + 


0 (P+ 
ay 


(Continued on page 81) 
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Stability and Control of Liquid Propellant Rocket Systems 


Yy. C. LEE,’ M. R. GORE,’ and C. C. ROSS* 


Aerojet Engineering Corporation, Azusa, Calif. 


Within the limitations of assumptions made in the 
anslysis presented herein, it appears that the problem of 
stability and control of bipropellant liquid rockets may 
be ‘nvestigated with respect to system dynamics, combus- 
tion time lag, and interaction of these two. For the gas- 
pressurized, bipropellant, liquid rocket system, stability 
criteria relating the system parameters have been estab- 
lished. Stability conditions for a typical pump-fed rocket 
are discussed. By the utilization of a suitable feedback 
control, it is theoretically demonstrated that bipropellant 
liquid rocket systems may be stabilized not only for system 
dynamics but also for any combustion time lag. 


Nomenclature 


area 
= characteristic velocity 

frictional coefficient 

specific heat at constant pressure 

= specific heat at constant volume 
propellant line diameter 

= mass moment of inertia of turbine wheel 
characteristic length of thrust chamber 
propellant line length 

= equivalent propellant line length 

= pressure 

Laplace operator 

= gas constant of chamber gases 

= chamber gas temperature 

= time 

chamber volume 

= velocity 

w = weight rate of flow 

6 = turbine speed 

p = density 

7 = combustion time lag 

w = frequency 


Il 


~~ 


ll 


(Bar over a symbol denotes equilibrium conditions.) 


Subscripts 
= thrust chamber 
d = propellant pump discharge or propellant tank 
j = injector 
f = fuel 
ox = oxidizer 
t= chamber throat 


Introduction 


(A) INSTABILITY in liquid-propellant rocket systems 
may be defined as operation in which uncontrolled varia- 
tions in the system parameters occur. The rocket engineer is 
concerned with these variations because they may be of such 
magnitude as to cause poor performance of the propulsive 
unit or the vehicle which it propels. The rocket engineer 
would like to be able to design systems which are not only 
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stable but also are free from large-amplitude transient oscilla- 
tions. This paper suggests a method by which adequate 
stability and control of liquid propellant rocket systems may 
be achieved as part of a design procedure. The feasibility 
of such a method is demonstrated by the analysis presented 
herein. 

(B) To facilitate discussion and analysis of the problem 
of stability and control, two mechanisms of instability may be 
distinguished. The first may be defined as one that involves 
only the dynamic behavior of the rocket feed system and its 
components. A rocket system is composed of many physical 
components, each characterized by its own dynamic behavior. 
Pumps, turbine, thrust chamber, and hydraulic flow through 
propellant lines are incapable of instantaneous response to 
changes in operating conditions. They are all characterized 
by a finite delay time. Time delay gives rise to phase shift 
and amplitude variation between system output and input, 
which are a function of the frequency of excitation. A dis- 
turbance such as a pulsation in propellant flow or chamber 
pressure may be of sufficient amplitude to initiate sustained 
oscillations in an unstable system. 

(C) The second mechanism of instability is directly as- 
sociated with the combustion phenomena of propellants in 
the thrust chamber. This mechanism of instability is gener- 
ally referred to as combustion instability. There are two 
types of combustion instability. The first type involves pres- 
sure oscillations of frequencies low enough so that sound 
waves can traverse the thrust chamber several times during 
one oscillation. This type is generally referred to as ‘“chug- 
ging’’ and usually occurs in the frequency range from 40 to 
400 cps. The second type is associated with the motion of 
sound waves in the thrust chamber, and the observed fre- 
quencies of pressure oscillations are the characteristic fre- 
quencies of the various acoustic modes associated with a given 
chamber. This type is commonly called “screaming”’ in- 
stability and usually occurs in the frequency range from sev- 
eral hundred to several thousand eps. It is with the first 
type of instability, or chugging, that this paper is primarily 
concerned. 

(D) The forcing function, or the driving mechanism, for 
the pressure oscillation of both chugging and screaming in- 
stability is associated with the atomization, mixing, and reac- 
tion processes of the propellants in the combustion chamber. 
The combustion mechanism depends on the fact that there is 
a great deal of energy available in the unreacted propellant 
components in the thrust chamber, and that a finite time is 
required to release this energy; or, to put it differently, there 
is a delay between the time when liquid fuel and oxidizer are 
injected into the thrust chamber and the time when these 
propellants are converted into gaseous combustion products 
at the steady-state chamber pressure and temperature. In- 
stantaneous variations of any system parameter can change 
the instantaneous rate of energy release. If, for example, 
there is a pressure variation, and if the delay is of the proper 
magnitude, then the instantaneous increase or decrease of 
energy release due to the pressure change will be in the correct 
phase relationship with the resulting pressure oscillations, 
and these oscillations will increase in amplitude. Therefore, 
it appears that the rocket engineer is confronted with in- 
stability problems associated with the dynamic characteristics 
of the feed system and with the injection-combustion proc- 
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esses of the thrust chamber and the interaction of these two. 

(E) Increasing attention has been given recently to the 
problem of stability and control of liquid-propellant rockets. 
Particularly, the problem of combustion instability in liquid- 
propellant rockets has been exhaustively treated recently by 
Dr. M. Summerfield (1),4 Dr. L. Crocco (2), and Dr. H. S. 
Tsien (3). Basically, the key to the theory lies in the concept 
of a time lag which exists during the process of propellant 
injection and transformation into high-temperature gases. 
Summerfield assumes that the time lag is constant, and Crocco 
postulates that the time lag is dependent upon chamber pres- 
sure. The assumption of a constant time lag implies that 
the system cannot become unstable or self-excited if the pro- 
pellant flow is maintained constant. On the other hand, if 
the time lag is dependent upon chamber pressure, instability 
or self-excitation can result, even if the propellant flow is 
maintained constant. This represents the important and 
basic difference between Crocco’s and Summerfield’s work. 
Tsien, in his work, treated the same problem and attained 
the same results as Crocco by the use of an artifice of the 
Satche diagram (3). He further amplified his treatment of 
the problem by demonstrating the theoretical feasibility of 
achieving stabilization by means of feedback control in the 
propellant circuit. 

(F) The first part of the present paper essentially extends 
Summerfield’s work from a monopropellant to a bipropellant 
system for a simple gas-pressurized liquid-propellant rocket. 
The analysis of this system lends itself well to the servo- 
mechanism approach and may be treated most conveniently 
by the use of operational methods and the Nyquist criteria. 
Results of the analysis are given as stability boundaries re- 
lating the critical combustion time lag and critical frequency 
with the system parameters. The second part of this paper 
deals with a stability analysis of a rocket system, having fuel 
and oxidizer turbopumping plants. Both calibrated and 
feedback control systems are analyzed, and the effects of both 
system dynamics and combustion lag time are considered. 


2 Gas-Pressurized Bipropellant System 


(A) The gas-pressurized bipropellant rocket system is a 
simple system in which constant propellant-tank pressures 
are maintained by means of a pressure regulator, and pro- 
pellants enter the thrust chamber through lines connecting 
the propellant tanks to the chamber. Fig. 1 is a schematic 
diagram of such a system. 


OXIDIZER 
TANK 


REGULATOR 


THRUST 
CHAMBER 


TANK PRESSURE 
ORIFICE 


MIXTURE RATIO 
ORIFICE 


FIG. 1 BIPROPELLANT GAS-PRESSURIZED LIQUID ROCKET SYSTEM 
( SCHEMATIC) 


(B) To facilitate the analysis of this system the following 
assumptions have been made: (a) A constant lag time, r, 
exists between the time of propellant injection and time of 
equilibrium chamber pressure. (b) Propellant compressi- 
bility is negligible. (c) Propellant line elasticity is negligible. 
(d) System capacitance, such as trapped air in the Jines or 
injector manifold, is negligible. 

(C) The basic equations governing a monopropellant 
system have been developed by Summerfield (1). Those fora 
bipropellant system are presented as follows: 


‘ Numbers in parentheses refer to References on page 81. 
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(a) The flow of oxidizer and fuel to the thrust chamber at 
a time ‘—r must be equal to the rate at which the propellants 
are exhausted from the thrust chamber plus the rate at which 


they are stored in the chamber. Mathematically 
dP, 
Woz(t—7) + = 1} 


From energy Pan the external work done on t!ie 
propellant per unit time must be equal to the change in kinetic 
energy plus losses. For the fuel circuit 


— PAs Vi, = Va? + 


Vig Aun dt 


4 
since 

AgVa = A;V; = AsV; 
The above*equations may be simplified to 


1/l dw, 
+ Ay + Aj, 


Similarly, for the oxidizer circuit 


€osWox* 1 Lies ) dwWox 
ZporGA jor? 


A 


A jes Ais ’ 

(b) Inspection of Equations [1], [2], and [3] indicates that 
the system governed by them is analogous to a simple fee- 
back servomechanism. In fact, the system may be repre- 
sented by a block diagram, shown as Fig. 2. The inputs to 
the system are the tank pressures, Pa, and Pa,,. The output 
of the system is chamber pressure, P,, which is compared wit!) 
the input pressures. The error or difference between input 
and chamber pressures governs the propellant flow to the 
thrust chamber as described by Equations [2] and [3}. 
Since Equations [2] and [3] are nonlinear, small variations 
about an equilibrium point may be investigated. This is 
accomplished by linearization. The linearized system of 

equations results in 


Pa. — Po = 


where 


€f 


A Ve 


2 = 
Py Pe = (Pa ~ Pe) 
1 ic 
Pe = (Paw — wor + Worx 
oz joz 


These equations may be aieliste in operator form as the 
transfer functions appearing in Fig. 2. 


E 

3 T,P+l 
Wox* g- P; 

2 


FIG. 2. BIPROPELLANT GAS-PRESSURIZED LIQUID ROCKET SYSTEM 
(BLOCK DIAGRAM) 
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G = P. Kie — 
; Woz + Ws (Tip + 1) 
G, = = Pe [5] 
G; Wox K; 


Pas — Po” (Tp +1) 
where 7), Ts, ete., are time constants, and 
She kL* 


. 
K Ag 
k =° 
Ce 
2(Py — P.) 
ly 
k Wor 
2(Pas — Pe) 
— 29(Pacz — 
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From Equations [5], the system equation is 


GiG2P ay + GiGsP doz 
[1 + + G;)] 


The characteristic equation, F(p), is 


Ko(Tsp + 

+ (ip + + IT +1) [6] 
where 
Ky = Ki(K2 + K;) = + Dea) + Boa) ai + 

7, Ke 3 + Qor 


(ec) 1, T2, T3, and 7; are time constants of the system. 
7, represents the time constant of the thrust chamber, 7, 
represents the time constant of the fuel circuit, and 7’; the 
time constant of the oxidizer circuit. 7, is a time constant 
associated with the coupling of the fuel and oxidizer circuits. 
Ko is the static gain of the system. 

(d) The fundamental criterion for the stability of this 
system is that F(p) have no zeros in the right half of the com- 
plex p-plane; and system stability may be investigated by 
applying Nyquist’s criteria to Re~*® where 


Ko(T.p + — 7? 
(Tip + 1I)(T2p + 1)(Tsp + 1) 


Re-t¢ 


The system is unstable if R > 1 when the phase angle, ¢ = 
180 deg. The system is neutrally stable if the amplitude, R = 
| when ¢ = 180 deg. It is to be noted that for a monopro- 
pellant system or for equal oxidizer and fuel tank pressures, 
Ko, the static gain of the system, reduces to 

Ko = 
2(Pa — 
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and Equation [7] reduces to 

(Tip + 1)(Tep + 1) 
This is a direct result of the fact that 7, the coupling time 
constant, is equal to T, and 73, the fuel and oxidizer circuit 
time constants, respectively. 

(e) Based upon the criteria of stability, the following cases 
may be considered. 


= 


Case I: 7=0 


For r = 0, the stability criteria are due entirely to system 
dynamics. The system is always stable for all values of Ko, 
Ts, and because: 

First: If Ky <1, this system is stable for all values of 7}, 
T», T:, and Since T, = + ‘Ke > Ks, is 
always smaller than either T, or 73. Therefore as p = tw 
takes on all values of w, the amplitude R is always smaller 
than unity. A phase shift of 180 deg is possible only as 
w—> ©andR—0. 

Second: If Ky >1, the amplitude R can be greater than 
unity at some frequencies. But a phase shift of 180 deg is 
possible only asw + © and R +0. 


Case IT: 7 ¥0 


For this case, the stability criteria due to the interaction of 
combustion time lag and system dynamics may be investi- 
gated. 

First: If Ko <1, this #ystem is stable for all values of 7), _ 
T», T3, and T;. Since 7 is always smaller than either 7, or 
T3, and e~*” only contributes to phase shift, the amplitude 
R can never be greater than unity at any frequency. It is 
to be noted that this result is the same as that of a monopro- 
pellant system, investigated by Summerfield (1). 

Second: If Ky > 1, there always exists some 7 regardless 
of values of 7,, Ts, T3, and 7 to cause a phase shift of 180 
deg when R = 1. 

(f) For the general case, it is desirable to establish the 
relationship which must exist among the system parameters 
in order that stable operation may be possible. The effect of 
combustion lag time, 7, upon system stability may be best 
investigated by a polar plot of Re~‘*. This plot is known 
as the Nyquist plot (4). As an example, Fig. 3 is a Nyquist 
plot of a system, having the following parameters: 


L* = 30 in. C* = 5100 fps 
°F = 1.275 
ay 
7 = 3.52 
T, = kL*/C* = 0.00114 sec 
1 0.000472 
72 sec 
29 Pa — P- 
= —— + ——— ) = 1.19 
2(Wor + \Paor — Pe Pa; — Pe 
90° 
M.R=352 
U=30in 
C’=5100 ft/sec 
a 
T3* 000472 
T,=.00114 
CURVE (d) 
T=.02 
270° 
FIG. 3. NYQUIST PLOT, TYPICAL GAS-PRESSURIZED LIQUID ROCKET 
SYSTEM 
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It is evident from an inspection of curve (a) that stability 
exists for any value of Ky if + is assumed to be zero. For a 
fixed value of gain, curves (b), (c), and (d) demonstrate the 
effect that increasing lag time has on system stability. With 
7 equal to 0.005, the system is stable as indicated by curve 
(b). In the case of curve (d), for 7 equal to 0.02 the system 
is unstable. For 7 equal to 0.01, as shown by curve (c), the 
system is neutrally stable. It is of interest to note that the 
vector Re‘? “spins” forever about the origin due to the ef- 
fect of e~”?. The system response contains an infinite num- 
ber of discrete higher frequencies as the vector “‘spins’’ every 
2x radians and intersects with the real axis. The higher the 
frequency, the smaller the amplitude, indicating increasing 
attenuation with increasing frequency. However, it is pos- 
sible that the system may respond to the higher frequencies, 
provided that the excitation, or forcing function, possesses 
high enough energy and is in phase with these modes. 

(g) Neutral stability is characterized by a unique rela- 
tionship among the system parameters and leads to the con- 
cept of a stability boundary. With all other parameters 
fixed, the relationship which must exist between any two pa- 
rameters in order for Re~** to pass through the —1 point may 
be established. The locus of all points satisfying this relation- 
ship is termed a stability boundary. If variations in more than 
two parameters are considered, families of stability boundaries 
may be found. For example, Figs. 4, 5, and 6 show a family 
of stability boundaries relating Ky and 7 as Ky is varied for 
the system parameters given above. Fig. 4 shows a family of 


stability boundary curves for a,,/a, equal to 0.5 and for vari- 


ous values of #(Pa, — P.)/(P a. — The region below the 
stability boundary curve represents the stable region, and that 
above the boundary the unstable. For typical rocket design, 
the parameter 7(Pa, — P.)/(Pa, — P.) is essentially equal 
to 7, the equilibrium mixture ratio. It appears that for this 
particular example, increasing the equilibrium mixture ratio 
tends to increase the region of stability for a given Ko. Or 
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FIG. 4 STABILITY BOUNDARY PLOT, a@z/ay = 0.5 
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FIG. 5 STABILITY BOUNDARY PLOT, a@r/ay = 1.275 


for a given combustion time lag, higher values of Ky may be 
tolerated as 7 increases. Figs. 5 and 6 illustrate the effect 
of @,,/a, on stability. Comparison of Figs. 4, 5, and 6 indi- 
cates that the higher the value of a,,/a,, the larger the sta- 
bility region. 

(h) The values of the parameters defining a point upon 
stability boundary are such that the characteristic equation 


(Tip + 1)(T~p + 1)(Tsp + 1) 
must be satisfied. To each point on the stability boundary 
there corresponds a frequency at which sustained oscillations 
occur, and this may be related to the system parameters if 
p = iw. And for small values of zw, the characteristic equa- 
tion may be written 


Ko(1 + tw7s)(1 — 


(1 + + iw7,)(1 + iwT's) 


This equation provides an analytical means of locating points 
on the stability boundary. It can be shown that for the bi- 
propellant system, the stability boundary is defined by the 
system parameters and the critical time lag, 7*, as 


7: 
Ko+1 
+ 7,7; + — Kyr*1 


And the critical frequency at the boundary is 


Ky + 1 
T1T2 + T1T3 + T2T3 — Kor*T 


For the monopropellant case, these equations reduce to 


T, = Kor* 
Ky + 1 
and (w*) TT. 
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FIG. 6 STABILITY BOUNDARY PLOT, aoz/ay = 3.0 


where =7T;=T, 


This is the condition attained by Summerfield (1). However, 
it ean be shown that the approximation e~7”? = 1 — rp, for 
small values of rw, is relatively poor in practical cases. From 
a Nyquist plot, the frequency at which the amplitude is 
unity may readily be obtained. And, at this frequency the 
additional phase shift resulting in marginal stability may be 
calculated. Since 7 contributes only phase shift 


w*r* = 180° — oi 
where ¢ is defined by 


(Tip + + 1)(Tsp + 1) 


To illustrate this point, the following example is given. For 
2 bipropellant system having the following parameters 


0.00114 sec 


1 => 
= 0.000472 sec 
7; = 0.00212 sec 
7; = 0.000835 
Ko 1.6 


The approximation yields the critical time lag, 7*, of 0.0019 
sec and the critical frequency, w*, of 215 eps. For an equiva- 
lent monopropellant system, having the same value of Ko, 
the time constant of the feed circuit is equal to T, + 73 or 
().00259 sec, the thrust-chamber time constant, 7, being the 
same. The critical time lag and critical frequency may be 
calculated 


0.00114 + 0.00259 


r* = 16 = 0.0023 sec 
(1.6 + 1) X 10° 
(*) 1.14 X 2.59 
or w* = 150 eps 


It is seen that 7*w* is not small. On the other hand, if the 
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complete Equation [8’] is used, for the bipropellant system, 
the critical frequency and time lag from Fig. 5 may be calcu- 
lated as 80 eps and 0.004 sec, respectively. 


3 Aspects of Stability and Control of a Pump- 
Fed Rocket 


(A) A pump-fed rocket is one in which an oxidizer pump 
and fuel pump are employed to supply propellants to the 
thrust chamber. There are many types of pump drives. 
They differ only in the source of energy supply to the turbine. 
The turbine may derive its power from high-pressure air or 
hot combustion gases obtained from reaction of propellants. 
A typical pump-drive system is one in which the hot combus- 
tion gases from the thrust chamber may be employed to sup- 
ply power to the turbine. This system is of particular inter- 
est from the standpoint of stability and control, because of 
the regenerative feature of such a system. The system is re- 
generative because the thrust chamber supplies combustion 
gases to drive the turbopumps, which in turn supply propel- 
lants that react in the thrust chamber and generate combus- 
tion gases to drive the turbopumps. The analysis presented 
herein will demonstrate that a system which may be unstable 
due to system dynamics and combustion time lag can be 
stabilized by the use of chamber-pressure feedback control. 

(B) The analysis of the simple, gas-pressurized, bipro- 
pellant, liquid rocket system can readily be extended to a 
more complex system such as a regenerative pump-fed sys- 
tem. The rudimentary system is one in which an orifice that 
gives a calibrated hydraulic pressure drop in the hydraulic 
circuit is used to establish a thrust level. Fig. 7 shows a block 
diagram of such a calibrated system. In Fig. 7, G; represents 
the linearized turbine transfer function. If it is assumed 
that the drag, D, on the rotating parts is proportional to the 
square of the angular velocity, the equation for the turbine 
torque may be written as 

T = 16 + D(6)? 
If it is assumed that turbine torque is essentially proportional 
to chamber pressure for small variations of chamber pressures, 
then 
KrP. = 16 + D(6)? 


And the linearized turbine transfer function may be written 


as 


6 =Ks 
where 
_ Kr 
K; = D and 7; = D 


G, and G; represent the linearized transfer functions of fuel 


8 _ Ks 
5 T P+! Re 
WwW. 
' TP 
Pao K; 
3_BP+t 
Ks R 
FIG. 7 REGENERATIVE CALIBRATED TURBOPUMPING SYSTEM 


(BLOCK DIAGRAM) 
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> 
P, 


and oxidizer pumps. Since the pump discharge pressure is 
essentially proportional to the pump speed squared, the 
linearized pump-transfer functions may be written as 

Gy = = Kp, = Ke 
Similarly 
Gy = = Ki 


Therefore, the linearized equations governing the system are 


P, 6 
Wy + Wor Gs P. 
Pa 


The solution of these equations yields as a characteristic equa- 
tion 


1 + G,[G;(1 — + — GsGe)] = 0....... 


For symmetry, it is assumed that G; = G,; Equation [9] 
then reduces to 
1 + G,(G, + G;:)(1 — GsGe) = 0 
or 
1+ Koe*?(Tsp + 1) 
(Tip + 1)(T2xp + 1)(Tsp + 1) 
(Tsp + 1) 
In accordance with the Nyquist criteria, conditions of sta- 
bility may be investigated by examining the behavior of 
Ree~** about the —1 point, where 
Koe~7?(T sp + 1) 
(Tip + 1)(T2p + 1)(Tsp + 1) 
Kako + 1) 
(1 — (; ti 
(Tsp + 1) {11] 


It is to be noted that Equation [11] is similar to Equation [7] 
except for the additional term 


=0 [10) 


(Tsp + 1) 


~ (—7?— 


which results from the inclusion of the turbine and pumps in 
the system. The effect of the addition of the turbine and 
pumps upon previously established conditions for stability of 
a gas-pressurized system may be investigated. There are 
three distinct cases: 


Case KsKe =] 
In this case, Equation [11] may be reduced to 
(Tsp + 1) Tsp 
(Tip + 1)(T2p + 1)(Tsp + 1) (Tsp + 1) 
and the effect of T'sp(7T'sp + 1) is to introduce a low-frequency 
phase lead. High-frequency behavior of the gas-pressurized 
system is not affected by the inclusion of turbine and pumps. 


The criteria for stability are essentially those established for 
the simple gas-pressurized system. 


Case, IT: <1 


The conditions for stability are similar to those established 
for the gas-pressurized system except that conditions must 
now be imposed upon Ky(1 — A;5K¢) rather than upon Ko 
alone. The conditions for stable operation are somewhat less 
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stringent than the gas-pressurized system, since (1 — K;K,) 
is always <1. Not only is the phase lag caused by (Tsp + !) 
more than compensated by the phase lead due to 7;))/ 
(1 — K;K¢s) + 1, but also a higher value of Ko is permissible. 


Case IIT: >1 


In this case, the system is a regenerative feedback on. 
For r = 0, necessary but not sufficient conditions for stability 
are that 

|Ko(1 — KsKe)| <1 
and 


|Ko[T; + — KsKe)]| < 7: + + T3 + Ts 


(A) Ina system defined by a given set of parameters, it s 
desirable to have some means of determining the maximu:. 
value of 7 which the system may tolerate without instabilit 
due to combustion lag time. For this purpose, Equation [‘| ° 
may be rewritten as 


(Tip + 1)(Top + 1)(Tsp + 1)(T'sp + 1) + 
ro(1 K;Ke)(T sp + 1) + = |) 
It is possible to define 


Gp) = eve [ 


where 


g(p) = Ko(1 — KsKe)(Tsp + 1) (ee + 1) 


= (Tip + 1)(T2p + 1)(Tsp + 1)(Tsp + 1) 


The plot of G(p) as p follows a contour about the right half o: 
p-plane is known as a Satche diagram. A general discussio1) 
of the Satche diagram is included in Ref. (3). The vecto: 
e~*” maps into the unit circle; the points of intersection 0! 
[—C\(p)/ge(p)] with this unit circle may be investigated, 
and the critical combustion lag time, 7*, may be calculated 
from the relationship 


where 9* is the phase angle introduced by e~7? at the point 
of intersection. 


FIG. 8 SATCHE DIAGRAM, CALIBRATED CONTROL SYSTEM 


(B) Fig. 3 is a Satche diagram for a system having the 
following parameters 


T, = 1.14 X 10-3 K; = 2.44 X 1072 

T: = 6.26 X 10-3 Ke = 76 

= = K,(K. + K;) = 3.31 
(27's + = 2.35 X 10-3 

7, = Ke + Ki = 2.33 X 10-3 K; = 2.94 X 10-3 

T; = 05 K, = 1.12 X 10-3 


This system is dynamically unstable because AsK, >1, and 
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the inequality of Case III is not satisfied. The condition for 
stability for all values of 7 is that the curve [—gi(p)/g2(p) | 
must never intersect the unit circle. This condition suggests 
the use of feedback not only to improve the system dynamic 
performance but also to make stable operation possible re- 
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FIG. 9 FEEDBACK CONTROL SYSTEM (BLOCK DIAGRAM) 


gardless of the magnitude of 7. If a feedback control system 
of the form shown in Fig. 9 is utilized, the system equations 
become 


Pay 
=" wy + w Pan + AP, = Pa 
Pa. 
G, = —P. G; = Pan + AP; = Pao 
G; = P. P; = Pr 
6 AP 
G; = P. G. => P, 
From which, for G, = G; 
Fe GAG» + G;)G, 


Po 1+ — G.G)(G: + 
and the characteristic equation is 
GAGs G2) G, +1+ Gl G;Gs) (G2 + G;) =0 


1 + G,(G. + G;){G. + (1 — = 0 


which, upon substitution of the proper transfer functions, be- 
comes 


Ko(Tsp + 1)e77? 


+ 1)(Txp + + 1) 
> 
(1 — K;Ke) =~ 
G. + ) {12] 
(Tsp + 1) 


It is to be noted that Equation [12] is similar to Equations 
|6] and [10] except that the open-loop response of the simple 
vas-pressurized system is now modified by terms representing 
not only the turbine and pumps, but also the controller, 
represented by G.. 
If Equation [12] is rewritten 
(Tip + 1)(T2p + 1)(T3p + 1)(T;p + 1) + 

+1) 


[t is now possible to define 


where 
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= Ko Tip + nfe. (Tsp +1) + 


KsKe ks i) 


g'(p) = (Tip + 1)(T2p + 1)(Tsp + 1) (Tsp + 


(C) The value G, may be selected so as to warp the plot of 
in such fashion that intersection with the 
unit circle does not occur. For example, examination of the 
typical system shown in Fig. 8 indicates that G. should be 
selected to alter only the low-frequency behavior of this sys- 
tem. In the frequency range from 0 to 550 rad/sec, G-. should 
be such that the amplitude of [—g’2(p)/g’:(p)] is larger than 
unity. The selection of G. would depend upon a particular 
system. One important consideration would be that the 
selection of G, should be examined with respect to its prac- 
ticability. Thus it is theoretically demonstrated that a 
liquid propellant rocket system may be stabilized com- 
pletely independently of combustion time lag.® 


5 In his recent paper (5), Professor Frank E. Marble of Cali- 
fornia Institute of Technology provides in considerable detail an 
explanation of how one would determine the proper transfer func- 
tion for the feedback controller, showing it is possible to stabilize 
the system for all combustion time lags. 
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The Velocity of Sound in Some Rocket Propellant Liquids 


GEORGE G. KRETSCHMAR'! 
U. S. Naval Ordnance Test Station, Inyokern, China Lake, Calif. 


The velocity of sound in a number of rocket propellant 
liquids was obtained with ultrasonic frequencies. A quartz 
crystal driver was used in an ultrasonic interferometer. 
The interferometer was made with stainless steel parts, 
which permitted measurements to be made on concen- 
trated nitric acid. Data were obtained for various samples 
of nitric acid, hydrazine, and the hydrocarbon mixtures, 
JP-3 and JP-4. 


T IS well known that early work on the velocity of sound 
goes back to the time of Newton who showed that the 
velocity is proportional to [elastic force/density]'/?, and 
from this obtained a formula for the velocity of sound in air 


c = [pressure/density] '/2 


This formula gives a value which is too low to check with 
experimental results. Laplace, about 1803, explained the 
discrepancy by pointing out that the rapidity of the pressure 
transformations does not permit temperature equalization 
to take place, so that Newton’s formula should be modified 
by multiplying the velocity by [y]'/* where y is the ratio 
of the two specific heats (1).2_ Newton’s formula should then 
be written 
c = [y pressure/density] '/: 


The earliest tests of sound velocity in liquids were confined 
to long-distance tests in water (2). Also, considerable work 
was done in the early part of the present century using sta- 
tionary waves in tubes, following the procedures of the well- 
known Kundt tube method as used for sound velocity meas- 
urement in gases (3). 

There are considerable data on sound velocity in liquids 
available in the literature. There are, however, some notable 
gaps, particularly for some liquids used as rocket propellants. 
It was proposed that some of these missing data be obtained. 
The first attempts to obtain such data were made using audio- 
frequencies with simple equipment of the Kundt tube type. 
It was found advisable later, because of difficulties described 
below, to abandon the experiments using audiofrequencies, 
and to build entirely different apparatus, making use of ultra- 
sonic frequencies. 

The general formula for the propagation velocity of a 
plane wave in a fluid is (4) 


c= 


For the case of a liquid medium, the elastic constant K is 
the volume elasticity or bulk modulus. By definition, K = 
V -dp/dV, or it may be written in the reciprocal form, the com- 
pressibility, 8,as1/V-dV/dp. Thatis 


= [K/p]'/2 = [V/p-dp/dV]'/2 = [1/p8]'/2 


If the isothermal compressibility is used, one must introduce 
the factor y, the ratio of the two specific heats, just as in the 
corresponding case for the velocity of a longitudinal wave in 
agas. Thus, the velocity of sound in a liquid is 


v = [/pBis]'/2 or [1/pBea] '/2 
Received September 19, 1952. 


1 Physicist, Research Department, Michelson Laboratory. 
? Numbers in parentheses refer to References on page 84. 
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where 8,, is the isothermal compressibility, and 8,, is the adi:- 
batic compressibility. Most liquids are, however, so near] 
incompressible that y may be considered as unity. Ther: 
are, however, some liquids having considerable compress'- 
bility, and some situations in which it may be advisable t» 
take y into account. 


Experiments 


The Kundt tube method, using audiofrequencies, appeare:' 
to be a natural and attractive method at the beginning o: 
this work because the necessary equipment was readily availa- 
ble. There were tubes of various materials and lengths, « 
good audio oscillator, an amplifier, and a heavy-duty loud- 
speaker. The method chosen was to make use of a constant! 
length tube with a variable applied frequency. The output 
of the oscillator was amplified and used to drive the speaker 
which was directly connected to a small piston fitted into onc 
end of the tube. At the other end of the tube was fitted «a 
piston containing a thin diaphragm, behind which was 
mounted a pair of barium titanate disks immersed in oil. The 
barium titanate disks, connected in push-pull, acting through 
a preamplifier and oscilloscope, made a sensitive pickup 
system. 

First experiments with the Kundt tube equipment revealed 
a wealth of resonance points, most of which did not corre- 
spond to simple longitudinal modes of the liquid contained 
inthe tube. It was very apparent that many spurious modes 
were excited. Probably many of these were due to transverse 
and longitudinal modes of the tube itself. An attempt to 
reduce these by using a thick-walled tube of stainless steel 
was unsuccessful. Also, surrounding the glass tube with a 


FIG. 1 
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FIG. 2 ULTRASONIC INTERFEROMETER. 
DISASSEMBLED VIEW 
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thick layer of wet sand did not sufficiently attenuate the 
unwanted modes. 
Distilled water was used in order to have a liquid in which 
the sound velocity was accurately known. Working back 
from known values, it was possible to recognize the resonance 
points corresponding to the liquid column modes. But, 
without this help it was difficult if not impossible, to make a 
logical interpretation of the data.* Plotting of a complete 


FIG, 3. INTERFEROMETER CROSS SECTION 


set of data gave a completely ambiguous result, for one could 
select a number of sets of points, each of which would 
give a straight line relationship between the frequency and the 
number of half wave lengths, but each corresponded to a 
different value of the velocity of sound in the liquid. 

It was at this point that the decision was made to drop the 
Kundt tube method and to build an ultrasonic interfer- 
ometer, making use of a quartz crystal as the driving element 
(5, 6). Such an equipment makes use of standing waves 
just as in the Kundt tube. But, by making use of a rigid 
and accurately fitted moving plate and a good micrometer, 
« much higher degree of precision may be attained than with 
the Kundt tube, in spite of the much smaller size of the parts. 
There is, also, the additional advantage that a very small 
amount of liquid is required, and also the whole interfer- 
ometer may be immersed in a constant temperature bath for 
accurate temperature control. 

The instrument, as constructed, is shown in Fig. 1. On the 
left is the tube through which air pressure can be admitted 
below the pistonlike crystal holder. This enables one to 
compress the crystal slightly against the diaphragm. Fig. 2 
is a disassembled view in which the quartz crystal appears 
as the solid disk on the left front. The transparent parts 
are insulators, machined from lucite. The piston can be 
seen projecting from the base cylinder. To the right are the 
stainless-steel cylinder and the cover which contains the mova- 
ble plate and the micrometer head. 


3 Data have, however, been obtained by the method of reso- 
nances for a liquid contained in a large glass sphere. Ina private 
communication, Mr. David Lieberman, of Aerojet Engineering 
Corporation, has informed the author that he obtained a value 
for the velocity of sound in nitric acid of 1.25 & 10° em/sec by 
this method. 
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Fig. 3 is a detailed cross section of the interferometer. 
The stainless steel cylinder holds a 5-mil thick stainless-steel 
diaphragm against a shoulder in the base cylinder along the 
section AA, The crystal is just under and in contact with the 
diaphragm. The chamber is made tight with teflon gaskets 
at AA and BB; also a Vinylite O-ring around the stem of the 
movable plate prevents the escape of corrosive gases or 
fumes and effectively seals off the liquid compartment. The 
seals are made more secure and crystal contact is assured by 
spreading a thin layer of silicone grease on the gaskets and 
over the face of the crystal. 

The crystal used was an X-cut quartz disk about 2 cm in 
diameter and about 6 mm thick. The resonant frequency 
was 512 ke. It was driven by a single 6L6 tuned amplifier 
stage which was in turn excited by the output of a General 
Radio Type 700-A wide range beat frequency oscillator. 

Fig. 4 is a circuit diagram of the amplifier stage. The con- 
stants have been adjusted so that the tube operates with a 
plate current of approximately 50 milliamperes. In use there 
are slight variations of current in the plate circuit and this 
effect is made use of to determine the resonant positions. 
In order to follow these variations accurately, thus getting 
the positions of resonance with precision, a microammeter 
is introduced into the plate circuit, most of the current 
being neutralized by a bucking voltage across the microam- 
meter. This part of the circuit is shown in more detail in 
Fig. 5. The microammeter is paralleled with an equivalent 
resistance of about 600 ohms so that the circuit may be first 
adjusted so that the microammeter will be on scale. The en- 
tire circuit was well shielded and the alum’num box grounded, 
but it was found necessary to leave the dry cell, used for pro- 
ducing the neutralizing voltage, outside of the box. 

The frequency was determined and controlled with the aid 
of a Signal Corps Type BC 221 AH frequency meter. This 
was loosely coupled to the amplifier with a few turns wound 
around the crystal cable. Fig. 6 shows the interferometer 
in use with the required electronic gear. The voltage of the 
rectifier and the oscillator was stabilized with a Sorensen 
voltage regulator. The frequency meter can be seen at the 
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FIG. 4 CIRCUIT DIAGRAM OF 6L6 TUNED AMPLIFIER STAGE 
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FIG. 5 NEUTRALIZING CIRCUIT AND MICROMETER INDICATOR 
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FIG. 6 INTERFEROMETER CONTAINED IN CONSTANT TEMPERATURE 
BATH; CONNECTED TO ACCESSORY ELECTRONIC EQUIPMENT 


top of the photograph, and the coupling wire can also be 
seen. It will be observed, also, that the interferometer is 
immersed in a constant temperature bath. The pressure 
applied to the crystal is obtained from the air line and is 
controlled by a regulator and leak. A pressure of a few inches 
of mercury was needed for optimum crystal contact. The 
frequency meter was itself crystal-controlled with a standard 
100-ke crystal, and the crystal was checked by beating with 
a 100-ke standard frequency signal available in the labora- 
tory. 


Results 


Data were obtained on a number of rocket propellant 
liquids. A condensation of these data is given in Table 1. 
Specific gravity determinations were done on all the liquids 
studied, and these are included in the table. More data 
were taken on hydrazine than on other liquids studied be- 


TABLE 1 AVERAGE VALUES OF THE VELOCITY OF SOUND IN SOME 
ROCKET PROPELLANT LIQUIDS AT 512 Kc 
(Commercially available materials) 


Measured Measured 
density sound 
Liquid Composition at 25 C velocity at 
25 C, m/sec 
Hydrazine (96+ '/2)%anhydrous 1.0096 2059 + 6 
RFNA? HNO; 77% min. 1.585 1379.5 +1.4 
NO, 20% approx. 
H:O 2'/. max. 
RFN A? HNO; 90% min. 1.479 1340.4 + 3 
H,0 10% max. 
WFNA HNO; 97.08% 1.494 1239.2+0.5 
NO, 0.22% 
1.91% 
WFNA HNO; 97.52% 1.497 1239.7 +0.5 
NO, 0.28% 
2.11% 
JP-3 
Samplea MIL-F-5624A° 0.7682 1233.8+2 
JP-3 
Sample b MIL-F-5624A 0.7590 1196.54 1 
JP-4 
Samplea MIL-F-5624A 0.7757 «1242.24 1 
JP-4 
Sample b MIL-F5624A 0.7762 1242.3+1 


“ Reagent grade, A.C.S. Specifications, Code 1122. 

> Reagent grade, A.C.S. Specifications, Code 1121. 

© Specifications may be obtained from Air Materiel Command, 
Wright-Patterson Air Force Base, Dayton, Ohio. 
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cause of the lack of sharpness of the resonances, with conse- 
quent greater uncertainty in wave-length measurements. 
The errors are estimated from the consistency of the data and 
the number of runs made on each liquid. All of the data 
were taken close to the resonant frequency of 512ke per secon. 

As a check on the over-all accuracy of the interferomete,, 
measurements were made on distilled water at 25 C. The 
average of three measurements made is 1498.22 m/sev. 
This may be compared with the recently published value «f 
Barthel and Nalle, which is 1496.2 + 0.2 m/sec, and that «f 
Willard, which is 1498 + 1 m/sec at 25 C (7, 8). 
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ARS 1953 Junior Award Competition 


Papers are now being invited for consideration for 
the 1953 Junior Award, to be presented at the ARS 8th 
Annual Convention to be held in December 1953. 

The winning paper will be judged mainly on the basis 
of content which should reflect original thought and ef- 
fort. The age of authors of papers should not exceed 25 
years. 

Please refer to page 57, this issue of the JouRNAL, for 
acceptable scope of subject matter. Manuscripts are 
to be prepared in accordance with style instructions 
on inside back cover. 

Papers must be received no later than October 1, 
1953, and should be clearly marked, ‘‘Submitted for 
Junior Award Competition.” 

Send papers to: The Secretary, American Rocket | 
Society, 29 West 39th Street, New York 18, N. Y. 
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On Maximum Evaporation Rates of Liguid Droplets in 
Rocket Motors 


S. PENNER? 


Daniel and Florence Guggenheim Jet Propulsion Center, California Institute of Technology, Pasadena, Calif. 


Upper limits have been estimated for the rate of evapo- 
ration of small liquid droplets in representative rocket 
combustion chambers. The droplets are assumed to be 
isothermal at all times. The droplet temperature as a 
function of time is determined by an appropriate heat 
balance. The calculations are useful in determining the 
significance of inelastic collisions between liquid droplets 
fer complete combustion. Radiant heat transfer to mov- 
ing liquid droplets is considered briefly. 


I Introduction 


4 XPERIMENTAL studies of liquid fuel rockets have 

Ms shown that the performance of many bipropellant sys- 
tems is dependent on the type of injection system employed. 
The results indicate the importance of (a) the size of the 
liquid droplets at the point of injection, and (6) the relative 
directions of injection of the separate streams of fuel and 
oxidant. 

These experimental findings suggest a combustion mecha- 
nism in liquid fuel rockets which may be represented sche- 
matically as follows: 


Injection 
of 
Fuel and Oxidant 
A 
B 
Collisions 


Between Streams Vaporization 


of Liquid Droplets 


Liquid-Phase Reactions 
Gas-Phase Reactions 


Vaporization 


Gas-Phase Reactions 


Combustion Products 


Mechanisms A and B generally occur concurrently in a motor 
chamber. For any given set of conditions, either mecha- 
nism A or mechanism B may be of greater importance. 

An indication regarding the relative importance of mecha- 
nisms A and B may be obtained by determining whether liquid 
droplets are in a given motor chamber for a sufficiently long 
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period of time to vaporize completely without the occurrence 
of exothermic liquid-phase reactions. If the time required 
for the complete vaporization of the liquid droplets is greater 
than the time spent by the liquid droplets in the motor cham- 
ber, then liquid-phase reactions must occur, provided com- 
plete combustion takes place in the rocket chamber. It is 
therefore of interest to determine the minimum time required 
for the evaporation of liquid droplets in a rocket chamber.* 

The evaporation of liquids and solids has been studied both 
experimentally and theoretically (1-11). Extended theo- 
retical investigations of vaporization under idealized condi- 
tions have been carried out by Fuchs (5). A general treat- 
ment for the calculation of droplet lifetimes in combustion 
chambers has been described by Wentzel (6). Application of 
the results of these studies on evaporation to the special con- 
ditions existing in a rocket chamber cannot readily be made. 
The treatment of Wentzel (6), although applicable in prin- 
ciple, requires knowledge of a number of empirically de- 
termined parameters which are not generally available. 
Furthermore, there are indications that the usual laws which 
govern the evaporation of liquid droplets may not be valid for 
droplets of exceedingly small size. For example, Probert (1) 
indicates that the rate of evaporation of very small droplets is 
proportional to the diameter of the evaporating droplet rather 
than to the square of the diameter. Studies on compara- 
tively large drops (0.5 mm in radius) (4) led to the conclusion 
that the treatment of Fuchs (5) for the rate of evaporation of 
drops is essentially correct. The work of Fuchs, in so far as 
it applies to very low Reynolds numbers, has been criticized 
recently (11). Finally, a study of evaporation rates for fall- 
ing and burning droplets by Topps (12) shows a complicated 
dependence for evaporation rate on size and Reynolds num- 
ber. 

In view of the obvious complexity of the problem, no at- 
tempt was made to carry through a complete treatment of 
the rate of evaporation of liquid droplets in a rocket chamber. 
Instead, an approximate treatment has been employed which 
is designed to estimate a reasonable upper limit for the rate 
of droplet vaporization in a rocket chamber. A discussion 
of the maximum and minimum rates of evaporation, under 
conditions somewhat different from those of interest for the 
present study, has been given previously (13). 

The results of the present investigation emphasize the im- 
portance of droplet size. For example, at least 1.58 X 107% 
sec is required before the radius of an aniline droplet decreases 
from 5 X 107-8 em to 2.5 X 107% em. However, more than 
1 X 10? see is required before an aniline droplet initially 


3 It should be noted particularly that the present investigation 
was undertaken with the modest objective of approximating 
minimum lifetimes of droplets rather than with the idea of ob- 
taining a realistic estimate of the complicated evaporation 
processes mabe’ in heterogeneous combustion. For recently 

published reports on heterogeneous combustion see, for example, 

E. C. Topps, J. Inst. Petrol., vol. 37, 535 (1951 ) Go Ak. 
Godsave, National Gas Turbine Establishment Reports Nos. R. 
66, March 1950, and R. 87, April 1951, Pyestock, Hants, England; 
D. B. Spalding, Fuel, vol. 29, 2, 25, 1950; vol. 30, 121, 1951. 

4 Numbers in parentheses refer to References on page 88. 
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1 X 10-? cm in diameter decreases to one half of its original 
diameter. The rate of evaporation is, of course, dependent 
on the nature of the evaporating compound. Thus octane is 
found to evaporate more rapidly than aniline under similar 
experimental conditions. The transfer of radiant energy to a 
moving liquid droplet is discussed. The analysis shows that 
the maximum possible temperature rise associated with the 
absorption of radiant energy increases as the droplet diameter 
is decreased. 


II Outline of Theory 


The rate of droplet vaporization depends upon the nature 
of the evaporating liquid and upon the heat transferred to the 
liquid droplets after injection into the combustion chamber. 
The minimum life of a liquid droplet in a rocket chamber may 
be calculated by using an upper limit for the heat transfer and 
also for the isothermal rate of decrease of droplet radius with 
time. 

The order of magnitude of the upper limit for the isothermal 
rate of decrease of droplet radius with time can be calculated 
from the Knudsen equation, which may be written in the 
form (14-16) 


—dr/dt = (p./pi(M/2xRT)°/® {1] 


where r is radius of the evaporating droplet, ¢ time, ps satu- 
rated vapor pressure of the evaporating liquid at a given tem- 
perature, o; density of the evaporating liquid, M molecular 
weight of the evaporating compound, 7’ temperature, R gas 
constant. An upper limit for the heat transfer to the liquid 
droplets is obtained by maximizing the product hAT in the 
relation 


where Q is rate of heat transfer to a spherical droplet of radius 
r, h is over-all coefficient of heat transfer to a liquid droplet 
evaluated at the mean film temperature, A7’ is mean tempera- 
ture difference between the gases surrounding the liquid 
droplet and the liquid droplet itself. 

The order of magnitude of the value of h may be ascer- 
tained from an empirical relation determined as the result of 
heat transfer studies between air and spheres (17, p. 237),5 viz., 


(2hr/ky) = [3] 


where ky is thermal conductivity of the hot gases at tempera- 
ture 7,, G = mass rate of flow of air past the liquid droplets, 
us = viscosity of hot gases at the temperature T;, T; = T + 
0.5 (T. — T) where T is the droplet temperature (assumed to 
be uniform), 7, = temperature of the hot gases. At a chamber 
temperature of 3000° K and a pressure of 300 psi, it is found 
that Equation [3] leads to values of 2hr/k, of about 4 for the 
largest droplets (r = 10~? em in radius), even if the droplets 
move at 100 fps relative to the gases in the chamber. For 
very small droplets, on the other hand, 0.33 (2rG/yy)%® <2. 
However, the theoretical lower limit for 2hr/ky is 2 in a stag- 
nant medium unless the spheres have diameters of the same 
order of magnitude as the mean free path in the fluid. There- 
fore, it is evident that 2hr/k, for small droplets does not differ 
greatly from 2, and considering the fact that the gases at the 
head end of the rocket chamber move, on the average, in the 
direction of the droplets,* it appears likely that even for r = 
10-* em, the approximation h = k;/r will yield a reasonable 
estimate for heat transfer to the droplet.? Thus we set 


5 It should be noted that the conclusion reached by using 
Equation [3] is the same as that obtained from the elegant treat- 
ment on heat transfer to spheres for very low Reynolds numbers 
given by Kronig and Bruijsten (11). 

6 Representative estimates of gas velocities near the injector 
are of the order of 50 fps. 

7 This last conclusion does not make allowance for the possible 
occurrence of turbulent motion near the injector. 
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FIG. | TEMPERATURE-VS.-TIME AND RADIUS-VS.-TIME DATA FO% 
AN ANILINE DROPLET IN A NITRIC-ACID-ANILINE ROCKET MOTO, 
(T, = 3000° K, To. = 300° K, ro = 5 X 1075 cm). 


AT = 3000 — T withh = ky/r........... [3a 


The rate of temperature rise with time, (dT /dt),, of 
spherical droplet receiving the energy Q is 


(dT /dt), = Q(4rr3/3) = 8hAT/rpCi..... {4 


where C; is the heat capacity per unit weight of the evaporat- 
ing liquid. The liquid droplet is cooled as the result of evapo- 
ration. The rate of decrease of temperature with time asso- 
ciated with the cooling effect of evaporation is given accord- 
ing to the relation 


—(dT /dt), = (8AM vap/CiMr)( —dr/dt)......... [5 


where AHva, is the heat of vaporization of the liquid pe: 
mole, and (—dr/dt) is given by Equation [1]. The net rate 
of change of temperature with time is therefore 


dT /dt = — (84H vap/rMC.i)( —dr/dt)... [6 


where use has been made of Equation [8a]. 

Reference to Equation [6] indicates that the temperature 
of an evaporating liquid droplet changes continuously with 
time. Since the first term on the right-hand side of Equation 
[6] varies nearly inversely with the square of the radius, 
whereas the second term varies inversely as the first power of 
the radius, it is apparent that the droplet temperature will 
rise rapidly for very small droplets. For moderately large 
droplet sizes, Equation [6] is found to lead to quasi-steady- 
state temperatures. 


Ill Calculation of Minimum Lifetimes of 
Droplets 


Equations [1] and [6] can be solved numerically without 
difficulty. A representative calculation is outlined below 
for the evaporation of aniline in a nitric-acid-aniline rocket 
operating at a chamber temperature of 3000° K. 

The aniline droplets are chosen to have an initial diameter 
of 5 X 10-*em. They are injected at an initial temperature 
of 300° K with an average velocity of 100 fps into a rocket 
chamber 1 ft in length. The droplet velocity is assumed to 
remain constant. Thus the droplets remain in the rocket 
chamber for 0.01 sec. It is therefore of interest to determine 
the temperature and radius of the aniline droplets as a func- 
tion of time for 0.01 sec. 

The gaseous combustion products in the nitric-acid aniline 
rocket consist predominantly of NO, CO, COs, H.O, and N». 
Approximate calculations for these components lead to a 
value for yy of about 0.157 lb/hr-ft (p. 411 of Ref. 17) at the 
temperature T, = 300 + 0.5 (3000 — 300) = 1650° K. For 
the specified composition, it is found (p. 415 of Ref. 17) that 
0.76. Since Cp 0.50 Btu/Ib-°F (p. 406 of Ref. 
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Fl.. 2. TEMPERATURE-VS.-TIME AND RADIUS-VS.-TIME DATA FOR 
A ANILINE DROPLET IN A NITRIC ACID-ANILINE ROCKET MOTOR 
(T. = 3000° K, 7) = 300° K, ro = 1 X 107? cM). 


17), it follows that ky ~ 0.104 Btu/hr-ft-°F = 4.30 107! 
cx|/sec-em-°K. 

For approximate calculations the numerical values of ky, 
p;. C;, and AHy., are assumed to be independent of tem- 
perature. Since the droplet temperature does not vary 
greatly with time, the error introduced by this simplification 
is not large. Thus, 7 = 3000° K, p ~ 1.02 gm/em’, C; 
0.53 eal/gm-°C, M = 93.12 gm/mole, AHva, ~ 12,000 
(18). 

Introduction of the numerical values into Equation [6] 
leads to the desiied result 


di'/dt = 2.88 X 10-3[(3000 — 7)/r?] — (730/r)(—dr/dt). . [6a] 


The values of —dr/dt may be calculated directly from Equation 
||] by use of measured values for the vapor pressure of aniline. 
For numerical work it is convenient to express —dr/dt, as a 
first approximation, by the relation log (—dr/dt) > A — 
BT, where A and B are nearly independent of tempera- 
ture. For aniline it is found that 


log (—dr/dt) ~ 7.09 — 2.60 X 10°7?......... [7] 


Temperature-time and radius-time curves can be constructed 
from Equations [6a] and [7] by using a straightforward 
iteration procedure. Representative results for aniline, 
hydrazine, and octane are shown in Figs. 1 through 5. The 
same numerical value for ky was employed for all of the caleu- 
lations. 


IV Results 


The results plotted in Figs. 1 to 5, inclusive, exhibit the 
expected behavior. The temperature of the liquid droplets 
increases rapidly at first, then reaches a quasi-steady state, 
and finally rises very rapidly after the droplet radius has be- 
come sufficiently small. The droplet radius decreases con- 
tinuously with time, very slowly at first, and then at an ever- 
increasing rate as the droplet temperature rises. From Figs. 
| to 5 it is possible to ascertain the minimum lifetimes of 
droplets for a number of different droplet sizes and liquids. 
The calculated results may be compared by evaluating the 
half life, ta/2, of the liquid droplets, i.e., the time required 
to decrease the droplet diameters to one half of their original 
values Relevant data are summarized in Table 1. 

Reference to the results summarized in Table 1 shows that 
an aniline droplet with an initial radius of 1 X 10~* em can- 
not evaporate in 0.01 sec unless heat transfer to the liquid 
droplet is considerably greater than that calculated from 

8 The numerical value of t(1/.) is always greater than one half 
he time required for complete evaporation of a given droplet. 
Calculation of the time required for complete evaporation is 
somewhat laborious. 
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FIG. 3. TEMPERATURE-VS.-TIME AND RADIUS-VS.-TIME DATA FOR 
AN OCTANE DROPLET SURROUNDED BY HOT GASES AT 3000° K 
(To = 300° K, ro = 5 X 10-%cm). 
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FIG. 4 TEMPERATURE-VS.-TIME AND RADIUS-VS.-TIME DATA FOR 
AN OCTANE DROPLET SURROUNDED BY HOT GASES AT 3000° K 
(To = 300° kK, ro = 1X 10-2 cm). 
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FIG. 5 TEMPERATURE-VS.-TIME AND RADIUS-VS.-TIME DATA FOR 
A HYDRAZINE DROPLET SURROUNDED BY HOT GASES AT 3000° K 
(T> = 300° K, ro = 1 X 1072 cm). 


Equations [2] and [3]. Since only lower limits for the drop- 
let lifetimes have been calculated, it is not possible to draw 
definite conclusions for the smaller droplet sizes for aniline. 
The evaporation of octane is seen to proceed more rapidly 
and at lower droplet temperatures than the evaporation of 
aniline. Hydrazine appears to evaporate somewhat more 
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TABLE 1 MINIMUM HALF-LIFE OF LIQUID DROPLETS 
(T, = 3000° K, 7) = 300° K, h = ky/r, ky = 4.3 X 10-4 cal/sec- 


em-°K) 
Compound rp X 10°(em) = (1/2) X 102(sec) (1/2) (°K) 
Aniline 5 0.158 392 
10 0.560 370 
Octane 5 0.052 352 
10 0.188 333 
Hydrazine 10 1.07 317 


@ t(,/2) = time required for the droplet diameter to decrease to 
one half of its original value. 
> 7 (1/2) = droplet temperature at ¢(1/2). 


slowly than aniline. Because of the uncertainties which 
exist with regard to droplet size distribution in rocket motors,° 
it is not possible at this time to draw definite conclusions with 
regard to the importance of liquid phase reactions in given 
rocket motors. However, it is evidently desirable to extend 
the calculations presented in this discussion to larger droplets 
by using the known size and velocity distributions in com- 
bustion chambers in conjunction with realistic heat transfer 
equations. 


V_ Radiation Heating of Liquid Droplets 


Equations for the heat transfer by radiation to a liquid 
droplet moving along the axis of a rocket motor chamber 
have been derived previously (19). It was shown that the 
radiant energy, I7, received per square centimeter of surface 
during the residence time ¢, in the rocket chamber, is given 
by the approximate relation 


= + wif, dn dt fu exp {—k, pla? + 
(L — vt)?]Q/2)} Jada + {l—- exp[—k,e(a? + 
2*)(1/2) act [8] 


where a is radius of cylindrical motor chamber (em); L 
length of chamber (cm); J°, radiant energy (cal/cm?-sec) 
emitted by a blackbody radiator in the wave-length region 
between \ and A + dad, the blackbody being at the same 
temperature as the radiating gases surrounding the moving 
liquid droplet; ¢ is time (sec); t, residence time of the liquid 
droplet in the motor chamber (sec); p density of the hot 
radiating gases (gm/cm*); k, spectral mass absorptivity of 
the gases (em?/gm); v linear flow velocity of the moving 
liquid droplets (cm/sec). 

Approximate calculations of heat transfer by radiation to 
moving liquid droplets have been carried out for the case 
where a = 7 cm, L = 30cm, T, = 3000° K, f°, dd = 115 
cal/em*-sec, t; = 0.01 sec, kx = 40 cm?/gm, p = 2 X 1073 
gm/cem’, v = 3 X 10% cm/sec. Evaluation of J7 from Equa- 
tion [8] for the specified numerical values leads to the result 
I, = 0.733 cal/em?. The maximum possible radiant energy, 
I’, incident on a liquid droplet surrounded by blackbody 
radiators for 0.01 sec, is I’r = 1.15 cal/em? if the blackbody 
radiators are at 3000° K. This is only about 1.6 times the 
value calculated from Equation [8]. However, the calcu- 
lated value of Jr may not be significant, since the assumed 
value of kk = 40 em?/gm may be too high for a liquid fuel 
rocket operating at a pressure of 300 psi (19). 

The temperature rise AJ of a spherical liquid droplet of 
radius r, caused by incident radiation J7) cal/em? in the wave- 
length region between \ and \ + dd is given by the rela- 
tion (19) 


= (3/rCip:) — exp(—2K,r)]da...... [9] 


® Dr. J. M. Schmidt, of JPL, CalTech, estimates the radius of 
aniline droplets at the point of injection into a nitric-acid-aniline 
rocket as probably lying between 5 X 10% and 5 X 107? cm. 
Our analysis does not apply to the largest droplets. 
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where K) is the absorption coefficient (in em~') of the liquid 
droplet for light in the wave-length region between A and \ — 
dy. Evaluation of AT, by use of Equation [9] can be per- 
formed only if K, is known as a function of wave lengtl. 
Experimental determination of K, could be carried out in the 
same manner as was done previously for double-base rocket 
propellants (20). 

If an average value K may be used for Ky, then AT’ is 
given by the relation 


ATr = 3I7[1 — [10 


For very small droplets and weak absorbers [1 — exp(—2Kr 
=~ 2Kr. Therefore 


The radiation temperature rise of small droplets, which ar: 
not blackbody absorbers, is therefore approximately propor- 


tional to K. An upper limit may be estimated for AT b: 
treating the absorbing liquid droplet as a blackbody. Ih 
this case 

= = 2.2r—!, °K......... [lla 


if Cp, 1 cal/em*— °K, and = 0.733 cal/em?. 

It can be seen from Equation [lla] that the maximum pos- 
sible temperature rise of liquid droplets caused by the ab- 
sorption of radiant energy increases as the droplet size de 
creases. If Equation [lla] holds, then ATp,,,. ~ 300° h 
forr = 5 X 10~*em, and — 300° K for r = 5 X 
em. In order to evaluate the relative importance of radiant 
heat transfer, we may compare the coefficient AT7'g/t, with 
(dT /dt), as given by Equation [4]. For very small droplets it 
appears that the equilibrium transfer of radiant energy in 
rocket motor is probably not very significant. However, 
the quantitative effects of chemiluminescent radiation may 
well be important and merit extensive further study. 
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Combustion Studies in Rocket Motors 


KURT BERMAN' and SAMUEL H. CHENEY, JR.’ 


General Electric Company, Malta Test Station, Ballston Spa, N. Y. 


investigation of combustion instability was continued 
with 3-in. diam, 1200-lb thrust rocket motors. The study 
was carried out both with normal developmental motors 
and the slit-window motor. Pictures of shock propagation 
axially along the motor are presented. The reactants 
used were ethyl alcohol and liquid oxygen. 


Introduction 


N A PREVIOUS paper (1),* * experiments with a rocket 

motor having a full-length observation window were de- 
scribed. Special emphasis was devoted to the presentation of 
plotographs showing combustion under unstable operation. 
It is the purpose of this paper to describe some recent experi- 
ments both with the window-motor assembly and normal 
developmental water-cooled motors which present further 
data on the phenomenon of instability. The reactants used 
in all experiments were ethyl alcohol and liquid oxygen. 


Water-Cooled Body Experiments 


As was pointed out in (1), two.general types of instability 
have been observed with 3-in. diam motors: (a) A high- 
frequency type of instability of about 1400-1700 eps, and 
(b) a low-frequency type of about 200-300 cps. As men- 
tioned in (1), 3-in. modified impinging head assemblies tend 
toward high-frequency instability as the chamber pressure is 
lowered. In order to study the effect of chamber length on 
the frequency of oscillation incurred in this type of in- 
stability, a modified impinging injector head was assembled 
with a 41-deg convergent half angle, 1.75-in. throat diam, 
and 15-deg divergent half-angle nozzle and various lengths of 
3-in. diam cylindrical water-cooled chambers. The normal 
instrumentation was supplemented by a GM-type Dynagage 
pressure pickup‘ mounted in a boss inserted in the °/j¢-in. 
diam igniter hole of the head. The diaphragm of the gage 
was located approximately 2'/2 in. from the head face. With 
the various assemblies, the chamber pressure was lowered suf- 
ficiently (to approximately 270 psig) to induce the unstable 
operating condition. 

Two general modifications of high-frequency chamber pres- 
sure oscillations were encountered: (a) The type shown in 
Fig. la which has been arbitrarily designated as “‘sinusoidal- 
type,’ and (b) the oscillations illustrated in Fig. 1b which have 
been designated as “shock-type.’”’ The latter type was in- 
variably accompanied by large chamber heat-transfer rates 
(1.5 to 2.5 times normal). Since a run is usually terminated 
us soon as a predetermined heat-transfer rate is reached, this 
caused the shutdown of most tests in which instability was 
characterized by shock-type oscillations. In one test the 
motor was permitted to run under unstable conditions for 29 
sec. Although there was considerable mechanical shock on 
the motor assembly, as evidenced by the fact that several 
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FIG. la DYNAGAGE CHAMBER PRESSURE RECORD OF “SINUSOIDAL’ 
OSCILLATIONS 


FIG. 1b DYNAGAGE CHAMBER PRESSURE RECORD OF “‘SHOCK- 
TYPE’’ OSCILLATIONS 


gage connections were broken, and the body heat-transfer 
rate was quite high (>2 Btu/in.? sec), the motor was un- 
harmed. The nozzle heat transfer rate during this run was 
about 25 per cent higher than normal. It appears therefore 
that the high heat-transfer rates obtained with shock-type 
oscillations are not caused by local hot spots. 

Since the Dynagages are calibrated statically for pressure 
and are temperature-sensitive, all pressure amplitude meas- 
urements must be accepted with some degree of caution. 
The peak-to-peak amplitude of the shock-type oscillations 
ranged from 200 to 500 psi. 

Sinusoidal-type oscillations were generally much smaller in 
amplitude, usually not much larger than 100 psi peak-to- 
peak. No abnormal heat-transfer rate seemed to be asso- 
ciated with sinusoidal-type oscillations. Shock-type in- 
stability was always preceded by periods of sinusoidal oscilla- 
tions. A Dynagage record of the transition period is shown 
in Fig. 2. 

The frequency of the oscillations varied with chamber 
length. The predominant components of the frequency spec- 
trum were the fundamental and a second harmonic. Sinu- 
soidal oscillations were considerably lower in frequency than 
the shock type. For example, with the 21%/;-in. assembly, 
the frequency of the sinusoidal oscillations was 830 to 860 
cps,'while for the shock-type, the frequencies ranged from 930 


FIG. 2. DYNAGAGE CHAMBER PRESSURE RECORD OF TRANSITION 
FROM SINUSOIDAL TO SHOCK-TYPE OSCILLATIONS 
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FIG. 3 FUNDAMENTAL FREQUENCY OF SHOCK-TYPE CHAMBER 
PRESSURE OSCILLATIONS VERSUS CYLINDRICAL CHAMBER LENGTH 
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FIG. 4 SLIT PHOTOGRAPH OF JET AT UNSTABLE OPERATION 


to about 1000 eps. in Fig. 3 the full-stage “shock-type”’ 
fundamental frequency is plotted as a function of the cylin- 
drical chamber length. 

A General Radio Type 65] continuous strip camera was 
used in some experiments to photograph the rocket motor jet. 
The camera viewed the jet through a slit and the film motion 
was perpendicular to the jet axis. Photographs of the dia- 
monds during stable and unstable operation are illustrated in 
Figs. 4 and 5. The frequency of the diamond oscillation cor- 
responded to that of the fundamental chamber oscillation 
frequency recorded by the Dynagage. 

As the operating chamber pressure for a fixed motor as- 
sembly is lowered, the propellant flow rate and hence the head 
pressure drop decreases. To determine whether the onset of 
instability is primarily a function of chamber pressure or head 
pressure drop, experiments were performed keeping the mass 
flow through the motor constant and changing the chamber 
pressure by using nozzles of various throat diameters. Of 
course the gas velocity in the chamber changes also as the 
throat to chamber area is varied. This factor is believed to be 
of minor importance as a stability parameter in the present 
tests. 

Motor assemblies consisted of a modified impinging head, 
a 15'/2-in. long, 4!/2-in. diam water-cooled body, and nozzles 
of 2.5, 3.0, 3.5, and 3.75-in. throat diam. The 4!/2-in. diam 
assembly was used, since insufficient 3-in. diam assemblies 
were available to extend the experiments over a wide range of 
chamber pressure. In Fig. 6 the chamber pressure is plotted 
versus the total flow rate of propellants. Runs in which inter- 
mittent instability was experienced are indicated. Although 
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the severity of the instability decreased with increasing throat 
diameter, it is interesting to note that the instability region 
begins when the total flow rate is below a definite value (about 
9 |b/see in this specific instance). This seems to indicate 
that head pressure drop is a more important factor in deter- 
mining stable operation with the modified impinging hea: 
than is the chamber pressure. 


Window Motor Tests 


Modified Impinging Head, 41-Deg Convergent Nozzle, 17- 
In. Chamber Tests 


Photographs of unstable combustion obtained with the 10- 
in. long window-motor chamber were presented in (1). In 
order to continue the investigation at various chamber lengths, 
the window motor was attached to standard 3-in. diam wate - 
cooled bodies. To date, tests have been conducted with tw» 
lengths of composite chamber: one about 17 in., the other ap- 
proximately 21 in. long. Two different injector heads wer: 
employed—a modified impinging and a conventional imping- 
ing head of the types described in (1). The effect of the con 
vergence angle of the nozzle was introduced as a variable by 
using 1.75-in. throat diam nozzles having either a 10 or a 41- 
deg convergent section. <A few tests were also conducted wit): 
a “stove-pipe motor”’ (i.e., one having no nozzle). The equip- 
ment and procedures were otherwise identical to those de- 
scribed in (1), except that the instrumentation was supple- 
mented by the use of Dynagages to measure chamber pressure 
and, at times, alcohol line-pressure oscillation frequencies. 
Timing marks of either 120 eps or 240 cps were provided on 
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FIG. 5 SLIT PHOTOGRAPH OF JET AT STABLE OPERATION 


the film edge by means of an externally excited spark gap. 

Since the slit window extended only along the window motor 
portion of the composite chamber, the sole means of studying 
the flow field along the entire length of the cylindrical motor 
section was to exchange the chamber sections so that the win- 
dow motor would at alternate times form either the upstream 
or downstream portion of the chamber. 

Figs. 7a-7e are a strip film record showing the transition 
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FIG. 6 CHAMBER PRESSURE VERSUS PROPELLANT FLOW RATE OF 
4'/.-IN. DIAM ROCKET MOTORS WITH VARIOUS THROAT DIAMETERS 
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FIG. 7b 


Figs. 7a, 7B, 7c 


R\DIATION PATTERN OF TRANSITION FROM SINUSOIDAL TO SHOCK- 

T) PE OSCILLATION. MODIFIED IMPINGING HEAD, 41° CONVERGENT 

NozZLE, 165/s-IN. LONG CHAMBER WITH WINDOW-MOTOR FORMING 

U’STREAM PORTION. SLIT LENGTH, SHOCK FREQUENCY, 
1200 crs 


FIG. 8a FIG. 8b 


Figs. 8a, 8B, 8c 


RADIATION PATTERN OF TRANSITION FROM SINUSOIDAL TO SHOCK- 

TYPE OSCILLATION, MODIFIED IMPINGING HEAD, 41° CONVERGENT 

NOZZLE, 173/s-IN. LONG CHAMBER WITH WINDOW-MOTOR FORMING 

DOWNSTREAM PORTION. SLIT LENGTH, 8?!/3.-IN.; SHOCK FRE- 

QUENCY, 1140 CPS; NOZZLE ENTRANCE | IN. DOWNSTREAM OF 
SLITEND. FILM SPEED ABOUT 40 FT/SEC 


motor, is shown in Fig. 9. The black streaks are drawn in for 
reference purposes. 


from sinusoidal to shock-type oscillations with a 16*/s-in. 


long chamber, the modified impinging head, and a 41-deg con- 

. radis attern o sinusoidal oscillations, as shown 
vergent half-angle nozzle. The window motor formed the attern idalo 

transition period near the nozzle, obtained with the reas- 
ber alone and shown in Figs. 10a, 10b of (1). In both cases 


sembled motor having the window motor forming the down- 
stream portion of the chamber. The nozzle end of the slit 
ended roughly 1 in. upstream of the nozzle entrance. Since it 


one observes pressure waves propagating upstream. The 
path of the pulse is made evident by the sharp deceleration 
which the gas stream undergoes immediately upon passage 


Was necessary to insert a special adapter between the bodies, a Te: : 
the total cylindrical chamber length was 17°/; in. Otherwise of the wave. ’ In addition, the radiation intensity “ f the — 
the assembly and operating conditions in the two tests were wasener aie behind the pulse. No downstream moving pulses 
identical. The chamber pressure was approximately 275 “— discernible. _The transition to shock oscillations nee 
psig. A sketch of the composite flow field, based on the within a short time interval as upstream moving pulses be- 

come more sharply defined and, finally, the first evidence of a 


photographs obtained at the two locations of the window : cme é 
downstream moving pulse can be observed. Within a fraction 


of a second, a fully developed shock-wave radiation pattern 

becomes established. The frequency of the oscillations ob- 

S| Nozze enp tained from the film record and Dynagages agree. For the 
/ 16°/;-in. long motor the sinusoidal oscillation frequency is 
1080-1140 eps, while the shock-type frequency ranges from 
1200-1230 eps. Concurrent with the reflection of each 
shock wave at the injector, a nonluminous region commences, 
which moves downstream at the existing gas velocity. This 
probably represents a cool mixture, consisting of an off-ratio 
fuel-oxidizer composition. It may be noticed that the radia- 
tion intensity emitted by the gas (except for the dark region), 
immediately after passage of the shock, increases considera- 
bly. This is probably caused by the higher pressure and 
temperature existing behind the shock front. In addition, it 
is quite feasible that the shock by its impulsive action ini- 
° Sy é tiates further combustion. The dark nonluminous region, 
on the other hand, seems to undergo no radiation increase as 

e1G. 9 COMPOSITE SKETCH OF SHOCK-WAVE RADIATION PATTERN as it passes through the shock. The upsteam-moving pulse 
SEEN IN WINDOW-MOTOR PHOTOGRAPHS is quite strong, as evidenced by the fact that the gas velocity 
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at times is actually reversed in direction immediately behind 
the shock. 

The shock velocity, within the accuracy of measurement, 
appears constant in the regions bounded by the nonluminous 
S streak. After passage through the nonluminous region the 
See shock has a higher velocity. 
= As described in (1), the tangent of the slope angle enters 
into the calculation of the velocity, and since the angles of the 
shock path approach 90 deg, a comparatively small error in 


angle measurement will cause a large error in the velocity calcu- 
lation. The measurements are probably no better than + 10°;. 
After reflection from the injector, the shock starts downstreain 
at an absolute velocity of about 3900 ft/sec. After passage 
through the nonluminous region, its velocity is approximately 
4700 ft/sec. The absolute shock velocity after returning 
from the nozzle is about 2600 ft/sec. After passing through 
two nonluminous regions, the shock approaches the injector 
at a velocity of approximately 3300 ft/sec. 
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FIG. 10 VELOCITY DISTRIBUTION IN FLOW FIELD AHEAD AND BEHIND SHOCK FRONT 
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The strip film records furnish the velocity distribution of the 
flow field. Assuming nonviscous flow theory, it is possible 
therefore to determine the pressure gradients existing in the 
motor. The partial derivative of the pressure with respect to 
the distance z along the axis of the slit is determined from the 
momentum equation 
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FIG. 11 VELOCITY DISTRIBUTION IN FLOW FIELD AHEAD AND 


BEHIND SHOCK FRONT 
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where 


p = density at any point x at time ¢ 
u = flow velocity at any point x at time ¢ 
p = static pressure at any point z at time / 


The partial derivative of pressure with respect to time is de- 
rived by using the continuity equation 


Por 
where it has been assumed that a p — p relationship exists. 
Since -V Op/Op is the sonic velocity a, Equation [2] can be re- 


written, after substituting the value 0p/dz from Equation [1], 
in the form 


In Figs. 10 and 11 the velocities at selected points of the flow 
field are indicated. In region I (Fig. 9), i.e., in the region be- 
tween the downstream moving shock and the returning 
shock, u (Ou/Ox) is positive, while Ow/Ot is negative. Calcu- 
lating the value of the gradients from the rather crude grid, 
it appears that Ow/Ot and u (Ow/dx) have almost identical mag- 
nitudes. Since 1/p (Op/Ox) is therefore calculated as the dif- 
ference of two large numbers, the possible error of the cal- 
culated value is very large. The calculations seem to indicate 
that 1/p (Op/dx) isso small that it is difficult to say what the 
sign of the pressure gradient is. Since the indicated flow 
velocity is always subsonic, Equation [3] shows that O0p/Ot is 
always negative and therefore the pressure at a fixed point x 
decreases with time. In region II, 0u/ot, Ou/Ox, and u are 
positive. Therefore, Op/dx is negative. Assuming a sonic 
value of a = 3000 ft/sec, it appears that Op/Odt is also nega- 
tive. In region III, 0u/Odt and Ou/dx are greater than zero; 
however, the gas velocity changes from a negative value, 
immediately after the passage of the shock, to a positive value. 
Since the time interval is so short, it becomes very difficult to 
establish an adequate grid to calculate the derivatives. It 
seems, though, that 0p/Ox is positive immediately after the 
passage of the shock, and then becomes negative. Op/Ot 
appears to be negative. Qualitatively, at least, it appears 
that the static pressure of the gas stream behind the downward 
moving shock decreases with time. The upstream moving 
shock is trailed by a rarefaction. 

The shock strength of the pulses is another quantity which 
can be approximated from the velocity gradient across the 
shock. Considering the pressure pulse as a steady one-di- 
mensional normal shock, and neglecting any heat release 
caused by combustion, the pressure ratio p2/p; can be caleu- 
lated at any given point. Using the normal shock equation 


(3) 
y¥-1 a*? 


Y specific heat ratio of combustion gases 
a* critical speed of sound = /UiU:2 
U 
U 


velocity relative to shock of undisturbed gas im- 
mediately before entrance into shock front 

velocity relative to shock of gas immediately after 
passage through shock 


ll 


~: = pressure of undisturbed gas immediately before en- 
trance into shock front 
p. = pressure of gas immediately after passage through 


shock front 
Assuming y = 1.3, it is found that the shock strength p2/p, 
is approximately 1.5. 
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FIG. 12. RADIATION PATTERN OBTAINED WITH ASSEMBLY CON- 

SISTING OF MODIFIED IMPINGING INJECTOR, 10° CONVERGENT 

NOZZLE, AND 17-IN. LONG CHAMBER. WINDOW-MOTOR FORMS UP- 
STREAM PORTION OF CHAMBER. FILM SPEED ABOUT 40 FT/SEC 


Modified Impinging Head, 10-Deg Nozzle, 17-In. 
Long Body Test 


The 41-deg convergent nozzle was replaced with a 10-deg 
convergent, 1.75-in. throat diam nozzle. The operating 
conditions were made similar to those of the 41-deg convergent 
nozzle assembly when the shock-type oscillations occurred. 
The run was smooth. The chamber pressure was then 
lowered about 20 psi and the test repeated. Although the 
run was unstable, no shock-type oscillations were in evidence. 
A photograph of the combustion pattern is shown in Fig. 12. 
The frequency of the oscillations is about 1000 eps. 


Conventional Impinging Head, 41-Deg Nozzle, 
17-In. Long Body 


Pictures of the radiation pattern obtained with the con- 
ventional impinging head and a 10-in. long window assembly 
were presented in (1). Oscillations of 240-340 cps were ob- 
served then during the first oxygen-full aleohol period. 

Five test runs were made with the 17-in. long assembly; 
two with the window motor at the nozzle end, and three with 
the window motor at the injector end of the composite 
chamber. When the chamber pressure was lower than about 
310 psig, the combustion became very unstable throughout 
the run. Severe shock-type oscillations were evidenced, pre- 
ceded by sinusoidal-type oscillations. For the 17%/s-in. as- 
sembly, the sinusoidal frequency was about 800-850 eps while 
the shock type was approximately 1100 eps; while for the 
16°/:-in. assembly the shock-type frequency was about 1150 


Mornin and « 


FIG. 13 RADIATION PATTERN OF SINUSOIDAL OSCILLATIONS OB- 

TAINED WITH ASSEMBLY CONSISTING OF CONVENTIONAL IMPINGING 

HEAD, 41° CONVERGENT NOZZLE, AND 17-IN. LONG CHAMBER. 

WINDOW-MOTOR FORMS UPSTREAM PORTION OF CHAMBER. FILM 
SPEED ABOUT 75 FT/SEC 


Fig. 14a Window-Motor Forming Upstream Portion of 
Chamber 


eps. Figs. 13, 14a, and 14b represent strip film records of the 
sinusoidal and shock-type oscillation, respectively. It may 
be noticed that the radiation pattern is not as uniform as the 
one obtained with the modified impinging head. The non- 
luminous region seems to drift with time. It appears that 
with each succeeding cycle of shock reflection at the injector, 
the nonluminous region seems to commence at a different lo- 
‘ation, and thus there appears to be a phase shift between th: 
motion of the nonluminous region and the shock path unti! 
finally several of the dark streaks almost merge with eac!, 
other to form an extended nonluminous region. Then th: 
cycle repeats itself. Its approximate period is 200-300 cps 


Conventional Impinging Head, 10-Deg Nozzle 
17-In. Long Motor 


All runs with this assembly were unstable. However, thi 
instability was of the sinusoidal type and had a frequency 0 
about 800 eps. The instability was somewhat random, i.e. 
it disappeared for a while and then started up again. A) 
example of the radiation pattern is shown in Fig.15 


21-In. Long Motor Assemblies 


The radiation patterns obtained with the 21-in. assemblie- 
had the general appearance described for the 17-in. as- 
semblies except that the shock oscillation frequency was 
lower (950-1000 cps). However, intermittent shock-type 
oscillations were obtained with the 10-deg nozzle assembly. 


**Stove-Pipe’? Assembly with Conventional Im- 
pinging Head 


Two test runs were made with a 2I-in. long “stove-pipe” 


assembly at chamber pressures of about 125 psig. The 
flow rates at these chamber pressures are equal to those for a 
conventional 1.75-in. throat diam motor operating at a pres- 
sure of 300 psi. Both runs were unstable on ful] stage opera- 
tion. The Dynagage chamber pressure record showed sinu- 
soidal-type oscillations at a frequency of 650 eps. A Dyna- 
gage had also been placed in the alcohol line, about */, ft up- 
stream of the injector inlet. The frequency of oscillations 
agreed with the chamber pressure oscillations; however, the 
alcoho! oscillations were of considerably higher amplitude. A 
portion of the pressure record is reproduced in Fig. 16. The 
strip -film record is shown in Fig. 17. The radiation pattern is 
composed of regions of weak and strong radiation which alter- 
nate at a frequency of about 650 cps. There is no evidence of 
any shock-type oscillation. 


Summary 


The present study seems to confirm that high-frequency in- 
stability in small-diameter motors corresponds to the propa- 
gation of finite pressure disturbances axially along the 
motor. In (1), it was postulated that the unsteady combus- 
tion cycle is initiated by the rapid ignition of a mass of com- 
bustible mixture near the injector head. It was conjectured 


Fig. 14b Window-Motor Forming Downstream Portion of 
Chamber 


RADIATION PATTERN OF SHOCK-TYPE OSCILLATIONS OBTAINED WITH ASSEMBLY CONSISTING OF CONVENTIONAL IMPINGING HEAD, 41° 
CONVERGENT NOZZLE, AND 17-IN. LONG CHAMBER 
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that the lower chamber pressure would increase the quench- 
ine distance, and hence, stable flame holding between the in- 
jector jets would become impossible. The present study ap- 
pears to indicate that these assumptions are incorrect. The 
fact that the initiation of unstable operation is dictated more 
by the head pressure drop than the chamber pressure is a 
peint which definitely invalidates the quenching-distance 
consideration, although it does in no way depreciate the funda- 
mental assumption that the low-velocity gas region existing 
between the liquid fuel and oxidizer jets, in the vicinity of the 
injector holes, is the flame-stabilizing or holding zone of a 
rocket motor. 

It is also pertinent to notice that in case of relatively weak 
oscillations (i.e., sinusoidal oscillations), the upstream-moving 
piuse is always evident while a downstream moving pulse is 
hardly discernible. It almost appears as if the sonic nozzle 
acts not only as a reflector for the pulse but also as an ampli- 
fier. On reaching and reflecting from the injector the pres- 
sure pulse might be expected to influence the mass flow into 
the motor as well as the oxidizer to fuel ratio, if the two pro- 
pellant components have different properties such as density 
aad sonic velocity. As was pointed out in the discussion of 
tie strip-film records, the nonluminous gas region, commenc- 
ing with each reflection of a pulse from the injector, is believed 
to be such an off-ratio, cold-gas region. At the entrance to 
the nozzle, this cold gas has therefore a different Mach num- 
her from the warmer gas which preceded it, since no sharp 
velocity gradient (in time) 0u/O is apparent at the nozzle end 
of the strip film. The possibility of pressure pulse formation 
ciused Ly this unsteady gas Mach number entrance condition 
into the nozzle is a topic which merits further analysis. 

The shock waves are apparently not established instanta- 
neously, as might be expected from a detonation, but build up 
gradually from finite pulses which gain in strength. Whether 
the energy gain is the result of heat addition by combustion or 
a result of amplification caused by the nozzle is a question 
which for the moment cannot be answered. It has, however, 
been found that, generally, shock formation is much more 
prevalent with longer motor chambers than shorter ones. It 
was found impossible to make unstable a 5*/,-in. long water- 
cooled assembly, consisting of the modified impinging head and 
the 41-deg convergent nozzle. With a 10-in. long assembly, 
as reported in (1), only sinusoidal oscillations were obtained. 


end of 


rig. 15 RADIATION PATTERN OBTAINED WITH CONVENTIONAL 

IMPINGING INJECTOR, 10° CONVERGENT NOZZLE AND 17-IN. LONG 

CHAMBER. WINDOW-MOTOR FORMED DOWNSTREAM PORTION OF 
ASSEMBLY. FILM SPEED ABOUT 40 FT/SEC 


The importance of the shape of the convergence section of 
the nozzle on the severity of instability is illustrated by the 
data. This possibility has been recognized previously (4). 
Unfortunately our present equipment did not permit the 
photographie investigation of the phenomenon occurring in 
the nozzle. The reflection of shock waves from solid surfaces 
has only recently been studied intensively, especially by Bleak- 
ney (5) and associates at Princeton University. If the results 
obtained from two-dimensional studies can be used as an 
analogy for three-dimensional studies, it appears that, in 
general, the mode of reflection is a function of the shock 
strength and the angle of incidence, which is defined as the 
angle between the normals to the shock wave and the wall. 
From the data presented by Bleakney and White (6), it ap- 
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FIG. 16 DYNAGAGE PRESSURE RECORD OBTAINED WITH 21-IN. 
LONG “‘STOVE-PIPE’’ MOTOR 


FIG. 17 RADIATION PATTERN OBTAINED WITH 2I1-IN. LONG 

“STOVE-PIPE’’ MOTOR WITH CONVENTIONAL IMPINGING INJECTOR 

HEAD. WINDOW-MOTOR FORMED UPSTREAM PORTION OF CHAMBER 
ASSEMBLY. FILM SPEED ABOUT 40 Fr/SEC 


pears that with the 41-deg nozzle one would probably get a 
“regular’’ reflection for the shock strength obtained in the 
present investigation, while for the 10-deg nozzle it would be 
of the Mach type. The detailed mechanism of the shock re- 
flection at the sonic nozzle, especially under conditions of 
nonuniform flow, is a problem which requires experimental 
investigation. 


Conclusion 


It has been established that high-frequency combustion in- 
stability in small-diameter motors is the result of pressure 
pulse propagation along the length of the motor. The fre- 
quency of oscillation is a function of the cylindrical chamber 
length of the motor. The severity of unstable operation is in- 
creased by low injector head pressure drop, long chamber 
lengths, and large nozzle convergent angle. In the case of 
very large-scale oscillation, the radiation pattern obtained on 
the strip film shows the existence of shock waves. The 
flow field ahead and behind the shock front has been meas- 
ured. 
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DISCUSSION 


By SIN-I CHENG® 


Daniel and Florence Guggenheim Jet Propulsion Center, 
Princeton University, Princeton, N. J. 


Dr. Berman and Mr. Cheney are to be congratulated for 
their important experimental information in confirming the 
basic postulate of the linearized, small-perturbation analysis 
for the stability of combustion system in rockets. 

There has been considerable dispute on the feasibility of 
making linearized stability analysis for the high-frequency 
oscillations encountered in rough burning rocket motors. 
After Berman’s presentation at the ARS 7th Annual Conven- 
tion, a written comment was read in which the shock-type 
oscillation is considered as a definite indication of nonlinearity 
and cannot be analyzed by linearized treatment. It seems 
that, by “rough burning,’’ we mean combustion with large 
amplitude pressure oscillations whose nonlinearity is obvious. 
The existence of the shock-type oscillation does not add any 
new information in this respect. In fact, any measurable or 
observable quantities must be of finite magnitude and cannot 
be treated with linearized analysis in the strict sense. It is, 
therefore, apparent that the properties of the large pressure 
oscillations in rough burning rockets cannot be analyzed on a 
linearized basis. However, this is not what one is most in- 
terested in at the present time. One knows rough burning is 
undesirable and one wants to avoid it. Any knowledge of the 
behavior of the nonlinear oscillations in rockets and how they 
cause damage is, of course, very valuable. But one is more 
interested in finding why such nonlinear oscillations can occur 
in rockets which seem to be designed for smooth operation, 
and how they can develop in rockets that burn smoothly in 
the early stage of operation and then become rough. Once 
one knows why and how, one can seek to avoid or suppress it. 

It does not seem likely that oscillations of large amplitudes 
are introduced into the combustion system or burned gas flow 
system from without. It is, therefore, postulated that os- 
cillations of large amplitudes are the result of rapid amplifica- 
tion of unstable small perturbations of certain frequency 
range which invariably are present in the system as ‘random 
noise.’”’ These random noises are sufficiently small so that the 
stability of these small random disturbances in the system can 
be determined by linearized analysis based on certain excita- 
tion mechanisms. If disturbances in certain frequency range 
become unstable when the mechanisms of excitation over- 
come the mechanisms of damping, they will grow in magni- 
tude with time. When these disturbances are sufficiently 
amplified, new damping or exciting mechanisms of nonlinear 
nature gradually come into play, and finally a definite finite 
magnitude is reached when damping balances excitation. 
Based upon such a postulate as to the origin and the course of 
development of the nonlinear oscillations, the instability of 
small disturbances becomes a necessary condition prior to the 
existence of oscillations of large amplitudes or rough combus- 
tion. The best justification of the contemplated course of 
development is to be able to record experimentally the os- 
cillations during the intermediate stages, and to follow the 
course of amplification that leads to the final stage of large 
amplitude oscillations. This experimental justification is 
found in the present work of Berman and Cheney. In their 


Received December 30, 1952. 
5 Assistant Professor of Aeronautical Engineering. 


96 


experiments it is consistently found that shock-type oscill:- 
tions are always preceded by sinusoidal-type oscillations 
which are of considerably smaller amplitudes. The authors 
are cautious in giving the absolute magnitudes of the two types 
of oscillations, but it is reasonable to believe that the relative 
magnitudes of the two types of oscillations are qualitatively 
correct. There are several radiation patterns of transition, 
but it is rather unfortunate that Berman and Cheney have not 
given more detailed records of the transition processes. .\ 
series of such records following the development more closely , 
such as those given by Schubauer and Skramstad on the am- 
plifieation of pulses in unstable laminar boundary laye), 
would be very impressive. 

It is found that shock-type oscillation has higher frequency 
than the preceding sinusoidal oscillation. This is not sur- 
prising because of its larger propagation speed which result. 
in a smaller characteristic time for the wave to travel tho 
length of the combustion chamber. It does not seem likel. 
that this change of frequency has anything to do with the basi: 
mechanism of instability. 

In a recent theoretical investigation (7, 8)® on the high- 
frequency combustion instability of the intrinsic type in 
liquid propellant rockets, it is found that the ratio of the 
length /s.» of the subsonic portion of the de Laval nozzle to 
the combustion chamber length L is an important parameter 
involved in the boundary condition at the exit of the com- 
bustion chamber. If the ratio /:,,»/L is reduced while othe: 
parameters remain unchanged, the rocket is found to be mor 
unstable. This theoretical deduction is confirmed by Ber- 
man and Cheney’s experimental results, that both larger 
combustion chamber length and larger nozzle converging ang] 
are destabilizing. The theory also predicts that the funda- 
mental mode is the most unstable mode in the series of higher 
frequency oscillations, and that each higher mode is mor 
stable than the preceding lower modes. Experiments pre- 
sented in this paper show that the fundamental and the 
second mode are predominating. The qualitative trends oi 
the theoretica] deduction agree with the experimental results. 
The presence of other modes at finite magnitudes may well be 
due to the nonlinear coupling of oscillations of different fre- 
quencies when the fundamental and the second mode are suf- 
ficiently amplified. 

Berman and Cheney put forth several interesting ideas. 
The nonluminous gas region commencing with each reflection of 
a pulse from the injector end seems to indicate the importance 
of the difference in the responses of the fuel and the oxidizer 
lines to the pressure change at the injector end. Whether 
this fact has fundamental importance in the stability behavior 
of small disturbances in the system seems to be worth further 
investigation. In an effort to explain the “superficially 
stronger” upstream moving waves reflected from the nozzle 
in the radiation streak records, they also propose the excita- 
tion mechanism, that “the sonic nozzle acts not only as a re- 
flector but also as an amplifier.”’ Since no combustion is sup- 
posed to take place in the nozzle, and the steady flow in the 
nozzle is considered to be isoenergetic, the pressure wave can- 
not possibly gain energy in the course of its reflection from 
the nozzle. The reflected wave must not be stronger than 
the incoming wave in order to conform to the energy princi- 
ple. In fact, normally incident compression shock is re- 
flected as expansion from open end and as weaker shock from 
a closed wall. The three-dimensional reflection is much more 
complicated but will not create energy in the course of reflec- 
tion. The fact that the reflected wave moving in upstream 
direction appears to be stronger in the radiation streak is 
probably due to a simple geometrical reason which becomes 
more significant when the angle between the radiation streak 
and the time axis becomes larger. If this angle is initially 

6 Numbers in parentheses refer to the References at end of this 
discussion. 
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Technical Notes 


Comment on “An Approximate Theory 
of Porous, Sweat, or Film Cooling With 
Reactive Fluids” by L. Crocco’ 


ANDREW E. ABRAMSON? 


National Advisory Committee for Aeronautics 
Lewis Flight Propulsion Laboratory, Cleveland, Ohio 


- px paper by L. Crocco presents a very interesting ap- 
- proach to a rather complicated mechanism. Evaluation 
of the validity of the simplifying assumptions used in the 
analysis remains for future work; however, the assumptions 
aj/pear to be in line with conventional assumptions used for 
analysis of this type. 

Whereas the paper is applicable to various types of cooling, 
it was reviewed only with reference to liquid film cooling, as 
data for this type of cooling are available from experiments 
conducted at the NACA Lewis Flight Propulsion Laboratory. 
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Coolant flow per circumferential length, |b/sec/ft 
(a) 4-inch diameter smooth-surface tube 


FIG. 1. VARIATION OF LIQUID-COOLED LENGTH WITH WATER FLOW 
PER UNIT CIRCUMFERENCE 


Gas Gas mass 

temperature velocity Reynolds 

Legend (°F) (Ib/ (see) (sg ft)) number 

& 900 57.9 8.18 X 105 

ie) 900 70.2 9.90 

—_—— 900 42.5 6.00) NACA 

—_-— 1400 45.3 5.25} RM 
4.70) E50F19 


—--— 1600 43.6 
—— Theoretical results corresponding to experi- 
mental conditions 


1 JOURNAL OF THE AMERICAN ROCKET Society, vol. 22, No- 
vember-December 1952, pp. 331-338. 
2 Aeronautical Research Scientist. 


The assumption of a smooth liquid surface limits the range 
of liquid coolant flow rates to which the results of this analy- 
sis can be applied. The results of a visual investigation of 
annular liquid flow with concurrent gas flow, NACA RM 
E51013 (now unclassified), indicate that a smooth liquid sur- 
face is encountered at only the relatively low liquid flow rates. 
It is questionable whether liquor or coolant flow-rates of this 
magnitude are practical for film cooling application to rocket 
engines. 

The expression used for the evaluation of the frictiort coeffi- 
cient (i.e., Cry = 0.046 Re~.*) does not appear to be valid for 
the flow system encountered in liquid film cooling. Investi- 
gations of the pressure drop encountered in a two-phase two- 
component flow reported by Lockhardt and Martinelli (in 
Chemical Engineering Progress, January 1949) and by Bergelin 
(in Chemical Engineering, May 1949) indicate the friction 
coefficient would be larger than given by the previous relation. 

A comparison of the theoretical results for the inert coolant 
with experimental results obtained from a gasoline-air burner 
test rig is given in Fig. 1. The experimental results are taken 
in part from NACA RM E50F19 (now unclassified), and the 
remainder from data to be published in the near future. The 
theory predicts a longer liquid cooled length than observed 
experimentally, and gives the best agreement at the low 
coolant flow rates which correspond to a smooth liquid sur- 
face. In applying the theory, the term hy — ho, which appears 
in Equation [18] of Croeco’s paper, was evaluated as the 
change in enthalpy of the coolant vapor in going from its 
boiling point to the free stream temperature. The paper does 
not indicate clearly whether the term hy, — ho is referred to the 
main gas stream or the coolant vapor. The above interpreta- 
tion seems the most logical and gives a closer agreement with 
the experimental results. 
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GROVE 
SELF-ACTUATED 


CONTROLS AND REGULATORS 
for 


High pressure liquids and gases 


Pressure Reducing 
Back Pressure 
Relief and By-Pass 
Manual Valves 


Special light weight units 
for air-borne service. 


Special materials for corrosive 
or low-temperature applications. 


GROVE CONTROLS, INC. 
6529 Hollis Street 


Emeryville California 


WANTED 


ENGINEERS and 
SCIENTISTS 


Unusual opportunities for outstanding and 
experienced men. 


These top positions involve preliminary and produc- 
tion design in advanced military aircraft and special 
weapons, including guided missiles. 


IMMEDIATE POSITIONS INCLUDE: 


% Boundary layer research scientists . . . 
Electronic project engineers... 
Electronic instrumentation engineers... 
% Radar engineers... 
% Flight-test engineers... 
% Thermodynamicists ... 
% Servo-mechanists . . . 
% Electro-mechanical designers . .. 
% Electrical installation designers . .. 
Excellent location in Southern California. Gener- 
ous allowance for travel expenses. 
Write today for complete information on these essen- 
tial, long-term positions. Please include résumé of 


your experience & training. Address inquiry to Di- 
rector of Engineering, 


NORTHROP AIRCRAFT, INC. 


1041 E. Broadway, Hawthorne (Los Angeles County), Cal. 


Combustion Studies in Rocket Motors 


(Continued from page 96) 


60 deg, a given velocity change of the particle that could turn 
the velocity vector in the radiation record for 15 deg in tiie 
upstream direction can only turn the velocity vector for 7'/; 
deg in the downstream direction. While an upstream shock 
will almost cause flow reversal, a downstream shock of equal 
strength can hardly turn the velocity vector through (5 
deg. In extreme cases, if the initial angle is as big as §8)- 
deg, a shock which is strong enough to cause flow reversal if 
moving upstream may hardly be discernible when it is moviiig 
downstream. 
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Jet Propulsion News 


C. F. WARNER, Purdue University, Associate Editor 
with the assistance of W. G. BOHL 


Rockets 


COMPLETELY redesigned and re-engineered Viking 
high-altitude rocket, one of a series built by The Glenn 
L. Martin Company for the Naval Research Laboratory, 
equalled the previous single-stage altitude mark of 135.6 miles 
wien launched at the White Sands Proving Ground, Las 
Cruces, N. Mex. Shorter by some six feet and of slightly 
greater diameter, with completely new control, fuel, and elec- 
tronics systems, Viking No. 9’s most readily distinguishable 
external feature was smaller, triangular fins. The power 
plant, made by Reaction Motors, Inc., was of the same type 
used on previous Vikings, liquid oxygen and alcohol, with a 
20,000-Ib thrust. Fig. 1 shows a Viking just after take-off. 
White cloud at right is steam from water in pit under launch- 
ing stand, vaporized by tremendous heat from rocket engine. 
Up to the time of this firing, the altitude record for single- 
stage rockets was held by Viking No. 7, which was also 
launched at WhiteSands on August 7, 1951. A WAC Corporal 
launched from the nose of a German V-2 in a two-stage opera- 
tion has reached about 250 miles. 

Vikings are being used by NRL scientists to extend knowl- 
edge of (1) the physical state of the earth’s atmosphere to as 
great a height as possible, the temperatures, pressures, densi- 
ties and composition being principal objects of study; (2) 
the fundamental nature and properties of the ionosphere and 
those processes which lead to its formation; (3) solar and 
terrestrial radiation; and (4) the physics of high-energy parti- 
cles by high-altitude cosmic ray studies. This information is 
expected to be of practical benefit to radio communication, 
meteorology, and the guided missile program, as well as to the 
basic physical sciences. Four additional rockets have been 
ordered by NRL from The Glenn L. Martin Company. 


+ 


SOME details of the rockets currently being developed by 
Oerlikon in Switzerland have recently been revealed. In ad- 
dition to their first powder rockets, which have a diameter of 
3.15 in. and a length of 24 in., they are producing several ver- 
sions of a 2.2-lb contact detonating rocket. All of these 
solid rockets may be fired from fighter and bomber aircraft, or 
from ground-stationed antiaircraft units. 

A longer range rocket with a proximity fuse is also under 
development. A small turbine driven by the exhausted 
rocket gases powers an electric generator. Current from the 
generator supplies power to a miniature radar instrument 
which continuously measures the distance of the rocket from 
the target and ignites the warhead at the correct instant. 

One of the Oerlikon guided missiles, guided to the target 
by the guide-beam method, has four small wings and is pow- 
ered by a nitric acid motor. These antiaircraft missiles have 
a length of 16.4 ft, a combustion cut-off speed of 2460 fps, 
and are said to reach an altitude of 66,000 ft. The gross 
weight of the missile, including a 44-lb war head, is 550 Ib. 


Aircraft 


AN Air Force North American F-86D Sabrejet interceptor 
powered by a GE turbojet was recently flown to a new official 


(Photograph courtesy of The Glenn L. Martin Company 


FIG. 1 VIKING NO. 9 RISING FROM LAUNCHING STAND 


world’s speed record of nearly 700 mph by Air Foree Capt. J. 
Slade Nash. The F-86Dcombineshigh speed with all-seeing ra- 
darand new electronic equipment to enable it to operate in any 
kind of weather. The “all-weather” engine makes possible 
the successful completion of missions under severe icing con- 
ditions. ‘Icing up” has been a serious problem in jet opera- 
tion. Supercooled water droplets, striking the air inlet at the 
nose of the engine, form into ice crystals and, within seconds, 
can block passage of air into the engine. This problem has 
been overcome by giving the engines “hot noses.’”’ Air, at 
temperatures high enough to melt any ice crystals, is bled 
from the compressor and piped to passages in the nose. 


+ + + 


THE Republie Aviation Corporation’s XF-91, powered by 
a GE J-47 turbojet and a Reaction Motors rocket engine, re- 
cently flew faster than sound in level flight. Thus the XF-91 
became the first U. 8. plane designed for warfare to achieve a 
level flight speed of greater than Mach 1. The necessary ad 
ditional thrust was supplied by the rocket motor. 

Republic also recently delivered its first production model 
of the F-84F Thunderstreak fighter bomber to the Air Force. 
The Thunderstreak, capable of speeds over 700 mph, is 
powered by the Curtiss-Wright built British Sapphire engine. 
The F-84F is said to be able to carry a small atom bomb to a 
range of more than 2000 miles. 


++ + 


THE Swedish Air Force claims that its new two-seated 
fighter, the Saab-32, powered by a Rolls-Royce Avon, has a 
top speed of 700 mph. The Saab-32 is armed with rocket pro- 


Epiror’s Nore: The information reported in this Section has been selected from approved news releases originating with the 
Department of Defense, private manufacturers, universities, etc., and from published news accounts in journals and newspapers. The 
reports are considered generally reliable, although no attempt has been made to verify them in detail. 
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jectiles, bombs, and cannon. The fighter has a span of 41 ft 
and an over-all length of about 46 ft. 


+ 


THE new Handley Page crescent-wing jet bomber, the 
H.P. 80, powered by four Sapphire turbojet engines, has been 
ordered into quantity production for the British Bomber 
Command. No details as to performance or dimensions are 
available. However, the maker claims that no other bomber 
can fly as far or as fast with its bomb load. 


~+ + 


THE world’s first delta-wing seaplane, the XF2Y-1 “Sea 
Dart,” a jet-powered experimental aircraft built for the Navy 
by Consolidated Vultee Aircraft Corporation, is shown in Fig. 
2 as it began preflight trials on San Diego Bay. Two Westing- 
house J-34 turbojet engines power the present Sea Dart, but it 
is thought that two J-46’s will be used eventually. It is re- 
ported that ten additional Sea Darts have been ordered at a 
reported cost of $2,250,000 each. 

To provide better rough-water landing and take-off charac- 
teristics, the XF2Y-1 is equipped with retractable hydro-skis. 
This marks the first application of the hydro-ski to a combat 
type aircraft in this country and, according to Navy officials, 
may well be the first such application in the world. The new 
hydro-ski will enable high-speed aircraft to operate from 
water, grass, ice, and snow. The NACA, Navy, Air Force, 
the Research Division of All-American Airways, and the Edo 
Corporation developed the hydro-ski in a six-year joint re- 
search program. The Canadian bush pilots originated the 
idea when they used their winter snow skis to land on water 
during periods of melting ice. The watertight fuselage of the 


FIG. 2 
CONVAIR “SEA DART” 
WATER-BASED INTERCEPTOR 


(Photograph courtesy of Convair) 


Convair XF2Y-1 enables it to float on the water. However, 
the hydro-skis take all of the water impact shock during land- 
ing and take-off. At take-off, the extended hydro-skis are 
under water as the XF2Y-1 begins its run (Fig. 3). The 
plane’s cambered surface provides hydrodynamic lift and 
lifts the fuselage out of the water; then the lower V-shaped 
surfaces of the skis provide the lift. The final take-off run is 
made on the skis. When landing, the extended skis take tlie 
initial water impact and carry the plane as it loses speed until 
it sinks and floats on its fuselage. 

The use of the hydro-ski will enable high-speed coastal de- 
fense aircraft to operate without the necessity of constructiiig 
expensive runways. When fitted with small integral whee's, 
the hydro-skis will enable an aircraft to land upon the water 
and taxi up a beach. This type of landing gear would be of 
great value to aircraft operating from an invasion beach hea:|. 


Turbojet Engines 


SWEDEN’S newest turbojet, the 9000-lb thrust Dover:,, 
will soon go into production at the underground plant «f 
Svenska Flaygmotor AB near Trollhéttan. The Dovern h:s 
undergone extensive flight tests in a modified Avro Lancaste «. 
This new turbojet was developed from the 3300-lb thru-t 
Skuten. Other more powerful engines are currently under 
development by Svenska Turbin-fabriks AB Ljungstréia 
(STAL) of Finspong. 


+ 

THE Société National d’Etude et de Construction de 
Moteurs d’Aviation (SNECMA) of Paris has been flight- 
testing a “jet deviator” which turns the turbojet engine thrust 


FIG. 3 
“SEA DART” 
AT START OF 
TAKE-OFF RUN 


(Photograph courtesy of Convair) 
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through an angle of approximately 120 degrees toward the 
front of an aircraft. This forward-directed thrust effects a 
braking action which permits landing in much shorter dis- 
tances. Although design and performance data are secret, 
SNECMA claims that they hope to be able to land a jet plane 
without the use of wheel brakes. It is reported that the de- 
viator consists of from eight to ten small jet nozzles mounted 
in « circle in the engine nacelle. A mechanical device deflects 
a portion of the engine exhaust gases to these nozzles which 
turns the gases forward. So far the deviator has only been 
tested on a Vampire. 


+ 
TOP officials of the British deHavilland Engine Company, 
Ltd., recently inspected the General Electric Company’s jet- 
engine manufacturing facilities at Evendale, Ohio, in con- 
junction with a “full interchange” of aircraft engine experi- 
ence and knowledge agreement between the two companies. 
Tlie agreement between deHavilland and the General Electric 
Company provides for an exchange of both current and future 
information in the aircraft gas-turbine field with the excep- 
tion of aircraft nuclear propulsion. The deHavilland Com- 
puny has developed such engines as the Goblin, in the 3000- 
3500 Ib of thrust category, and the Ghost, rated at 5000 lb, 
which powers the jet-propelled transport, the Comet. 


Jet Assisted Take-Off 


IT IS reported that the Flying Tiger Air Line is seeking 
C.A.A. approval to install Marboré II turbojets on all of its 
C-46 aircraft. Flying Tigers carried out the first American 
tests of the engine as a standby unit to provide extra 
thrust on take-off and to serve as an emergency power plant in 
the event of failure of one of the regular engines. It is pro- 
duced in the United States by Continental Aviation and En- 
gineering Corporation, under license from the developers, 
Société Turbomeca of Bordes, France. 

The engine, 82 in. long and 25 in. in diam, is mounted under 
the fuselage between the wheels in a streamlined fairing. It 
delivers 880-lb thrust at sea level, and at 5000 ft it provides 
the C-46 with a 100 per cent power increase for climb. The 
engine uses the same fuel as the regular piston engine. 


Surface Cooling at High Speeds 


A STUDY of the air-resistance heating effect encountered 
by supersonic aircraft and missiles has been under way for 
some time at the NACA Ames Aeronautical Laboratory. 
Major contributions to the knowledge of aerodynamic heating 
are expected from tests of the Bell X-2 stainless-steel super- 
sonic aircraft. The tests of the X-2 will provide a check on 
the theoretical and scale model work. 

Aerodynamic heating, directly related to the boundary 
layer surrounding the aircraft, increases with speed. It has 
been estimated that this boundary layer may reach a tem- 
perature of 500 F for a plane flying 2000 mph at altitude, and 
1200 F at 3300 mph. The boundary-layer friction may also 
account for approximately one half of the aircraft’s drag at 
Mach 3. During a flight of short duration, the heat generated 
may not have time to penetrate beyond the surface skin of the 
aircraft or missile. The German V-2 is an example of high- 
speed short-duration flight. Aircraft designed for prolonged 
high-speed flight willneed some type of cooling. Water cooling 
may be practicable in some cases. The water can be evapo- 
rated at the hot skin of the missile or aircraft and vented over- 
hoard. 

It has recently been discovered that the regions of laminar 
flow in the boundary layer can be increased by cooling. Lami- 
nar flow reduces both skin friction and aerodynamic heat- 
ing. Experimental work at Ames is directed toward extend- 
ing the regions or areas of laminar flow at supersonic speeds. 
The amount of cooling required for laminar flow is of critical 
importance. At the present time, both the skin friction and 
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FIG. 4 INSTALLATION OF G. E. TAIL CONE THERMOCOUPLES 


the heat transfer rates have been investigated by the NACA 
at speeds up to Mach 3 for both laminar and turbulent flows. 


Thermocouple Controls 


A VITAL component of automatic electronic controls for 
F-86D jet interceptor planes is being manufactured by the 
General Electric Company’s Meter and Instrument Depart- 
ment. Thermocouples, mounted in the tail cone of the J47- 
GE-17 jet engine, transmit temperature changes to the elec- 
tronic control. The control system automatically regulates 
fuel flow to provide optimum performance of the jet power 
plant under all flight conditions. 

Sustained accurate measurement of the jet-engine exhaust 
gas temperature previously had been prevented by mechani- 
cal difficulties caused by the terrific buffeting of the gases 
which often approached the speed of sound. This problem 
was overcome by G-E engineers after five years of research 
and testing. The new thermocouple responds rapidly to tem- 
perature changes, is highly accurate, and retains its character- 
istics in spite of severe operating conditions. The couples 
(Fig. 4) operate over a temperature range, 1000 to 1900 F. 


New Facilities 


PURDUE University has begun construction of the Univer- 
sity-financed Jet Propulsion Laboratory No. 2. The Labora- 
tory, when completed, will contain three rocket-motor re- 
search cells and control rooms, a classroom, drafting room, in- 
strument laboratory, and a hydraulic laboratory. In addi- 
tion, the new building will house the administrative offices of 
the jet propulsion personnel, and a technical library. A 
wing, to be added to the building in the future, will house air- 
burning jet propulsion engine laboratories. 

At present, jet propulsion research at Purdue is being con- 
ducted in Rocket Laboratory No. 1 (housing four rocket- 
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’ motor test cells, control rooms, and a machine shop) and in 
the Gas Turbine Laboratory. The research projects are be- 
ing sponsored by ONR, Squid, and USAF. 


+ + 


REACTION Motors, Inc., in co-operation with the Navy 
Department, announced that it is undertaking a $3,500,000 to 
$4,000,000 construction program to provide new engineering, 
research, laboratory, and administrative facilities. The new 
facilities, which will include a 60,000-sq ft machine shop, will 
be financed in part by the Company and in part by a govern- 
ment facilities contract. 

Reaction Motors’ new construction program continues a 
story of rapid growth and achievement for this rocket con- 
cern. RMI, the oldest rocket company in the United States, 
last year celebrated its tenth anniversary. RMI rocket en- 
gine applications include the famous Air Force (Bell) X-1, 
first supersonic airplane; the Navy (Douglas) D-558-II Sky- 
rocket, present world speed (1238 mph) and altitude (79,- 
494 ft) record holder for piloted aircraft; the Air Force (Re- 
public) XF-91, turbojet and rocket-powered interceptor air- 
plane which is the first combat-ready aircraft to travel at 
supersonic speeds; the Navy (Fairchild) Lark subsonic missile; 
and the Navy (Martin) Viking high-altitude sounding rocket, 
present world speed (4100 mph) and altitude (136 miles) 
record holder for single-stage missiles. 

Mr. Charles W. Newhall, Jr., executive vice-president of 
the Company, stated that plans for Reaction Motors’ new 
facilities call for work to be started this year with a completion 
date some time in late 1953. The new plant will be in the 
Rockaway-Lake Denmark vicinity. It will provide com- 
plete, modern buildings, laboratories, and equipment for all 
phases of the development and fabrication of rocket engines 
and related components, and for a variety of other aircraft 
and guided missile auxiliaries. In addition to Reaction Mo- 
tors’ new facilities, substantial additions and improvements 
are being made to the present test facilities leased and 
operated by the Company at the Naval Air Rocket Test 
Station, Lake Denmark, Dover, N. J. 


++ + 


IN A new fuel test laboratory just completed at Northrop, 
engineers can stand on the ground and observe fuel mecha- 
nisms of jet aircraft operating under conditions that in flight 
would whisk them in a matter of minutes from searing desert 
sea-level temperatures to the thin, cold atmosphere found at 
80,000 ft. Northrop built the fuel test facility because, in the 
current era of high-speed aircraft and guided missiles, it is 
easier and more economical to “bring” the weather into the 
laboratory than it is for engineers to take their elaborate test- 
ing equipment to high altitude. 

Housing several full-scale fuel system mockups, including 
one of the U. 8. Air Force’s “grand slam’ Scorpion F-89D 
rocket-armed all-weather interceptor, this facility provides 
reproduction of atmospheric conditions at any altitude up to 
80,000 ft above sea level within temperatures ranging from 
—85 degrees to +170 degrees. 


~++ + 


TESTING equipment ten miles up in the sky without ever 
leaving the ground is now possible in this new turbopump 
altitude test chamber at General Electric’s Aircraft Gas Tur- 
bine Division facilities here. This chamber is one of the 
many testing facilities used to test equipment under the sever- 
est conditions which could be encountered in actual flight. 
Simulated altitudes over 60,000 ft can be obtained, as well as 
simulated aircraft climb rates in excess of 20,000 fpm and 
climb rate angles as high as 60 degrees. Used primarily for 
testing reheat (afterburner) fuel pumps for jet engines under 
varied conditions and altitudes, this chamber has a volume of 
approximately 950 cu ft and can accommodate a 470-gal 
capacity fuel tank. Engineers can observe pump operation 
under all conditions of altitude by watching through ob- 
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servation ports provided in the sides of the chamber. The 
entire end of the chamber can be removed to permit the in- 
stallation of larger equipment to be tested under altitude con- 
conditions. 

+ + 

A NEW altitude test chamber—largest in the United 
States, and probably largest in the world—has recently been 
completed at the McDonnell Aircraft Corporation in *t. 
Louis, Mo. This addition to the company’s expanding e:- 
gineering facilities in the physical laboratory is designed io 
simulate altitudes up to 70,000 ft and has a temperature ranye 
of —100 F to +165 F. Relative humidity can be con- 
trolled from 15 per cent to 95 per centat dew point temperatures 
above 35 F. The unit will provide McDonnell enginee's 
with the exact answers to many of the functional problenis 
faced when aircraft equipment is operated under conditions of 
extreme heat, cold, altitude, or humidity. Despite the initi:! 
cost of the installation, considerable saving may be effected 
by “de-bugging” all equipment well ahead of the actual 
flight-test program. 

~+ + + 

THE world’s two most powerful motors and two smaller 
“giants” to drive the world’s largest rotating machine are 
nearing completion at the East Pittsburgh plant of the Wesi- 
inghouse Electric Corporation. The four motors—two rate: 
at 83,000 hp each, and the smaller pair rated at 25,000 hp 
each—will be installed tandem fashion in a rotating machine 
almost two-football-fields long to be used in new wind tunnels 
of the U.S. Air Force. 

In addition to the 216,000-hp four-motor drive, the gargan- 
tuan rotating machine will have five tremendous compressors, 
which are under construction at the Westinghouse plant in 
Sunnyvale, Calif. The entire unit is slated for installation in 
transonic and supersonic wind tunnels being built at the 
Arnold Engineering Development Center in Tullahoma, Tenn. 
The Center is under the command and supervision of the Air 
Research and Development Command, responsible for the 
over-all U.S. Air Force research and development program. 

Prior to this order from the Air Force, the most powerful 
motors were those built for the Grand Coulee Dam by Westing- 
house and rated at 65,000 hp each. The 83,000-hp giants 
each will stand 21'/2 ft high and weigh 225 tons. Their 122- 
ton rotors will spin at the dizzy rate of 600 rpm. Some 31 
miles of copper wire are used in their windings. From end to 
end it will measure about 500 ft, and the four-motor drive 
alone will weigh more than 1000 tons. 

~ + + + 

The Aluminum Company of America has been awarded « 
contract to install a 20,000-ton extrusion press under a lease 
agreement with the Air Force. The 20,000-ton press will be 
installed in the Alcoa Lafayette, Ind., plant where a 14,000- 
ton press is currently being installed. These presses will be 
used to extrude ribbed shapes for aircraft wing sections. The 
internal strains will be relieved and the shapes will be straight- 
ened by a 3,000,000-lb stretcher capable of straightening 
shapes having a cross section of 60 sq in. The larger press 
will be capable of extruding shapes weighing 2!/, tons. 


Theodore Von Karman to Lecture at 
Princeton 


The Department of Aeronautical Engineering of Princeton 
University announces that Professor Theodore von Karman, 
world-famous aerodynamicist, isscheduled to deliver three guest 
lectures on the subject, ““Aerothermodynamics and Combus- 
tion Theory,” on Tuesday, Wednesday, and Friday, March 24. 
25, and 27, 1953, at 8:00 p.m., in Frick Auditorium. The 
Daniel and Florence Guggenheim Jet Propulsion Center at 
Princeton is sponsoring the lectures. A general invitation to 
attend has been extended to those outside of Princeton Uni- 
versity who are interested in scientific aspects of cumbustion 
phenomena. 
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ocket Society News 


H. K. WILGUS, Associate Editor 


Final Report of the Ad Hoc Space Flight Committee of 
the American Rocket Society 


HIS Committee was authorized and 

constituted by the ARS Board of 
Directors at the meeting of February 4, 
1952. The basic task of the Committee 
was briefly outlined in the Minutes of the 
meeting as follows: 

The Ad Hoc Space Flight Committee is 
to investigate and recommend to the 
Board the ways in which the American 
Rocket Society can best serve its members 
an fulfill the aims outlined in its By- 
Laws for activity in the field of space 
flight, satellite space stations, astronautics, 
etc., and to define the extent to which the 
American Rocket Society should become 
active in these fields.”’ 

\ better understanding of the Commit- 
tec’s comments and recommendations 
requires a brief restatement of certain 
background information. 

Qur Society, as originally incorporated 
in April 1931, under the name of the 
“American Interplanetary Society,”’ was 
empowered to engage in the following ac- 
tivities: 

“The promotion of interest in and ex- 
perimentation toward interplanetary ex- 
peditions and travel; the mutual enlight- 
enment of its members in matters bearing 
on the astronautical, physical and other 
problems pertinent to man’s ultimate con- 
quest of space; the stimulation by expen- 
diture of funds and otherwise of American 
scientists toward a solution of the problems 
which at present bar the way toward tra- 
vel among the planets; the raising of 
funds for research and experimentation; 
and such other activities as the Society 
may from time to time deem necessary or 
valuable in connection with the general 
aim of hastening the day when interplane- 
tary travel shall become a reality.”’ 

In 1934 the present name of the Society 
was adopted, and in 1947 the charter pow- 
ers were changed to read as follows: 

“The purpose of the Society shall be to 
aid and encourage by all suitable means 
the development and application of the 
principle of jet propulsion as applied to 
rockets, aircraft, water and underwater 
craft and to all other appropriate and 
practical devices; to aid and encourage 
the development of the sciences and en- 
gineering techniques pertaining thereto, 
and to create increasingly wide interest in 
the field of jet propulsion and rocketry 
among both technicians and laymen, to 
the end that jet propulsion in all its vari- 
ous forms shall rapidly and permanently 
he developed for the ultimate good of man; 
the preparation, collection, correlation 
and dissemination, by publication or other- 
wise, of facts, information, articles, books, 
pamphlets and other literature pertaining 
to jet propulsion, rocketry and subjects 
relating thereto; the establishment of a 
library containing such literature for the 
information of members, scientists and 
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others to whom the privileges of such li- 
brary may be granted by the Society; the 
collaboration or affiliation with other or- 
ganizations, whether technical or other- 
wise, in any manner, and to any extent 
which, in the judgment of the Society or 
of the Board of Directors, will best aid in 
accomplishing its objectives; the raising 
of funds for research and experimentation, 
and such other activities as the Society or 
the Board of Directors may from time to 
time deem necessary or desirable in con- 
nection with the foregoing.” 

Throughout its deliberations the Com- 
mittee has borne in mind the basic pur- 
poses and powers of the Society. 

The great interest in the subject of space 
flight, which undoubtedly motivated the 
creation of the Society, was necessarily 
subordinated to the practical problems on 
hand attending the preparation for World 
War II, the conduct of the war, and the 
post var developments in a world still filled 
with turmoil. A real interest in the sub- 
ject of space flight, however, has always 
persisted among many members of the 
Society. This interest has gained mo- 
mentum during the past three years, and 
has been demonstrated by the attendance 
of Society delegates at the International 
Astronautical Congresses held in London 
(1951) and Stuttgart (1952); the adher- 
ence of the Society to the International 
Astronautical Federation; the direct ex- 
pressions of many members of the Society, 
such as the remarks of Commander R. C. 
Truax during the 1951 Honor’s Night 
Dinner, and so on. The Committee was 
constituted on a special and temporary 
basis, as a result of this interest. 

The work of the Committee was initi- 
ated by each member submitting proposed 
agendas for study and discussion. The 
agendas were then disseminated on a ran- 
dom basis to a third of the membership of 
the Society for comment and viewpoints. 
The replies evidenced uniform interest in 
the subject of space flight among the mem- 
bership of the Society. 

A formal meeting of the Committee was 
held in Washington, D. C., on May 17, 
1952, and the proceedings were recorded 
verbatim. Many informal meetings were 
held by groups of Committee members, 
and numerous consultations were held 
with leading scientists. A second formal 
meeting was held in Washington, D. C., 
on November 21, 1952. A third and final 
formal meeting was held at noon on Thurs- 
day, December 4, 1952, to consider com- 
ments of the Society membership on this 
report. 


FINDINGS 


With respect to technical and other sub- 
stantive matters, the Committee finds as 
follows: 

(a) Space flight is a basic challenge to 


the progress of mankind and its achieve- 
ment will ultimately prove to be worth the 
cost and effort required; but it is not pos- 
sible to spell out in hard economic terms 
at the present time the extent of the phys- 
ical and social returns that will flow from 
the conquest of space. We have insuffi- 
cient information at the present time on 
which to base considered judgment con- 
cerning the utility or feasibility of various 
widely publicized proposals advocating 
space ships as military weapons, as re- 
connaissance or bombing platforms, as 
missile guidance adjuncts, as_ television 
relays, as weather observation stations, for 
basic physical experiments, or for astro- 
physical observations. 

(b) Successful flight from the earth 
and return of satellite, lunar, and inter- 
planetary space ships will some day be 
technically feasible. Judgment is with- 
held on the technological feasibility of 
establishing large, inhabited, satellite 
stations. No opinion is offered as to the 
date when successful space flight might be 
achieved. The rate of progress will be 
determined to a great extent by the 
amount and quality of effort that is focused 
on the project, and perhaps to an even 
larger degree by unforeseen fundamental 
problems that may arise as the program 
proceeds. It is believed, however, that 
space flight of inhabited vehicles may be 
achieved during this century if mankind is 
in a position diligently to undertake the 
task, 

(c) The development of a rocket-pro- 
pelled space ship requires the successful 
solution of many basic problems involved 
in the development of long-range guided 
missiles. In this respect, a space flight 
program cannot be divorced from the cur- 
rent military program for guided missiles. 
A space flight project, if initiated under 
present world conditions, would require the 
efforts of many of the very engineers and 
scientists now engaged in the develop- 
ment of guided missiles, and consequently 
must be carefully evaluated in that light. 

(d) For the purpose of building up ex- 
perience and knowledge in an orderly fash- 
ion, the following projects should be em- 
phasized and supported, in the sequence 
indicated: High-altitude sounding rockets 
for scientific purposes; ground-to-ground 
unmanned guided vehicles; orbital (satel- 
lite) unmanned vehicles; manned satellite 
vehicles; manned, returnable, lunar, or 
interplanetary vehicles. The develop- 
ment of manned, long-range, ground-to- 
ground vehicles similar to present rocket- 
propelled aircraft is not considered a pre- 
requisite for achieving interplanetary flight, 
although data obtained from such proj- 
ects could serve as a useful adjunct to the 
technology of guided rocket vehicles. 
The development of a one-way rocket to 
the moon is considered a waste of effort, 
although interesting from a publicity and 
educational standpoint. Encouragement 
and support should be given to a high- 
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altitude rocket research program, such as 
that carried out by the Army, Navy, and 
Air Force, and co-ordinated by the Upper 
Atmosphere Research Panel. 

(e) Key factors in the achievement of 
space flight are recognized as money, 
competent personnel, a recognized long- 
range program, and continued public en- 
thusiasm and support. Key problems, as 
understood today, in the establishment of 
an unmanned satellite or space vehicle are 
(not in order of importance): Reliability 
of a vehicle’s mechanical and electronic 
systems, radiation hazards, aerodynamic 
heating upon re-entry into the earth’s 
atmosphere, automatic solution of naviga- 
tion problems, effects of prolonged weight- 
lessness upon humans, recovery of spent 
stages, communications, telemetering, and 
minimum weight designs of all compo- 
nents. 

RECOMMENDATIONS 

With regard to the position of the So- 
ciety and measures that it might take to 
promote the early achievement of space 
flight, the Committee recommends the 
following: 

(a) The Society should on appropriate 
occasions take a public stand in favor of 
realistic research and development looking 
toward the achievement of space flight. 
This stand should be based upon the set of 
beliefs outlined in the findings above. It 
should refrain from lending support to 
sensational projects motivated by a desire 
for publicity or private gain. 

(b) The Society should implement its 
stand by: 

Encouraging the presentation of high 
quality technical papers on space flight 
at National and Sectional meetings of the 
Society. 

Publication of contributed technical 
articles on all aspects of space flight in the 
JOURNAL. 

Encouraging scientific and technical ef- 
fort on space flight problems on the part of 
any group competent to perform such 
work. 

Co-operating with other nonprofit or- 
ganizations having sound programs along 
the lines indicated in the Findings above. 

Advocating the expenditure of govern- 
ment funds, if no adequate private financ- 
ing is available, to support work on the 
problems of space flight, to the extent that 
such effort is compatible with national 
security. 

Bringing such matters to the attention 
of the public through various media in a 
manner calculated to promote sound pub- 
lic thinking on the subject. 

Fostering international co-operation and 
exchange of ideas on the subject of space 
flight through the medium of the Inter- 
national Astronautical Federation, United 
Nations Educational, Scientific and Cul- 
tural Organization, and _ International 
Council of Scientific Unions. 

Maintaining a continuous program of 
co-operation with existing foreign rocket 
societies and encouragement of formation 
of such societies in countries where they do 
not at present exist. 

Establishing a permanent Society space- 
flight committee which will meet at regular 
intervals for the purpose of recommending 
to the National Board of Directors meth- 
ods of implementing the aforesaid objec- 
tives effectively. 
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(c) The Society should point out the 
need for a realistic survey of the technical 
problems standing in the way of space 
flight. The support of appropriate govern- 
mental agencies should be solicited in 
order to consider the extent of these prob- 
lems. A survey of this kind may serve 
to stimulate research in the required di- 
rections. At this stage of man’s knowl- 
edge the Committee believes that a study 
program of all phases of space flight, sup- 
ported by the government, possibly by the 
National Science Foundation, would pro- 
vide a sound basis for a vehicle construc- 
tion program in the least possible time. 
Such a study would tend to ‘‘clear the air’ 
of some of the sensational concepts en- 
gendered by unprincipled or uninformed 
fiction writers. Further, such a study 
would elicit support of competent pro- 
fessional engineers and scientists who at 
present are fearful of supporting projects 
which are incompletely thought out. At 
the same time, recognition of excellent 
work by United States and foreign experts 
would inspire even greater efforts. 

(d) The Society should not presume to 
conduct practical research or development 
under its own auspices on the subject of 
space flight, but rather to function as a 
“catalyst” for such activity on the part of 
private industry, universities, private re- 
search institutions, and governmental 
agencies. 

(e) The Society should refrain from en- 
dorsing any specific development proposal 
in preference to other competitive propos- 
als, but should offer a forum at its meetings 
and in its JouRNAL for the open discussion 
of each proposal. For the present, the 
Society should take the position that to 
decide upon a specific design is less im- 
portant than the promotion of broad tech- 


nical analysis of all the problems in-— 


volved. Similarly, the Society should not 
be associated with activities conducted pri- 
marily for the purpose of publicity, such 
as science fiction movies, television or 
radio shows of a fictional character, and 
the like, unless such efforts represent a 
sincere attempt to educate the public. 
The cost and difficulty inherent in the 
achievement of space flight must never be 
underplayed. 

(f) The Society should promote the 
view that, to the extent that world politi- 
cal conditions permit, and with due regard 
to United States security, the conquest of 
space should be carried out on an inter- 
national basis. The Society should re- 
iterate its adherence to the aims of the In- 
ternational Astronautical Federation as a 
specific step in this direction. 

Finally, the Committee is of the view 
that exploratory space flight as a specula- 
tive venture into the unknown should be 
regarded as the ultimate objective of a 
reasonable effort on the part of the United 
States and mankind as a whole. It is a 
logical function of the Society and its con- 
stituted committee to recommend ways of 
synthesizing existing unconnected activi- 
ties with the objective of directing such 
efforts toward a unified space flight pro- 
gram. 

Respectfully submitted: Andrew G. Ha- 
ley, Chairman; F. C. Durant, III, R. W. 
Porter, Milton W. Rosen, Kurt R. Stehl- 
ing, Franki van der Wal, John R. Young- 
quist. November 25, 1952. 


Separate Report by R. W. Porter 


While generally in agreement with the 
other members of the Committee on the 
foregoing report of November 25, I would 
like to see the conclusions more specifica! |v 
worded. For example, I believe that the 
following statement of policy should he 
adopted by the Society: “The Americ: 
Rocket Society recommends that explora- 
tory space flight as a speculative venture 
into the unknown should be regarded as 
the ultimate objective of a reasonable «- 
fort on the part of the United States of 
America, It further recommends th:t 
actual design or construction of satellite 
or space vehicles should not be under- 
taken until the present war emergency is 
resolved, and that such work should !:e 
regarded as one of several large projecis 
into which our national effort could |e 
directed when no longer urgently needed 
for military work.”’ 

In furtherance of the above policy, the 
following should be adopted as_ specific 
objectives of the Society. 

a Toseek and publish articles of sound 
scientific content on the subject of sate|- 
lite and space flight. 

b To encourage discussion of satellite 
and space flight in Section and Nation:l 
meetings of the Society. 

c To foster international co-operation 
and exchange of ideas on the subject of 
space flight through the medium of the 
IAF and UNESCO. 

d To carry on a continuing campaign 
of encouragement and support for the 
high-altitude rocket research program be- 
ing carried out by the Army, Navy, and 
Air Force, co-ordinated by the Upper At- 
mosphere Research Panel. In particular, 
to interest this Panel in problems related 
to space flight. 

e To urge the National Science Foun- 
dation to give serious consideration at 
this time to the possibility of establishing 
within its organization the initial phase of 
a long-range space-flight project, and if 
successful, to assist the Foundation in ob- 
taining financial support by all means at 
the Society’s disposal. 

f To establish a permanent  space- 
flight committee which will meet at reg- 
ular intervals for the purpose of imple- 
menting the above objectives effectively. 

The Society should not sponsor any in- 
dividual projects except as indicated 
above, nor should it seek association at this 
time with activities conducted primarily 
for the purpose of publicity, such as 
science fiction movies, TV or radio 
shows of a fictional character, and the like. 


Admiral Bolster Addresses 
Joint Meeting at North- 
eastern New York Section 


Or JANUARY 13, 1953, a meeting 
sponsored jointly by the ARS North- 
eastern New York Section and the 
Schenectady Chapter of the ASME, was 
held at Union College, Schenectady, N. Y. 
Over 200 persons attended to hear the 
guest speaker, Admiral Calvin M. Bolster, 
Chief of Naval Research. 
Admiral Bolster, in discussing ‘‘Aspects 
of Current Rocket Research,’ described 
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AT THE JOINT MEETING OF ARS NORTHEASTERN NEW YORK SECTION AND ASME SCHENECTADY 
CHAPTER, JANUARY 13, 1953 


Left to right: 
Schenectady, N. Y.; K. 


G. M. Ketchum, chairman, ASME Schenectady Chapter; 
Admiral Calvin M. Bolster, Chief of Naval Research, who was guest 


Board of Directors; 


Commdr. F. W. Kinsley, USN, Office of Inspector of Naval Machinery, 
Berman, vice-president, ARS Northeastern New York Section; 


R. W. Porter, ARS Nationat 


speaker; and F. M. Cooper, president, Northeastern New York Section. 


the workings of the Bureau of Naval 
Research, and recounted the latest de- 
velopments in jet-assisted take-offs. His 
talk was illustrated with films of ex- 
periments with jet-assisted aircraft, and 
concluded with movies of a recent Viking 
firing. 


Admiral Bolster, an active member of 
the ARS as well as of the Institute of the 
Aeronautical Sciences, is a noted aviator 
in both heavier and lighter-than-air craft. 
In 1947 he was chosen to the post of 
Deputy Chief of Naval Research, and was 
made Chief in 1951. 


Joint ARS-[AS Rocket Propulsion Session 


THE American Rocket Society, in co- 

operation with the Institute of the 
Aeronautical Sciences, held a_ technical 
session on rocket propulsion on the opening 
day of the IAS Annual Meeting, January 
26-29, 1953, at the Hotel Astor, New 
York, N.Y. Mr. F. C. Durant, III, presi- 
dent of the American Rocket Society, 
served as chairman of the session, which 
was attended by more than 130 persons. 
Three technical papers constituted the 
program. 

The first paper, “The Effect of Fluid 
Properties on the Spray Formed by Two 
Impinging Jets,” by R. J. Priem and M. F. 
Heidmann, both of NACA Lewis Flight 
Propulsion Laboratory, was presented by 
Mr. Priem. The paper described the four 
different spray shapes which were ob- 
tained using glycerol-water mixtures and 
hydrocarbon oils, and the conditions under 
which these sprays were formed, together 
with their apparent properties. A hypothe- 
sis was offered on the connection be- 
tween spray fluctuations and combustion 
instability in rocket motors. Excellent 
high-speed photographs of the spray 
patterns accompanied the talk. 

The second paper, “Earth Scanning 
Techniques for Orbital Rocket Vehicles,” 
by Kurt R. Stehling of Bell Aircraft 
Corporation, examined two possible meth- 
ods of utilizing an earth satellite vehicle 
to gather and transmit information about 
the earth’s surface. One of these methods 
utilized an optical objective lens image- 
converter system; the other, a microwave 
scanner. Assumed were a vehicle flight 
altitude of 500 miles, orbital velocity of 5 
miles per sec, a hypothetical target of 600- 
yd diam, and restricted design require- 
ments of a minimum weight satellite 
vehicle payload. 
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It was shown that the net receiving 
power of the optical-video system is 
limited by the image converter and availa- 
ble transmission band width, and not by 
the optical train. Effects of the atmos- 
phere, meteorological conditions, and 
albedo were also discussed. The approxi- 
mate resolving power, weight, and power 
and space requirements of a microwave 
radar scanning unit were outlined, and the 
advantages and disadvantages of this 
system were weighed. 

Mr. Stehling believed that, with the 
present state of radar seanning (which 
can be discussed on an unclassified status), 
the optical system has the weight and 
resolving power advantage, but its day- 
light dependence and atmospheric and 
cloud attenuation limitations are severe 


enough to indicate questionable feasility. 
Furthermore, the low resolving power of a 
radar unit, from the 500-mile altitude, and 
the attenuation losses tend to limit the 
system’s usefulness. 

The final paper, ‘‘Measurements in 
Rocket Engine Testing,’ was presented 
by Howland B. Jones, Jr., of the Potter 
Aeronautical Company. The — paper 
covered the application of conventional 
and special instrumentation equipment to 
the static testing of liquid propellant 
rocket engines. In his introduction, Mr. 
Jones discussed peculiarities of this testing 
field, including accuracy limits and need 
for transient as well as static measure- 
ments. The conventional quantities of 
pressure, flow, force, temperature, linear 
and angular displacement, and vibrational 
acceleration were then treated. Modi- 
fications to conventional measuring equip- 
ment and specially designed instruments 
were then described. Of particular in- 
terest were the descriptions of a multi- 
channe! pressure recording unit operating 
into a single cathode-ray tube, a combina- 
tion of a 4-gun cathode-ray oscillograph 
and tape recorder, and the use of a con- 
ventional sonic analyzer with associated 
equipment. 


Washington- Baltimore 
Section Hears Rosen 


ORE than 150 members and guests 

attended the Washington-Baltimore 
Section meeting on February 4, 1953, 
held at the Enoch Pratt Library, Balti- 
more, Md. Mr. Milton W. Rosen, Viking 
Project Chief, Naval Research Labora- 
tory, was the guest speaker. 

Mr. Rosen presented a talk on “High 
Altitude Sounding Rockets,” which was 
highlighted by the showing of slides. 
Following his speech, he showed sound 
movies of the development and firing of the 
Viking rocket. 

In the discussion after his speech, when 
asked his opinion of recent proposals to 
develop a satellite vehicle, Mr. Rosen’s 
reply was to the effect that we should 
crawl before we walk; that there are too 
many unknown and uninvestigated factors 
to warrant immediately going overboard 
with large expenditures. 


AT THE JOINT ARS-IAS ROCKET PROPULSION SESSION DURING THE IAS ANNUAL 


MEETING, JANUARY 26, 1953 


Left to right: F. C. Durant, III, ARS president, with the three engineers whose 
papers were presented: R. J. Priem, NACA Lewis Flight Propulsion Laboratory; 
Kurt R. Stehling, Bell Aircraft Corporation; and Howland B. Jones, Jr., Potter 


Aeronautical Company. 
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Technical Literature Digest 


H. S. SEIFERT, California Institute of Technology, Associate Editor 


Contrisutors: D. Altman F. C. Gunther 
D. I. Baker H. S. Seifert 
R. B. Canright F. H. Wright 


Book Reviews 


DEVELOPMENT OF THE GUIDED MISSILE, 
by K. W. Gatland, Philosophical 
Library, Inc., New York, N. Y., 1952, 
133 pp. $3.75. 


Reviewed by G. P. Surron 
North American Aviation, Inc. 


This book is a major revision of a series 
of articles published in the British maga- 
zine Flight in 1951. It contains separate 
chapters on antiaircraft (or antimissile) 
missiles, on air-to-air missiles, on super- 
sonic rocket flight, on high-altitude re- 
search rockets, on satellite and space ve- 
hicles, and on interplanetary flight. The 
book gives a good discussion of many of the 
technical features necessary for these 
specialized applications of rocket power 
and as such is helpful to acquaint the 
reader with typical technical and develop- 
mental problems and in outline form also 
with existing or contemplated solutions. 
For this purpose this book is recom- 
mended. 

Most of the examples used by the author 
to illustrate the various missile types and 
their technical features are taken from 
German World War II practice. Data 
and information on U. S. guided missiles 
are next most prominent, but there is 
relatively little additional mention of 
British missiles and other rocket devices. 

The very fact that the author has not 
been able to secure much information on 
British rockets and missiles (a fact about 
which he seems to be unhappy) is a tribute 
to the British security system. In sharp 
contrast, it is amazing to this reviewer 
how much information the author has 
been able to obtain on American develop- 
ments. 

The reader should be careful not to in- 
terpret the information of this book as 
being a picture of the actual state of the 
art of rocketry and guided missiles in the 
various countries. Some of the detailed 
technical information on various missiles is 
not correct, and the emphasis given to 
specific missiles is not properly weighted. 
Some of the missiles which are treated in 
considerable detail in this book are rela- 
tively unimportant test missiles or merely 
missile studies, whereas others, which 
are mentioned only in passing, are im- 
portant tactical weapons. 

The last two chapters on space and satel- 


lite rockets and interplanetary flight 
differ in tenor from the rest of the book 
because they deal not with existing pieces 
of equipment but only with proposed 
plans and future thinking. 


Comets AND METEOR Streams, by J. G. 
Porter, John Wiley & Sons, Inc., New 
York, N. Y., 1952, 123 pp. $5.25. 


Reviewed by 8. B. NicHoLson 
California Institute of Technology 


In this delightful little book, Dr. Porter 
tells about the celestial mechanics of 
bodies moving in the solar system subject 
to the gravitational pull of the sun and 
the planets. The problems of comets 
and meteor swarms are specifically dis- 
cussed, but the general laws of motion are 
first developed. The presentation is 
simple and direct, technical enough to be 
rigidly correct but readable even for those 
who lack the mathematics to follow all 
the formulas. Although most of the de- 
velopments follow the classical lines, some 
original departures and extensions are 
found. 

Those who expect a book in the Inter- 
national Astrophysics Series to deal mainly 
with the spectra of comets may be sur- 
prised to find only one page of the 123 
devoted to that subject. In addition to 
the orbital problems, the statistics of 
comets and meteors are discussed in de- 
tail. These statistics have direct ap- 
plication to the theories of the origin of 
comets and of their place in the solar 
system. 


BALLISTICS IN THE SEVENTEENTH CEN- 
tury, by A. R. Hall, Cambridge Univ. 
Press, N. Y., 1952, 186 pp. $4. 


Reviewed by H. J. Srewart 
California Institute of Technology 
Jet Propulsion Laboratory 


The author had a twofold purpose in 
writing this book. The first was to de- 
scribe the foundations of the science of 
exterior ballistics, and the second was to 
relate the development of this science to 
the state of technological development, 
of practical military science, and of the 
governmental organization with respect to 
research and development existing during 
the 17th century. In _ particular he 
investigated the suggestion that exterior 
ballistics was the first applied science 
whose development was stimulated by 


military requirements and whose resu'ts 
were applied in increasing the effective- 
ness of artillery. 

The author first presents the back- 
ground of the relations between govein- 
ment and industry and of the techniques of 
manufacturing heavy ordnance. He then 
describes the practical aspects of the «rt 
of gunnery in the 17th century. The 
principal portions of the book trace tiie 
development of the science of mechanics 
with particular emphasis on_ ballistics 
from the Aristotleian period, through the 
‘impetus’”’ theories, through the Galilean 
period, and to the Newtonian period 
during which the theory of ballistics 
was placed on a basis very close to the 
currently accepted theory. The conclu- 
sions are reached that it is inappropriate 
to consider the classical development of 
ballistics as an example of applied science, 
that the scientific interest ballistics 
was motivated primarily by the fact that 
the motion of a body through the air was 
one of the classical problems treated by 
Aristotle and was thus the most suitable 
test of any new theory, and that, even 
when the theory was placed on a firm 
foundation, the state of technological 
development was too primitive to permit 
effective applications. The idea that the 
theory should be useful was widespread 
but, in the 17th century, premature. 

This book should be of general interest 
to any engineer or scientist and of partic- 
ular interest to students of the history 
of science and to anvone concerned with 
the problems of military research and 
development. 


ENGINEERING MATERIALS, THEIR Me- 
CHANICAL PROPERTIES AND APPLICA- 
Tions, by J. Marin, Prentice-Hall, 
Inc., New York, N. Y., 1952, 491 pp. 
$8.70. 


Reviewed by C. E. LEvVor 
California Institute of Technology 
Jet Propulsion Laboratory 


The stated objective of the book is to 
provide the student with sufficient know!- 
edge of the behavior of stressed engineering 
materials to permit their intelligent selec- 
tion and use. The book is divided into 
three parts: Part one deals with mechani- 
cal properties, their definition, deter- 
mination, and utilization. Part two is a 
brief treatment of certain important mz- 
terials, emphasizing specific mechanical 


Eprtor’s Note: This collection of references is not intended to be comprehensive, but is rather a selection of the most significant and 


stimulating papers which have come to the attention of the contributors. 


The readers will understand that a considerable body of 


literature is unavailable because of security restrictions. We invite contributions to this department of references which have not come 
to our attention, as well as comment on how the department may better serve its function of providing leads to the jet propulsion applica- 


tions of many diverse fields of knowledge. 
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properties and the factors influencing 
these properties. Part three presents a 
description of materials-testing machines 
and strain gages. 

The first part of the book, which is in- 
tended as an undergraduate text, deals 
with general mechanical properties and 
coutains most of the basic information re- 
quired for an understanding of the strength 
of materials. It treats several phases of 
the subject not ordinarily covered in a 
tet of this type. For example, static 
properties under combined stress, fatigue 
properties, impact properties, and creep 
an | temperature properties are thoroughly 
covered. 

The section of the book covering specific 
properties of materials attempts to cover 
too broad a subject in a few short chapters 
to give really adequate information for 
selection and use of materials. It is 
probable that, by listing only typical or 
average mechanical properties of a few 
m:iterials, certain misconceptions may 
arise. An example of a_ possible mis- 
conception may be cited. On page 47 
of the book are listed strength-weight 
ratios of four materials: an aluminum 
alloy, a magnesium alloy, a steel alloy, 
and a glass laminated plastic. The 
strength-weight ratios for these materials 
as listed in the text show the aluminum 
alloy test (double that of the steel alloy), 
the magnesium alloy second, the laminated 
plastic third, and the steel a poor fourth. 
The design stresses selected for these ma- 
terials are, however, inconsistent. Actu- 
ally, if one inquires deeply enough into the 
strength of materials, it will be found that 
the strength-weight ratio is a function of 
the state of development of the group of 
materials. Compare, for example, a com- 
monly used steel alloy SAE 4130 heat- 
treated to 200,000 ultimate tensile strength 
aluminum alloy 75S heat-treated to 77,000 
UTS, and magnesium alloy AZ80A heat- 
treated to 47,000 UTS. The strength- 
weight ratio for each of the materials is 
200,000/0.284 = 7.05 X 10°, 77,000) 
0.101 = 7.62 X 105, and 47,000/0.065 = 
7.24 X 105, respectively, in units of lb per 
in.? divided by Ib per in.*. From this it is 
seen that the strength-weight ratio of most 
commonly used, high-strength materials is 
essentially the same. More often than 
not, other factors than strength-weight 
ratio govern selection of materials. 

In general, it is considered that Part two 
of the book treats the subject in a manner 
not ideally suited for engineering design 
use. It would seem preferable to intro- 
duce the student to the need for certain 
information and then to supply tables of 
basic data rather than a few typical 
examples. 

The bulk of Part one of the book can be 
recommended to rocket engineers. The 
remainder, although understandably writ- 
ten and containing much information, 
seems rather sketchy to be of use for de- 
sign purposes, 


Books 


Proceedings of the London Conference 
on Heat Transfer, Institution of Mechani- 
cal Engineers and American Society of 
Mechanical Engineers, 1952, 496 pp. 
£2 5s. 

Instrument Engineering, by 
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Draper et al., McGraw-Hill Book Co., 
Inc., N. Y., 1952, 269 pp. $6. 

The Design and Analysis of Experi- 
ments, by O. Kempthorne, John Wiley & 
Sons, Ine., N. Y., 1952, 631 pp. $8.50. 

Numerical Mathematical Analysis, by 
J. B. Scarborough, Johns Hopkins Press, 
Baltimore, 1950, 511 pp. $6. 


Rocket Engines 


Tip-Pod Rockets for F-89D Confirmed, 
unsigned, Western Aviation, vol. 32, 
Sept. 1952, pp. 12-13. 

Bell X-5; Variable Wing Sweep, un- 
signed, /ntervias, vol. VII, no. 8, 1952, 
p. 446. 

Rockets Push XF-91 Past Mach 1, un- 
signed, Aviation Week, vol. 57, Dee. 15, 
1952, p. 17. 

NACA’s Rockets Gather Data for 
Missile Design, by J. A. Shortal, SAF J., 
vol. 60, Nov. 1952, pp. 17-19. 

Russians Push Rocket Fighter Program, 
unsigned, Amer. Aviation, vol. 15, Nov. 10, 
1952, pp. 2-4. 

Snarler Boosts the Jets, unsigned, 
Aviation Age, vol. 17, Nov. 1952, p. 44. 

Jet Assisted Take-Off, by von 
Karman, /nteravia, vol. VII, no. 7, 1952, 
pp. 376-379. 


Jet Propulsion Engines 


The Turbojet, Ram-Jet and Rocket Air 
Force of the Future, by N. Macmillan, 
Aeronautics, vol. 27, Aug. 1952, pp. 24-31. 

Tests on a Working Model Ram Jet ina 
Supersonic Wind Tunnel, by J. R. Sing- 
ham, F. W. Pruden, and R. C. Tomlinson, 
Gt. Brit. Aero. Res. Counc., Reports and 
Memoranda 2568, 1951 (Nov. 1947): 
Brit. Info. Serv., N. Y., 16 pp., illus. 

Ramjet Favored for Mach 2-4 Range, 
unsigned, Aviation Week, vol. 57, Dee. 15, 
1952, pp. 21-28. 

A Theoretical and Exnerimental In- 
vestigation of the Feasibility of the Inter- 
mittent Ram-Jet Engine, by A. Kahsne 
et USN-Squid/TR-35, Princeton, 
Aug. 1, 1952. 

Pulse-Jet Engine Has 200-Hour Operst- 
ing Life, unsigned, Mach. Design, vol. 24, 
Oct. 1952, pp. 250-253. 

Can Pulsejets Simplify Army Logistics, 
by W. Coughlin, Aviation Week, Oct. 13, 
1952, pp. 21-28. 

The Shock Ignition Engine, by F. F. 
Rand, Jr., Aero. Engng. Rev., Oct. 1952, 
pp. 22-27, 62. 

Afterburners and Variable Nozzles, by 
R. Kress, Aereplane, vol. 82, Nov. 14, 
1952, pp. 665-666. 

Military Role for Ram-Jet Hiller, un- 
signed, Flight, vol. 37, Sept. 1952, p. 19. 

Principles and Application of Jet Power 
to Aircraft, by I. B. Lakowitz, Amer. 
Helicopter, vol. 24, Oct. 1952, pp. 6-7. 

Simulation and Linearized Dynamics of 
Gas Turbine Engines, by J. R. Ketchum 
and R. T. Craig, Nov. 1952, NACA TN 
2826. 

Protecting Aprons from Jet-Blast, un- 
signed, Flight, vol. 61, Oct. 3, 1952, p. 446. 


Heat Transfer and Fluid 
Flow 


A Comparative Study of Momentum 
Transfer, Heat Transfer, and Vapor 


Transfer. Part III. Free Convection, 
by C. 8. Yih, Colorado A & M College, 
Report No. 3, Feb. 1951, 32 pp. 

Heat-Transfer and Flow-Friction Char- 
acteristics of Some Compact Heat- 
Exchanger Surfaces. III, by A. L. 
London, W. M. Kays, and D. W. Johnson, 
Trans. ASME, vol. 74, Nov. 1952, pp. 
1167-1178. 

An Experimental Investigation of Con- 
vective Heat Transfer to Air from a Flat 
Plate with a Stepwise Discontinuous 
Surface Temperature, by 8S. Scesa and 
F. M. Sauer, Trans. ASME, vol. 74, 
Nov. 1952, pp. 1251-1255. 

Convective Heat-Transfer and Flow- 
Friction Behavior of Small Cylindrical 
Tubes: Circular and Rectangular Cross 
Sections, by W. M. Kays and A. L. Lon- 
don, Trans. ASME, vol. 74, Nov. 1952, 
pp. 1179-1189. 

Additional Measurements of Heat Con- 
ductivity of Nitrogen, Carbon Dioxide, 
and Mixtures, by F. G. Keys, Trans. 
ASME, vol. 74, Nov. 1952, pp. 1303-1306. 

Design of Apparatus for Determining 
Heat Transfer and Frictional Pressure 
Drop of Nitrie Acid Flowing Through a 
Heated Tube, by B. A. Reese and R. W. 
Graham, NACA RM 52D08, June 16, 
1952, 18 refs., illus., 61 pp. 

An Experiment on Nozzle Flow, by 
V. D. Naylor, Aircraft Engng., vol. 24, 
Nov. 1952, pp. 344-347. 

Convection Heat Transfer Coefficients 
in de Laval Nozzles, by S. J. Kline, 
Stanford Univ. ONR Contract N6-onr- 
251 T/OVI. TR HS-2, 1950. 

Heat Transfer in a Nozzle at Supersonic 
Speeds, by O. A. Saunders, Engineering. 
vol. 173, Aug. 29, 1952, pp. 281-284. 

On the Theory of the Ranque-Hilsch 
Cooling Effect, by J. J. van Deemter, 
Appl. sci. Res. (A), vol. 3, 1952, pp. 174- 
196. 

On the Temperature Distribution Along 
a Semi-Infinite Sweat-Cooled Plate, by N. 
Ness, J. Aero. Sci., vol. 19, 1952, pp. 
760-768. 

On Heat Transfer over a Sweat-Cooled 
Surface in Laminar Compressible Flow 
with a Pressure Gradient, by M. Mordu- 
chow, J. Aero. Sct., vol. 19, Oct. 1952, 
pp. 705-712. 

The Measurement of Skin Friction in a 
Turbulent Boundary Layer at a Mach 
Number of 2-5, Including the Effect of a 
Shock-Wave, by W. F. Cope, Proc. 
Roy. Soc. Lond. (A), vol. 215, Nov. 6, 
1952, pp. 84-99. 

A Free-Flight Investigation of the 
Possibility of High Reynolds Number 
Supersonic Laminar Boundary Layers, by 
J. Sternberg, J. Aero. Sci., vol. 19, Nov. 
1952, pp. 721-733. 


Combustion 


Activation Energies and Frequency 
Factors in the Thermal Decomposition of 
Paraffin Hydrocarbons, by M. G. Peard, 
Proc. Roy. Soc. Lend. (A), vol. 214, Oct. 9, 
1952, p. 471. 

Structural Influences in the Oxidation of 
Aliphatic Amines, by C. F. Cullis, Trans. 
Faraday Soc., Nov. 1952, p. 1023. 

The Oxidation, Decomposition, Ig- 
nition, and Detonation of Fuel Vapors and 
Gases, by R. O. King, Canad. J. Technol., 
vol. 30, Sept. 1952, pp. 222-257. 

The Water Gas Shift Reaction in Fuel 
Systems. I. The Water Gas Shift on 
Carbon, by D. G. Ingles, Trans. Faraday 
oa vol. 48, Part 8, Aug. 1952, pp. 706- 

12. 
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Pressure Limits of Flame Propagation 
of Pure Hydrocarbon Air Mixtures at 
Reduced Pressure, by A. E. Spakowski, 
NACA RM E52H15, 1952, 35 pp. 

Flame Stability of Liquid-Vapor Oxy- 
gen Mixture, by J. A. Browning and M. 
L. Thorpe, Project Squid TM DART-1, 
Feb. 11, 1952. 

The Temperature of the Hydrogen- 
Fluorine Flame, by R. H. Wilson, Jr., 
J. B. Conway, A. Engelbrecht, and A. V. 
Grosse, J. Amer. Chem. Soc., vol. 73, 1951, 
p. 5514. 

On the Stability of Laminar Flame 
Fronts, by M. Lessen, J. Aero. Sci., vol. 19, 
Dec. 1952, p. 852. 

Combustion Tunnel 
Second Quarterly Report, by H. W. 
Emmons, Harvard University, Combus- 
tion Tunnel Lab. QR-2, May-Aug. 1952. 

Investigation of Chemical Combustion 
Kinetics; Final Report, by R. A. Jermain, 
J. A. Pisani, and I. Rosen, Reaction 
Motors, Inc., Report RMI-442F, Nov. 17, 
1952. 

Spark Ignition; Effect of Molecular 
Structure, by H. F. Calcote, C. A. Gregory 
Jr., C. M. Barnett, and R. B. Gilmer, 
Indust. Engng. Chem., vol. 44, no. 11, 
Nov. 1952. 

Sur le passage par la vitesse du son au 
cours de la détente d’un gaz en com- 
bustion, by M. Serruys, Comptes rendus 
des séances de l’académie des sciences, vol. 
235, no. 21, Nov. 24, 1952, pp. 1280-1281. 

On the Solution of the Poisson-Boltz- 
mann Equation with Application to the 
Theory of Thermal Explosions, by P. L. 
Chambré, J. Chem. Phys., vol. 20, no. 11, 
Nov. 1952, pp. 1795-1797. 

The Effect of High-Frequency Sound 
Waves on an Air-Propane Flame, by 
C. J. Kippenhan and H. O. Croft, Trans. 
ASME, vol. 74, Nov. 1952, pp. 1151-1155. 

Study of Laws of Combustion of Col- 
loidal Powders (in French), by H. Muraour 
and G. Aunis, Mém. Artill. fr., vol.25, no. 1, 
1951, pp. 117-165. 


Laboratory ; 


Propellants, and 


Materials 


Solution to the Caking Problem of 
Ammonium Nitrate and Ammonium Ni- 
trate Explosives, by J. Whetstone, /ndust. 
Engng. Chem., vol. 44, Nov. 1952, pp. 
2663-2667. 

Kinetics of the Transformations of the 
Polymorphous Modifications of Ammo- 
nium Nitrate. II. Effect of the Previous 
Treatment on the Rate of Transformation 
NH,NO;IV~+NH,NO,III, by B. V. Ero- 
feev and N. I. Mitskevich, Zhur. Fiz. 
Khim., vol. 26, 1952, pp. 848-861. 

Behavior of Liquid Hydrocarbons with 
White Fuming Nitric Acid, by C. H. Trent 
and M. J. Zucrow, Indust. Engng. Chem., 
vol. 44, Nov. 1952, pp. 2668-2673. 

Determination of the Pore-Size Distri- 
oution in Porous Metal Bodies, by P. R. 
Marshall, Nature, vol. 170, Nov. 1952, 
p. 761. 

The Metallurgy of Uranium, by E. W. 
Colbeck, Metal Indus., vol. 81, Nov. 14, 
1952, pp. 387-389. 

A Half Century of Aluminum Develop- 
ments, by J. D. Edwards, J. Electrochem. 
Soc., vol. 99, Nov. 1952, pp. 298C-300. 

New Titanium-Boron Alloy Steel Shows 
Promise for Jets and Rockets, by J. L. 
Everhart, Materials and Methods, vol. 36, 
Sept. 1952, pp. 97-98 

Titanium; the New Light Metal, 
unsigned, Light Metals, vol. 15, Sept. 1952, 
pp. 292-293; Oct. 1952, pp. 318-319. 
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Fuels, 


Physical-Chemical 
Topics 


Calculation of Vibrational Relaxation 
Times in Gases, by R. N. Schwartz, 
Z. I. Slawsky, and K. F. Hersfeld, J. 
Chem. Phys., vol. 20, Oct. 1952, pp. 
1591-1599. 

Ionization Produced in Metallic Salts in 
Flames. III. Ionic Equilibria in Hydro- 
gen/Air Flames Containing Alkali Metal 
Salts, by H. Smith and T. M. Sugden, 
Proc. Roy. Soc. Lond. (A), vol. 211, 1952, 
pp. 31-58. 

Integrated Absorption for Vibration- 
Rotation Bands of COs, by D. Weber, R. 
J. Holm, and 8. 8S. Penner, J. Chem. Phys., 
vol. 20, Nov. 1952, p. 1820. 

Emissivity for CO. at Elevated Pres- 
sures, by R. J. Holm, D. Weber, and 8. 8. 
Penner, J. Appl. Phys., vol. 23, Nov. 1952, 
p. 1283. 

Infrared Emission and Absorption of 
Thermally Excited Carbon Dioxide, by 
R. H. Tourin, J. Chem. Phys., vol. 20, 
1952, p. 1651. 

Vibration Spectrum of Nitric Acid in the 
Infrared, by C. Fréjacques, Comptes 
rendus des séances de l’académie des sciences, 
vol. 234, 1952, pp. 1769-1770. 

Kinetics of Decomposition of Nitric 
Oxide at Elevated Temperatures. II. 
The Effect of Reaction Products and the 
Mechanism of Decomposition, by H. 
Wise and M. F. Frech, J. Chem. Phys., 
vol. 20, Nov. 1952, pp. 1724-1727. 

Flame Emission Spectrum of HBr, by 
B. W. Bullock et al., J. Chem. Phys., 
vol. 20, Nov. 1952, pp. 1808-1809. 

Separation of the Absorption Spectra of 
NO, and N2O, in the Range of 2400-5000A, 
by T. C. Hall, Jr., and F. E. Blacet, 
J. Chem. Phys., vol. 20, Nov. 1952, pp. 
1745-1749. 

Polytropic Exponent of Isentropic Ex- 
pansions in the Co-existence Region, by 
S. G. Reed, Jr., J. Chem. Phys., vol. 20, 
1952, pp. 533-534. 

An Equation of State Applicable to 
Gases at Densities Near That of the Solid 
and Temperatures Far Above the Critical, 
by T. L. Cottrell and S. Paterson, Proc. 
Roy. Soc. Lond. (A), vol. 213, 1952, pp. 
214-225. 

Thermodynamic Properties of Nitrogen 
as Function of Pressure and Temperature 
Between 0 and 6000 Atmospheres and 
—125 and +150°, by R. J. Lunbeck, A. 
Michels, and G. J. Wolkers, Appl. Sci. Res., 
vol. A3, 1952, pp. 197-210. 

The Condensation of Nitrogen in a 
Hypersonic Nozzle, by W. W. Willmarth, 
J. Appl. Phys., vol. 23, Oct. 1952, pp. 
1089-1095. 

Viscosity, Density, and Critical Con- 
stants of OF,, by R. Anderson et al., 
oo Chem., vol. 56, 1952, pp. 473- 


Vapor Pressure of Tetranitromethane, 
by G. Edwards, Trans. Faraday Soc., 
vol. 48, 1952, pp. 513-515. 

Vapor Pressure-Composition Measure- 
ments on Aqueous Hydrazine Solutions, by 
J.G. Burtle, Indust. Engng. Chem., vol. 44, 
1952, pp. 1675-1686. 

Vapor Pressures and Freezing Points of 
the System Nitrogen Tetroxide-Nitric 
Oxide, by A. G. Whittaker et al., J. Amer. 
Chem. Soc., vol. 74, Oct. 5, 1952, pp. 
4794-4797. 


Instrumentation and Ex- 
perimental Techniques 
Rocket Test Aid Measures Pressures, 


Week, vol. 57, Nov. 10, 
1952, p.8 

oie of Pressure Gages for Work 
in Ballistics, by S. Raynor, J. Appl. 
Mech., vol. 19, Dec. 1952, pp. 517-520. 

Short-Time Measurements (in Frencli), 
by P. Fayolle, Actes Coll. inter. Mécin. 
II, Publ. sci. tech. min. air, Paris, no. 250, 
1951, pp. 51-64. 

A Chronophotographic Method of Re- 
cording Trajectories of Fast Devices from 
a Distance of Several Kilometers (in 
French), by A. Sestier, Actes Coll. inter. 
Mécan. II, Publ. sci. tech. min. air, Paris, 
no. 250, 1951, pp. 65-72. 

Double-Flash Microsecond Silhouette 
Photography, by H. E. Edgerton, Ri». 
Sci. Instrum., vol. 23, Oct. 1952, pp. 
532-533. 

The Development of a Multi-fla-h 
Camera and Its Application to the 
Study of Liquid Jets, by M. A. Mahrous, 
Brit. J. Appl. Phys., vol. 3, Oct. 1952, 
pp. 329-331. 

Radiofrequency Mass ee for 
Upper Air Research, by J. W. Townsen:|, 
Jr., Rev. Sci. Instrum., vol. 25, Oct. 1952, 
pp. 538-541. 

Mass Spectrometer; New Technique 
for Analysis of Metals, unsigned, Metal 
Indus., vol 80, Aug. 22, 1952, pp. 149-150. 

Six-Ton Parachutes, unsigned, /nteraviv, 
vol. VII, no. 8, 1952, pp. 444-445. 

Instrument Servomechanisms, by W. |). 
Deerhake and A. C. Hall, USN ONR 
Project 24-DC-1, U. 8S. Department of 
Commerce, Office of Technical Services, 
Washington 25, D.C. 265 pp. 

Force-Reflecting Positional Servo- 
mechanism, by C. Goertz and PF. 
Bevilacqua, Nucleonics, vol. 10, Nov. 
1952, pp. 43-45. 

Fundamental Aspects of a Hydraulic- 
Type-Powered Flight-Control System, by 
H. M. DeGroff, Aero. Engng. Rev., vol. 71, 
Nov. 1952, pp. 42-55. 

Fundamentals of General-Purpose Re- 
mote Manipulators, by R. C. Goertz, 
Nucleonics, vol. 10, Nov. 1952, pp. 36-42 

Bearings for Gyroscope Rotors, un- 
signed, Engineer, vol. 194, Nov. 14, 1952, 
p. 651. 

Proc. Inst. Radio Engrs., Transistor 
Issue (50 papers), vol. 40, Nov. 1952. 

High Frequency Heating, by A. E. 
Williams, Engineer, vol. 194, Aug. 15, 
1952, pp. 215-219. 


Terrestrial Flight, Ballis- 
ties, and Vehicle Design 


Windtunnel is a Firing Range, unsigned, 
Aviation Week, vol. 56, Oct. 20, 1952, p. 35. 

Ballistics for Engineers (in Swedish), by 
S. Rydberg, Tekn. Tidskr., vol. 81, Dec. 
1951, pp. 1073-1083. 

Pilotless Aircraft, by R. Goldin, Aero. 
Engng. Rev., vol. 11, Nov. 1952, pp. 22-28. 

The First Guided Missile Was a British 
Invention, by A. M. Low, Flight, vol. 61, 
Oct. 3, 1952, pp. 436-438. 

British Guided Missiles, unsigned, E'n- 
gineering, vol. 173, Aug. 1, 1952, p. 137. 

Guided Missiles Blast Off New Oerlikon 
Arms Program, by Swiss correspondent, 
Aviation Age, vol. 17, Dec. 1952, pp. 25-27 

Missile Guidance System Outlined 
unsigned, Aviation Week, vol. 57, Nov. 24 
1952, pp. 60-61 

Shipping Unit Saves Missile from Shock. 
unsigned, Aviation Week, vol. 57, Nov. 24 
1952, p. 28 


ARS Journal 


by 
19 
Fl 
in 
Se 
tic 
Ki 
19 
De 
(w 
ri\ 
pp 
by 
Ne 
in 
Pa 
Pi 
10 
Ay 
Te 
3 
R: 
by 
Re 
: th 
J. 
; 
Ey 
Bs} 
Sir 
>) 
sig 
by 
Ne 
of 
Co 
me 
of 
nic 
. 43 
an 
pp 
‘ 


ous, 
B52, 


for 
anid, 
952, 


que 
etal 
150. 


vw, 


NR 

of 
ces, 


OV. 


Space Flight 


The Practical ye Ty to Astronautics, 


by J. Humphries, Flight, vol. 61, Oct. 24, 
1952, pp. 528-532. 

l'scape from Earth, by J. G. Vaeth, 
Flying, vol. 52, Dec. 1952, pp. 27-44. 

Ascent from Earth, by L. Grant, J. 
Spuce Flight, vol. 4, Sept. 1952, pp. 1-4. 

A Possible Method of Generating Power 
in Space, unsigned, J. Space Flight, vol. 4, 
Sept. 1952, p. 7. 

\ Preliminary Survey of the Construc- 
tioual Features of Space Station, by H. B. 
Keichum, J. Space Flight, vol. 4, Oct. 
1952, pp. 1-4. 

\on Braun Offers Plan for Station in 
Sp.ce, by Von Braun, Aviation Age, vol. 17 
Dev. 1952, pp. 61-63. 

interorbital Transport Techniques 
(with Special Reference to Solar De- 
rived Power), by H. Preston-Thomas, 
J. Brit. Interplan. Soc., vol. 11, July 1952, 
pp. 173-193. 

lights to the Major Planetary Systems, 
by H. B. Ketchum, J. Space Flight, vol. 4, 
Nov. 1952, pp. 1-3. 


Astrophysics, Aerophysics. 
and Atomic Physics 


Pressures, Densities, and Temperatures 
in the Upper Atmosphere, by The Rocket 
Panel, Harvard College Observatory, 
Phys. Rev., vol. 88, Dec. 1952, pp. 1027- 
10:32. 

Radiation Hazards in High-Altitude 
Aviation, by A. Tobias and E. R. Smith, 
Tech. Data Digest, vol. 17, Nov. 1952, pp. 
3-14. 

Deflection and Diffusion of a Light 
Ruy Passing Through a Boundary Layer, 
by H. W. Liepman, Douglas Aircraft 
Report SM 14397, May 16, 1952, p. 12. 

The Increase or Decrease of Mean-Flow 
nergy in Large-Scale Horizontal Flow in 
the Atmosphere, by G. W. Platzman, 
J. Meteor., vol. 9, Oct. 1952, pp. 347-358. 

High-Altitude Research, by FE. Burgess, 
Engineer, vol. 194, Sept. 19, 1952, pp. 
370-373. 

High-Altitude Research, by E. Burgess, 
Engineer, vol. 194, Sept. 12, 1952, pp. 
338-340. 

Meteorites and Cosmic Rays, by S. F. 
Singer, Nature, vol. 170, Nov. 1, 1952, pp. 
728-729. 

Reactor That Even Looks Safe, un- 
signed, Business Week, Nov. 8, 1952, p. 76. 


General Topics 


Reliability Analysis of Modern Weapons 
by M. M. Munk, Aero. Digest, vol. 63, 
Nov. 1952, pp. 101-102, 104, 106-108. 

Probabilitistic Logics and the Synthesis 
o! Reliable Organisms from Unreliable 
Components, delivered by J. von Neu- 
mann, noted by R. S. Pierce, Calif. Inst. 
of Technology Lectures, Jan. 1952, 51 pp. 

From Newton’s Note Book to Mam- 
moth Research Organizations, by T. von 
Karman, Interavia, vol. VII, 1952, pp. 
431-432. 

The University and Government Spon- 
sored Research, by V. L. Parsegian, Chem. 
and Engng. News, vol. 30, Oct. 27, 1952, 
pp. 4470-4472. 

Records of Research, by J. G. Jackson, 
/. Franklin Inst., vol. 254, Nov. 1952, 
pp. 355-368. 

Through the Curtain to Woomera, by 
}.. C. Shepherd, Aeroplane, vol. 82, Nov. 
7, 1952, pp. 630-633. 
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Recognized as outstanding producers of small parts and 


close tolerance work ... with plenty of Yankee ingenuity 
plus a wide open mind for the best ideas of others. 

In our modern plant we also have the automatic screw 
machines for producing quality parts up to 1%” diameter 
by 8%” long. 

Through constant research and development work 
in the interest of our sinessilines we have maintained our 
leadership as specialists in fine instrument screws for over 
sixty years... with a variety of products which cover 
every department of the screw machine products wenoen oF 


field. HON 


Specialists in Screw Machine Products Since 1893 


ALTHAM 


SCREW COMPANY 


76 Rumford Avenue, Waltham, Massachusetts 
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CURRAN ENGINEERING CO. 


A 


Manufacturer of 


MECHANICAL COMPONENTS 
from 
METALS, CERAMICS, AND PHENOLICS 


A 


Consultants and Specialists of 
ROCKET IGNITER ASSEMBLIES 
and 
LONGITUDINAL SHAPED CHARGE 
CUTTERS 
‘“CENGO" Process for 
HIGH TEMPERATURE-HIGH DIELECTRIC 
INSULATING OF METALLIC ASSEMBLIES 


4423 W. Jefferson Blvd. 
Los Angeles 16 California. 


ponents for 
automatic 


flight 


Accelerometers 
Aeroheads 


G. M. GIANNINI & CO., INC., Pasadena 1, California 


1946 


EX erving the rocket industry since 


“O” RINGS 
MANUFACTURED OF POLYETHYLENE, TEFLON, 
AND KELL-F TO A.N. SPECIFICATIONS 

HOSES 
FLEXIBLE HIGH PRESSURE — BRAIDED STAINLESS 
STEEL POLYETHYLENE LINED WITH A. N. END 
CONNECTIONS 

POLYETHYLENE SHEET, RODS, AND TUBES 
LARGE STOCK AVAILABLE FOR GASKETING AND 
CHEMICAL CARRIAGE 

VALVES 
STAINLESS STEEL CHECK VALVES—HIGH PRESSURE 
— POLYETHYLENE OR TEFLON POPPET 


ef SPECIALTIES COMPANY 


INDUSTRIAL SOUND 


SOUND CONTROL 


INDUSTRIAL 
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—_ 3348 East 14th Street, Los Angeles, California 
CONTROL @ INDUSTRIAL SOUND 9 
4 
rocket 
silenced by 
iSC mufflers’ 
@ In the laboratory, in the test cell, on the 7) 
airfield, hundreds of INDUSTRIAL SOUND «a 
CONTROL installations are subduing the . 
noise, heat and gas velocities generated dur- 
ing testing of the big jets. = 
Whatever your noise problem, ISC’s skilled > 
engineering, design, and installation “know = 
how” — gained through years of practical 
experience — can help you, and is at your 7) 
instant service. re) 
We welcome the challenge of the unusual 
problem. For full information c 
WRITE ISC TODAY 2 
Foreign Licensees: 9 
Cementation (Muffelite) Limited, London 
Les Travaux Souterrains, Paris a 
Industrial Sound Control, Inc. 
ndustrial Sound Control, inc. - 
45 Granby Street, Hartford, Conn. # 
2119 So. Sepulveda Bivd., Los Angeles, Calif. 0 
© 
@indDUSTRIAL SOUND CONTROL 
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AIRCRAFT | 
INSTRUMENTS 


At cousman 


Twenty-three years back, the Kollsman Instru- Aircraft Instruments and Controls * * Min- 
ment Corporation took its first forward steps by _iature AC Motors for Indicating and Remote 
devising a sensitive altimeter twenty times more Control Applications * * Optical Parts and 
responsive than preceding models. This instru- Optical Devices * * Radio Communications 
ment actually made possible the first completely and Navigation Equipment 

blind flight in the history of aviation. While current production is largely devoted to 
Today, this venturesome spirit continues tospark _—_ National Defense, the research facilities of the 
the design, development and manufacture of  Kollsman Instrument Corporation remain avail- 
products of precision and dependability in the _ able to scientists working toward the solution of 
four distinct yet allied fields of: instrumentation and control problems. 


SUBS'DIARY OF 


Standard coit PRODUCTS CO. INC. 


Marcu-Aprin 1953 


- - ----- ------- -- --- ----- 
; 
KOLLSMAN INSTRUMENT CORPORATION 
ELMHURST, NEW YORK 
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LIQUID OXYGEN 
GENERATORS 

An air-transportable 
plant for the separa- 
tion of high-purity oxy- 
gen from the air. 


LIQUEFIED-GAS 

STORAGE CONTAINERS 
Containers for the stor- 
age or transport of 
liquefied gases such as 
liquid hydrogen or oxy- 


a” 


LIQUID OXYGEN 

PUMPS 

inating unit for supply- ] mathematics, 
ing high-pressure oxy- > al economics 
gen gas from low-pres- 

age containers. 


HEAT EXCHANGERS 
A gas-to-gas heat ex- 
changer with excep- 
tional characteristics. 


HELIUM 
REFRIGERATORS 


the prevention of evap- 
oration loss in stored 
liquefied gases. 


AIR COOLERS 


A unit which provides 
cooled, compressed air 
which is free of dirt, 
oil, or entrained water. 


THERMODYNAMICS HEAT TRANSFER REFRIGERATION TO MINUS 456°F GAS LIQUEFACTION 
e VACUUM ENGINEERING + ELECTROMAGNETISM + MECHANICAL DESIGN + VIBRATION 
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chanical Division to provide industry 
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|| 
h our engineers in prototype 
high level of engineering skills. 
= staff is experienced in interpreting the 
| = __ ideas of industry and following 
through with the perfection of 
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